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Summary. 


One  of  the  major  problems  in  obtaining  the  surface 
friction  coefficient  is  that  accurate  velocity  measure¬ 
ments  must  be  obtained  in  the  boundary  layer.  There  are 
several  methods  of  achieving  this  and  it  is  hoped  that 
the  following  will  give  the  reader  a  better  understand¬ 
ing  as  to  the  applications  ana  disadvantages  of  these 
techniques . 


'■One  of  the  more  modern  approaches  for  evaluating  the 
surface  friction  coefficient  is  use  of  the  logarithmic 
inner  la w  which  has  been  used  in,  and  the  results  veri¬ 
fied,  using  the  Prandtl  power  law  technique. 
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velocity  vector. 

friction  velocity  Eq..  (3.14)  • 

u/u«  Eq.(8.31)« 
freestrean  velocity. 

U/u.  Eq. (8.33) • 

mean  velocity. 

average  velocity. 

velocity  vector. 

voltage;  velocity  vector. 

wake  function  Eq.(8.34);  Eq.(8.36). 

variable. 

distance  from  surface. 
yu*/i>  Eq.(8.32). 
variable. 


Greek  symbols. 


0  Clausers  parameter  Eq.(8.28). 

i  boundary  layer  thickness. 

£  displacement  thickness  Eq.(8.4). 

A  constant  in  Eq.(3*l)»  defect  thickness  Eq.(8.29) 

8  angle;  momentum  thickness  Eq.(8.6). 

A  wavelength;  local  skin  friction. 

p-  absolute  viscosity. 

V  kinematic  viscosity.  P’/p 

7T  3.1415327 

7T  Coles  wake  parameter. 

P  density. 

t:  shear  stress. 

63  angular  frequency. 


Dimensionless  groups. 

Of  skin-friction  coefficient. 

Pr  Proude  number  Eq.( 2.1.4). 

Re  Reynolds  number  Eq.(4«3). 

Reynolds  number  Eq.(8.10). 


xlii 


Subscripts 


oo  far  field, 

f  friction, 

o  wall . 

t  turbulent, 

i  incident, 

s  scattered. 

D  Doppler  frequency 
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HtELH-UKARY  CUTLET Z  AND  SURFACE  REPLICA. 


1.1  Introduction. 

The  aim  of  this  project  is  to  compare'  and  critically 
analize  the  methods  of  deriving  the  surface  friction 
coefficient  of  a  replica  of  a  badly  corroded  ships 
plate.  The  particular  replica  is  attached  to  the  bottom 
of  a  20m  flow  channel  which  can  simulate  various  flow 
situations. 

The  surface  friction  coefficient  is  obtained  from  the 
velocity  profiles  and  use  of  analytical  expressions. 
Boundary  layer  hot  film  anemometer  probes  were  used  for 
obtaining  velocity  profiles  and  the  accuracy  and 
reliability  of  such  measurements  were  to  be  checked 
using  a  Laser  Doppler  Anemometer. 

It  is  unfortunate  that  due  to  a  delay  outside  the  author's 
control  that  the  Laser  Doppler  equipment  was  not  obtained 
until  r.arch  1979  which  will  account  for  the  results 
from  the  LDA  being  rather  scanty. 

The  use  of  the  LDA  and  HFA  has  the  backing  of  SRC  and 
also  of  the  U.S.Navy. 

Further  investigation  using  LDA  equipment  is  envisaged 
as  a  research  project  sponsored  by  the  U.S.Navy,  B.S.R.A. 
and  Shell  Research  Co. 
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1.2  Replica  of  Ships  Kull. 

The  replica  of  the  ships  hull  was  supplied  by  Shell 
Research  Co.  It  is  approximately  6.1  metres  long  by 
30.5cm  wide.  Initially  a  polystyrene  negative  was  cast 
on  the  side  of  a  tanker  which  has  been  in  continuous 
service  for  25  years. 

As  is  customary  wita  the  Shell  t&nker  fleet,  the  ship 
has  been  dry-docked  every  three  years  and  the  hull  shot- 
blasted  and  repainted  at  such  regular  intervals.  The 
negative  was  taken  after  dry-docking  but  before  any 
treatment  was  applied  other  than  hosing  down  with  clean 
water  to  remove  slime. 

The  negative  was  then  brought  into  the  laboratory  and 
a  frame  made  the  width  of  the  channel.  The  negative  was 
placed  horizontally  in  the  frame  and  a  cast  made  from 
the  negative  in  an  Epoxy  resin,  care  was  needed  so  that 
the  pouring  of  the  resin  did  not  in  any  way  damage  the 
surface  profile  on  the  negative. 

Once  the  cast  was  dry  the  replica  was  removed  from  the 
mould.  After  the  negative  was  removed,  some  imperfections 
in  the  casting  were  revealed,  mainly  at  the  edges. 

These  were  made  good  with  glass-fibre  resin  after  the 
replica  was  placed  in  the  flow  channel. 
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1«3  Installation  of  cast  replica  in  the  Channel. 

The  cast  replica  is  approximately  25mm  thick  and 

j 

therefore  to  eleminate  any  disturbing  effects  caused 
by  the  abrupt  rise  in  xbe  channel  floor,  an  aluminium 
ramp  was  installed.  Also  to  eliminate  any  effects  from 
the  end  condition  a  similar  ramp  was  attached  at  the 
rear  of  the  model. 

The  model  was  now  sealed  along  its  edges  using  a  plumber 
paste  so  as  to  obviate  any  disturbances  caused  by  water 
passing  between  the  edges  of  the  model  and  the  glass 
sides  of  the  channel  or  the  occasional  liberation  of 
trapped  air  buocles. 

This  method  of  sealing  also  held  the  model  in  place 
on  the  bottom  of  the  flow  channel. 

After  carrying  out  certain  flow  tests  over  the  model  it 
was  found  that  the  initial  method  of  securring  the  model 
to  the  base  of  the  channel  was  totally  unsatisfactory. 
This  was  due  to  the  fact  that  even  at  moderate  flow 
velocities  of  the  order  of  0.5m/s  the  model  would  be 
parted  from  the  base  due  to  the  combined  buoyancy  effects 
and  the  viscous  drag  caused  by  the  flowing  medium. 

Once  the  model  raised  even  slightly  above  the  channel 
base  the  flowing  stream  lifted  the  model  plate  when  the 


dynamic  force  exceeded  its  apparent  weight.  Therefore 
a  more  secure  method  had  to  be  employed. 

The  next  method  used  was  to  wedge  the  model  between  the 
glass  sides  of  the  channel.  Small  wedges  in  hardwood 
v/ere  made,  and  this,  for  a  time  proved  very  successful. 
However,  this  method  after  a  period  of  time  again  proved 
unsatisfactory,  and  it  was  therefore  decided  to  bolt  the 
model  to  the  channel  base.  It  should  be  remembered  that 
uhe  first  two  methods  used  were  done  so  as  to 

a)  not  provide  any  obstructions  to  the  flow  so  as  to 
cause  an  interuption  on  the  surface,  and, 

b)  not  to  damage  the  surface  profile  in  any  way 
whatsoever. 

For  the  third  and  final  method  it  was  decided  to  bolt 
the  model  to  the  base  using  6-^-inch  BSF  countersunk 
bolts  spaced  out  and  placed  in  areas  wnere  the  surface 
was  fairly  smooth.  See  figs.  1.1.  and  1.3. 


Fig.  1.  1. 

Showing  positions  of  holes  relative  to  model. 
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SECTION  (A. A.) 


Showing  method  of  securing  model  to  channel  base." 

The  model  and  base  were  initially  drilled  using  i-inch 
Tapping  drill,  next  the  model,  and  only  the  model  was 
drilled  using  a  £incb  dia.  drill  and  the  hole  in  the 
base  tapped.  The  model  was  now  countersunk  so  as  to 
completely  hide  the  head  of  the  bolt  from  the  surface. 
The  remaining  area  between  the  bolt  head  and  surface 
was  then  filled  using  plumbers  paste.  It  should  be  noted 
tkac  the  fixing  bolts  were  all  placed  away  from  the 
centre  line  as  this  is  the  main  region  in  which  velocity 
measurements  were  to  be  taken. 
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20m  PLOW  CHANNEL  EQUIPMENT. 


2.1  The  Flow  System. 

The  existing  2Cm  flow  channel,  plate  1,  is  30.5cm  wide 
ana  30cn  deep. 


PLATE  1  SHOWING  20e  FLOW  CHANNEL. 


The  channel  can  be  served  using  two  main  pumps,  these 
will  be  referred  to  as  PI  and  F2  respectively.  In  these 
investigations  the  small-capacity  axial-flow  pump,  ?2 
was  never  used  due  to  faulty  bearing,  but  it  is  enough 
to  say  that  it  has  a  much  lower  delivery  than  PI  and 
was  the  original  pump  that  supplied  this  channel  with 
water . 
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PI  however,  is  a  fairly  recent  installation  and  has  the 
advantage  that  it  can  be  throttled  in  order  to  obtain 
very  low  flow  rates;  at  its  rated  discharge  it  will  give 
super-critical  flow  in  she  channel  with  a  suitable  slope. 

The  channel  has  the  capability  of  being  set  at  varying 
slopes  by  means  of  an  electrically  operated  tilting 
apparatus.  This  ensures  that  the  different  type  of  flow 
modes  can  be  obtained. 

Sub-critical  and  super-critical  flow  are  defined  by  the 
Froude  number,  SJ.:  For  values  of  FT  less  than  one,  the  flow 
is  sub-critical,  and  when  FP  is  greater  than  one,  the  flow 
is  super-critical.  The  Froude  number  can  be  defined  as 
the  modulus  of  the  inertia  force  divided  by  the  modulus 
of  the  gravity  force  which  gives  the  non-dimensional 
equation  2.1.1. 

Froude  No.  =  U* 

Ig  (2.1.1) 

However,  in  practice  it  is  often  convenient  to  use  the 
square  root  of  the  ratio  so  as  to  have  the  first  power 
of  the  velocity. 

This  is  quite  permissible:  equality  of:- 

TTS7*1  (2.1.2) 

implies  equality  of:- 
U* 

Ig  (2.1.3) 
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Equation  (2«1.2)  ia  the  more  usual  definition  of  Froude 
number  3J..  Again  this  yields  a  non-dimensional  equation. 


In  long,  open  channels  the  usual  form  of  toe  Froude 
number  is:- 

Fr  "  TgH7'*  (2.1.4) 

Where  U  is  the  velocity  of  the  flow,  g  is  gravitational 
acceleration,  and  h  is  the  mean  depth  of  flow. 

The  mean  depth  of  flow  being  defined  as  the  cross 
sectional  area  of  the  stream  divided  by  the  width  of 
the  channel. 

For  the  purpose  of  this  investigation  the  flow  was 
maintained  below  its  critical  value.  This  is  to  ensure 
that  any  flow  instabilities  are  reduced  to  a  minimum. 
Any  misalignment  of  the  glass  panels  which  make  up  the 
sides  of  the  channel,  or  undulations  in  the  channel 
floor,  would  magnify  any  instabilities  throughout  the 
flow  region  at  the  critical  and  super-critical  stage, 
whereas  at  tranquil  (i.e. sub-critical)  stage,  the 
effects  are  not  as  serious. 

Upon  installation  of  the  glass  panels  in  the  channel, 
care  was  taken  to  ensure  that  the  glass  plates  were 
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perfectly  aligned,  and  the  two  sides  made  parallel  to 
very  close  limits.  Any  gaps  in  the  joints  are  smoothed 
out  using  a  sealing  compound. 
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2.2  The  Channel  Apparatus* 

As  previously  mentioned ,  the  channel  is  approx.  20m 
long  by  30.5cm  wide  and  30cm  deep;  therefore  sufficien¬ 
tly  long  working  plate  models  can  be  accommodated . 
However,  allowances  have  got  to  be  made  for  the  inlet 
settling  length  and  exit  length  respectively,  bearing 
in  mind  the  unsuitable  disturbances  propagating  in  the 
deturbulating  tank  at  inlet  and  the  action  of  a  weir 
plate  at  the  exit. 

The  channel  is  supplied  with  clean  water  from  a  pit 
whose  capacity  is  approx.  7* 5m  by  means  of  two  pumps 
electrically  driven,  PI  and  P2  arranged  in  parallel. 

P2  is  an  axial  flow  extraction  pump,  and  the  pump  head 
is  submerged  in  the  water  pit.  The  pump  impeller  is 
some  4ft.  deep  and  is  connected  via  a  drive  shaft 
housed  in  rubber  water  cooled  bearings.  This  pump  is 
capable  of  delivering  20  l/s  at  2900  rev/min. 

However,  due  to  faulty  bearings,  this  pump  was  not  used 
in  the  following  tests. 

P2,  which  is  the  latest  addition  to  the  equipment,  so 
that  the  flow  in  the  channel  can  be  operated  in  either 
the  tranquil,  critical  or  super-critical  flow  mode. 

It  is  a  rotor  dynamic  pump  of  the  radial  flow  type. 
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This  pump  has  been  derated  by  reducing  the  impeller 
diameter  and  the  reduced  delivery  rate  is  some  60  l/s 
at  a  rotational  speed  of  $60  rev/min. 

The  pump  delivers  the  flow  to  the  channel  through  a 
150mm  bore  header  interconnecting  pipe  system,  and 
manually  operated  flow  control  valves.  The  flow  circuit 
incorporates  a  Dali  tube  meter  with  its  pressure 
tappings  connected  to  mercury-under-water  manometers. 

Two  manometers  are  connected  in  parallel:  a  vertical 
one  for  large  flow  rates  and  an  inclined  one  for  accurate 
measurement  of  low  flow  rate. 

Water  is  delivered  to  a  deturbulating  tank,  through  a 
submerged  inverted  T-piece  whose  ends  are  blocked  off; 
the  water  has  to  pass  through  a  multiplicity  of  small 
diameter  orifices  drilled  around  the  bottom  leg  of  the 
T-piece.  It  then  passes  through  a  honey-comb  type  flow 
straigbtener  laid  horizontally  over  the  supply  section 
in  order  to  minimise  swirl  and  damp  out  large  scale 
disturbances  produced  by  the  pipe  work  system. 

The  shape  of  the  deturbulatory  tank  is  such  as  to 
accelerate  the  flow  before  entering  the  channel. 

As  the  flow  passes  through  the  deturbulating  tank  into 
the  channel  it  passes  through  another  honey-comb  flow 
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straightener  whose  purpose  is  to  ensure  flow  approaching 
that  of  uniform  potential  at  the  channel  inlet. 

Without  the  honey-comb  sections,  the  swirl  produced  by 
the  combination  of  bends  in  the  pipe-line  would  be 
conveyed  to  the  channel  inlet  and  would  persist  for 
some  distance  downstream,  thus  causing  non-uniform, 
grossly-disturbed  flow  in  the  working  section  of  the 
channel.  This  situation  would  be  unacceptable  for  taking 
meaningful  velocity  measurements.  As  the  water  passes 
around  the  system  it  must  be  ensured  that  there  is  an 
adequate  reserve  in  the  supply  pit;  if  not,  air  en¬ 
trainment  and  possibly  cavitation  may  occur  in  the 
pump  which  would  mar  the  experimentation. 

The  channel  is  piveted  at  the  entry  section  and  has 
two  jacking  points  connected  to  an  automated  system, 
therefore  the  channel  slope  can  be  adjusted  in  order 
to  give  varying  conditions  in  conjunction  with  the 
outlet  weir  which  can  be  set  to  a  suitable  height  to 
hack  up  the  flow  as  required. 

It  should  be  noted  that  if  the  channel  exit  is  raised 
approximately  two  inches  then  the  weight  of  water  iD 
the  system  i.e.  channel  and  deturbulating  tank,  over¬ 
loads  the  electric  drive  actuates  the  tilting  jacks. 

The  channel  would  not  then  be  able  to  be  lowered. 
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Excessive  overload  could  cause  serious  damage  to  the 


torque  motor.  To  relieve  this  condition  the  channel 
should  preferrably  be  emptied  and  water  released  from 
the  detur'oulating  tank  through  a  drain  valve  at  the 
base,  before  gross  adjustments  are  made  to  the  slope 
of  the  channel. 


2.3  Filtration  System. 

When  using  this  system  it  is  important  that  the  water 
is  as  clean  and  clear  as  possible.  To  achieve  this, 
filtration  is  carried  out  on  a  by-pass  principle;  this 
entails  a  separate  system  whereby  water  is  pumped,  at 
a  constant  rate,  from  the  pit  through  a  small  centri¬ 
fugal  pump,  via  a  hose  into  a  large  cylindrical 
container.  The  water  discharges  by  gravity  back  into 
the  pit  through  a  series  of  filters  housed  one  on  top 
of  each  other. 

* 

Filtration  is  vitally  important  since  loosened  rust 
and  flakes  of  old  paint  from  the  pump,  and  accumulated 
grime  in  the  pipe-line  system  cause  discolouration  of 
the  water  and  the  carry  over  of  all  sorts  of  debris; 
the  resulting  situation  being  such  that  the  accuracy 
and  reliability  of  measurements  could  not  be  depended 
on. 

Also  with  equipment  as  delicate  as  the  Boundary  layer 
probes,  serious  damage  can  be  caused;  this  however 
will  be  discussed  later. 
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2.4  The  Pall  Flow  Tube. 

Situated  approximately  midway  between  the  pump  PI  and 
the  deturbulating  tank  is  a  Dali  Plow  Tube  measuring 
device.  This  unit  is  attached  via  two  bolted  flanges 
and  connected  in  the  main  straight  section  of  the 
interconecting  pipe-line. 

The  Dali  flow  tube  is  a  modified  venture  type  meter 
which  gives  a  measure  of  the  water  flow  rate  through 
the  channel  at  any  particular  set  of  flow  conditions. 

The  flow  meter  has  two  water-on-mercury  manometers 
connected;  one  being  vertical  for  large  flow  rates 
and  the  other  inclined  for  accurate  measure  of  low 
flow  rates.  Due  to  the  working  range  used  in  these 
investigations,  the  inclined  manometer  was  never  used. 

Once  the  pump  PI  was  running  and  water  was  flowing 
through  the  channel,  before  the  manometer  reading 
could  he  relied  upon,  the  system  required  bleeding. 

The  procedure  here  is  to  firstly  throttle  the  flow  using 
the  diaphragm  valve  at  the  deturbulating  tank  inlet 
in  order  to  create  a  back  pressure  in  the  pipe-line. 

This  is  necessary  because  to  ensure  that  air  is  not 
trapped  in  the  Dali  tube  the  two  vent  cocks  must  be 
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opened.  These  vent  cocks  are  on  the  high  pressure  and 
low  pressure  side  of  the  unit  respectively,  and  if  the 
flow  is  not  throttled,  then  air  could  be  drawn  in  through 
the  low  pressure  vent.  Once  all  the  air  was  exhailed 
through  the  vent  cocks,  and  water  was  coming  out,  these 
could  be  closed  and  the  manometer  bled  to  ensure  that 
the  system  was  working  correctly. 

The  diaphragm  valve  could  now  be  re-opened. 

It  is  important  to  note  that  the  diaphragm  valve  should 
always  be  fully  opened  when  the  system  is  started  up, 
this  is  to  ensure  that  the  pressure  exerted  from  the 
flow  as  it  initially  passes  through  the  pipe  system 
does  not  damage  the  diaphragm  plate. 

The  Dali  flow  tubes  main  function  is  that  it  can  be 
used  to  record  the  flow  rate  accurately  for  each 
particular  set  of  measurements  taken.  This  then  makes 
it  possible  to  reproduce  the  flow  condition  with  the 
variable  height  weir  plate  and  the  channel  slope, 
whereby  repeat  procedures  can  be  carried  out  with 


accuracy. 


2.5  Sunning  up  of  System. 

From  the  preceeding  passages  it  is  evident  that  a  range 
of  varying  conditions  can  he  obtained  by: 

a)  Control  of  valve  at  pump  exit. 

b)  Channel  slope. 

c)  Introduction  of  a  weir  plate  at  the  exit  of  the 
channel. 

d)  Adjustment  of  diaphragm  valve  at  entry  to 
deturbulating  tank. 

Therefore  not  only  can  the  conditions  be  controlled 
over  a  significant  range,  but  also  by  using  the  Dali 
Tube  flow  meter  a  given  set  of  flow  conditions  can  be 
repeated  thus  ensuring  that  results  can  be  checked 
and  verified. 
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EQUIFJIEL'IT  FOR  MEASUREMENT  CF  VELOCITY  PROFILES 


3.1  Measurement  of  Velocity  Profiles. 

Before  any  calculations  can  be  carried  out,  it  is 
necessary  to  obtain  velocity  profiles  of  the  fluid  as 
it  passes  over  the  replica  of  the  ship's  hull  surface. 

There  are  four  possible  methods  for  obtaining  velocity 
profiles  in  the  present  flow  system,  and  these  are: 

1)  Pilot-static  traverse. 

2)  Propellor  type  current  meter  (OTT  meter) 

3)  The  Constant  Temperature  Anemometer.  (C.T.A.) 

4)  The  Laser  Doppler  Anemometer.  (l.D.A.) 
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3«1«1  Pilot-Static  Traverse. 

The  pilot-static  tube  is  connected  to  a  water  manometer 
which  is  open  to  the  atmosphere.  The  difference  in  the 
static  and  stagnation  pressures  being  measured  in  terms 
of  height  difference  in  the  manometer  reading. 

The  pilot-static  tube  can  be  connected  to  a  carriage 
arrangement  which  allows  the  unit  to  move  along  the 
length  of  the  channel,  but  also  gives  movement  across 
the  channel.  The  height  of  the  pilot-tube  can  also  be 
adjusted  manually.  Therefore  this  system  allows  pilot- 
static  readings  to  be  taken  in  three  directions. 

This  particrilar  piece  of  apparatus  has  several  dis¬ 
advantages,  and  these  are: 

a)  Due  to  measurements  being  taken  in  a  liquid 
medium,  the  time  required  for  the  manometer  to 
settle  to  an  equilibrium  position  is  such  that 
it  renders  the  equipment  far  to  inefficient. 
However,  there  will  soon  be  a  pressure  transducer 
available  in  the  department  which,  when  connected 
to  the  pilot-static  tube  system,  will  give 
instantaneous  readings. 

At  present  the  method  employed  to  measure  the 
manometer  height  is  by  means  of  a  depth  gauge 
or  depth  micrometer. 
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b)  As  previously  mentioned  the  flow  system  incor¬ 
porates  a  filtration  unit  in  order  to  ensure 
the  water  is  clean.  This  filtration  unit  however, 
will  only  filter  the  larger  size  debris  and 
also  the  filtration  is  such  that  we  cannot 
ensure  that  100$  of  the  water  passes  through 
the  unit.  Therefore  the  water  still  has  impurities 
and  rust  deposits  floating  about,  all  of  which 
tend  to  block  up  the  pilot-static  system.  Thus, 
before  this  unit  can  be  used  successfully,  a 
means  of  purging  the  pilot  tube  at  regular  in¬ 
tervals  must  be  incorporated. 

By  the  introduction  of  the  Pressure  transducer  and  a 
purging  system,  the  pilot-static  tube  would  be  a  reliable 
method  for  velocity  measurement. 


3.1.2  laboratory  Propeller  Type  Current  I'eter. 

The  propeller  cype  current  meter  incorporates  a  propeller 
type  unit  which  comprises  the  meter  body  and  a  propeller 
attached  to  the  end. 


PLATE  2  SHOWING-  PH0PE110H  TYPE  CUEBEKT  L'ETEH . 


The  propeller  is  connected  to  a  spindle  which  is  sit¬ 
uated  in  two  ball  bearings.  Attached  to  the  spindle 
are  two  contact  springs,  these  in  turn  fora  the  elec¬ 
trical  connection  to  an  electric  contact  unit.  Proa 
the  electrical  contact  unit  a  connecting  cable  is 
attached  to  an  electrical  counter  which  registers  each 
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revolution  of  the  propeller.  The  electrical  counter 
unit  is  powered  bp  2  HP2  batteries  connected  in 
series  giving  an  electrical  potential  of  3  volts. 

The  maximum  voltage  the  unit  should  be  subjected  to  is 
6.  ’.Yhen  the  unit  is  enersed  in  a  flowing  medium,  the 
propeller  revolves.  The  equation  from  which  the  mean 
flow  velocity  can  be  calculated  is: 

v  =  Kn  +  a  (3.1.) 

Where  v  =  mean  velocity, 

and  n  =  number  of  revolutions  of  the  propeller/sec 
K  and  A  are  constants  depending  upon  propeller  used. 

The-  propellers  vary  in  size,  i.e.  diameter,  helix  angle 
and  the  correct  propeller  must  be  chosen  to  cover  the 
particular  velocity  range  to  be  used. 

In  this  particular  project,  the  propeller  used  was 
No.  135G9-2  type  Cl,  and  from  the  calibration  charts 
for  the  unit,  i.e.  Calibration  No.  2592  date  of 
calibration  12.10.77* 

The  equation  from  which  the  velocity  can  be  calculated 
directly  is: 

v  =  0.1021n  +  0.0206  (3.2.) 

for  the  range:  1 .677  <  n <  11.827. 
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The  method  employed  in  these  measurements  was  to  time 
the  number  of  revolutions  of  the  propeller  using  a 
stop  watch,  then  by  choosing  the  correct  range  of 
values  for  n  the  velocity  of  the  fluid  can  be  calcul¬ 
ated  directly  using  the  equation  above,  (3*2). 

The  period  of  test  was  60  seconds. 

This  unit  was  mainly  used  to  establish  that  there  is 
a  central  working  region  in  the  channel.  This  can  be 
seen  from  fig.  Di  in  appendix  D 

The  propeller  type  meter,  however,  cannot  be  used  with 
accuracy  for  velocity  measurements  near  to  a  boundary 
because  its  size  and  operating  characteristics  interupt 
the  flow  within  that  boundary.  Therefore  this  unit  was 
mainly  used  so  as  to  establish  the  working  region 
across  the  channel  and  also  the  local  velocity  in 
the  freestream. 


t 
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3.1.3  Constant  Temperature  Anemometer. 

The  constant  temperature  anemometer  used  in  the  -follow¬ 
ing  tests  is  of  the  Disa  Boundary  layer  hot  film  type. 

This  anemometer  consists  of  a  thin  electrically  heated 
Nickel  film  fused  to  a  gold  plated  quartz  fibre  for 
support.  The  fibre  plus  film  has  a  thin  insulating  co&t 
of  quartz,  and  is  held  between  two  arms  of  the  probe 
supports . 

The  probe  traversing  equipment  consists  of  a  guide 
tube  system  with  a  probe  connection  support  attached 
to  a  mounting  tube.  The  probe  connection  support  is 
squeezed  into  position  in  the  mounting  tube  and  the 
cable  from  the  probe  connection  runs  through  the  mount¬ 
ing  and  guide  tube  system.  The  cable  terminates  at 
a  B.N.C.  plug  at  the  top  of  the  guide  tube. 

Prom  the  E.N.C.  plug  is  connected  the  probe  cable 
which  is  attached  to  the  probe  B.N.C.  socket  on  the 
Disa  55KLO. 

Situated  above  the  551-110  is  a  551.101  main  unit.  This 
unit  consists  of  a  meter  for  measurement  cf  probe 
resistance;  a  decade  resistance  with  a  maximum  value 
of  99»99  Ohms  and  various  other  facilities,  i.e.  probe 
type  selector,  square  wave  frequency,  HP  filter  and 
gain  control. 
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The  55^10  unis  is  now  connected  to  a  55D15  linearizer 
so  that  the  signal  outputs  from  the  probe  when  operating 
can  be  linearized  to  read  VotfT  =  Kv  with  good 
approximation,  where  v  =  measured  flow  velocity. 

The  hot  film  probe  forms  one  arm  of  a  Wheatstone  bridge, 
the  other  arm  consisting  of  a  non  inductive  variable 
resistance  box.  The  system  operates  whereby  if  the 
probe  resistance  is  less  than  the  resistance  of  the 
adjustable  bridge  arm,  then  an  error  voltage  will  be 
produced  across  the  amplifier  inputs.  Therefore,  this 
error  voltage,  wiien  amplified  causes  the  probe  to  heat 
up  until  its  resistance  is  equal  to  that  as  previously 
set  on  the  variable  resistance  box. 

The  bridge  is  now  in  balance  again. 

By  placing  the  probe  in  a  flowing  medium,  and  switching 
to  operate,  the  flow  will  try  to  cool  the  probe.  This 
results  in  a  small  charge  in  the  probe  resistance  which, 
as  described  previously,  causes  an  error  voltage  across 
the  probe. 

A  servo-amplifier  connected  inside  the  equipment  feeds 
back  a  signal  causing  the  bridge  current  to  increase, 
and  the  probe  temperature  to  return  to  its  original 
value.  The  voltage  which  is  fed  to  the  bridge  to 


maintain  the  temperature  of  the  probe  can  now  be  related 
to  the  fluid  velocity  by  a  suitable  calibration.  The 
calibration  is  carried  out  using  a  calibration  tube 
which  is  attached  to  the  side  of  the  channel  and  will 
be  discussed  later. 
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3.1.4  Laser  Doppler  Anemometer. 

The  fourth  and  fin&l  method  for  measurement  of  flow 
velocity  is  by  means  of  using  the  L.D.A.  (Laser  L'oppler 
Anemometer . ) 


The  main  principle  of  operation  is  that  a  coherent 
source  of  laser  light  is  split  in  a  beam  splitter  to 
form  two  seams  which  cross  some  distance  from  the 
optical  unit.  This  distance  depends  on  the  focal  length 
of  the  lens  system  chosen.  As  the  two  beams  intersect, 
the  point  consists  of  alternate  light  and  dark  fringes. 


SECTION  THROUGH  MEASURING 

See  fig.  3-1.  VOLUME  SHOWING  FRINGE 

SPACING 


REAR 

REFLECTING 

MIRROR 


FIG.  3.1. 


SCHEMATIC  DIAGRAM  OF  OPTICAL  UNIT 
AND  FLOW  DIRECTION  ADAPTER. 


Therefore  any  particle  which  passes  through  the  fringe 
system  will  omit  light  pulses  at  a  frequency  which  is 
dependant  upon  the  particle  velocity.  This  frequency 
is  known  as  the  Doppler  frequency. 

The  scattered  light  is  now  detected  on  a  Photomultiplier 
and  the  information  relayed  to  the  input  of  a  Preamplifier, 
and  then  to  a  Frequency  Tracker.  The  Frequency  Tracker 
is  an  analoy  instrument  and  'tracks'  the  Doppler  frequency. 
This  is  achieved  "by  means  of  a  voltage  controlled 
oscillator  in  the  feed  back  loop.  The  controlling 
voltage  is  directly  proportional  to  the  Doppler  frequency 
and  hence  also  the  flow  velocity. 

The  Doppler  frequency  can  now  be  obtained  from  the 
analogue  meter  or  by  connecting  a  digital  voltmeter 
where  0-10V  corresponds  to  the  full  scale  meter  reading 
in  all  frequency  ranges. 

The  L.D.A.  will  be  explained  in  further  detail  later. 
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3.2  Summary  of  Velocity  Profile  Ilea  sure  cent  Ev^uipnent. 
From  the  preceding  sections,  3.1.1  to  3.1.4  it  is 
obvious  that  there  are  various  methods  available  for 
measurement  of  the  velocity  profile.  However,  in  order 
to  obtain  accurate  velocity  measurements  from  the 
surface  in  question,  i.e.  the  replica  of  the  ship's 
hull,  to  the  mainstream,  it  is  necessary  to  record 
velocities  in  the  boundary  layer.  Therefore  it  is  vital 
that  the  equipment  to  be  used  does  not  in  any  way 
disturb  the  flow. 

Any  disturbance  of  the  flow  regime  will  give  rise  to 
inaccuracies  and  irregularities  in  the  values  recorded. 
Therefore  before  any  measurements  can  be  taken,  the 
correct  equipment  must  be  used  and  for  these  investi¬ 
gations  the  boundary  layer  hot  film  probes  are  small 
enough  to  cause  minimum  flow  disturbance. 

The  pitot  tube  is  adequate  for  fairly  close  boundary 
layer  measurements,  but  the  physical  size  of  the  unit 
causes  disturbance  of  the  flow  oetween  the  pitot  tube 
and  the  surface,  therefore  for  close  boundary  layer 
work  this  type  of  equipment  is  not  recommended. 

The  propeller  type  current  meter,  due  tc  its  physical 
size,  restricts  its  use  to  anything  other  than  free- 
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stream  velocity  measurements. 


The  boundary  layer  hot  film  probe  is  an  ideal  piece  of 
equipment  for  velocity  measurements  close  to  the  surface 
4  due  to  the  physical  size  of  the  sensing  element; 
therefore  because  of  the  advantages  of  the  hot  film 
probe  over  the  other  type  of  equipment,  makes  it 
ideal  for  velocity  measurements  and  so  was  the  equipment 
used  in  the  investigations  carried  out. 

The  l.D.A.  which  was  to  be  available  December  -  January 
1978  -  1979,  did  not  arrive  until  March  1979*  This 
equipment  would  have  been  the  most  accurate  for 
velocity  measurements  because  the  two  beams  of  Laser 
light  passing  through  the  flowing  medium  do  not  in  any 
way  disturb  the  flow.  Therefore,  had  the  equipment 
been  available  earlier,  it  would  then  have  been 
possible  to  check  the  results  obtained  from  the  C.T.A. 
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4.1  Constant  Temperature  Anemometer . 

The  constant  temperature  anemometer  was  the  main  item 
of  equipment  used  for  velocity  measurements,  however, 
before  any  measurements  could  be  taken  the  probes 
working  resistance  had  to  be  determined. 

4.2  Initial  Setting  of  Controls. 

Before  the  equipment  is  switched  on  ,  and  the  probes 
working  resistance  calculated,  the  following  procedure 
should  be  followed. 

a)  Check  probe  cable  length  is  as  indicated  in  the 
probe  cable  length  window,  i.e.  5,  20  or  100  metres. 

b)  Set  Square  Wave  to  off. 

c)  Set  HP  Filter  to  1. 

d)  Set  Volts  to  1. 

e)  Set  Function  switch  to  'Stand  By'. 

f)  Set  Gain  to  1. 

g)  Set  Decade  Resistance  00.00. 

h)  Set  Probe  type  switch:  depending  on  type  of  probe 
to  be  used,  i.e.  film  or  wire. 
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When  the  power  has  been  applied,  the  pilot  lamp  on  the 
main  unit  will  shov/  light.  The  meter  needle  will  firstly 
deflect  to  the  left  to  below  zero  and  after  about  90 
seconds  will  move  slowly  to  the  right,  indicating  that 
the  various  internal  circuits  are  ready  for  operation. 

Before  beginning  any  measurements,  the  instrument  snould 
be  allowed  to  warm  up  for  15  minutes. 

The  probe  cable  can  now  be  plugged  into  the  probe  socket 
and  the  cable  and  support  terminated,  (short  circuited) 
using  a  shorting  probe.  By  turning  the  Function  switch 
to  the  Res.  Keas  position  this  causes  a  change  in 
position  of  the  meter  needle  which  is  then  set  to  cover 
the  red  mid-scale  mark  by  adjustment  of  the  Zero  Ohms 
potentiometer . 

The  functions  switch  is  now  returned  to  ’Stand  By'  and 
the  shorting  probe  replaced  by  the  probe  to  be  used  for 
measurements . 

Plate  3  shows  the  boundary  layer  probe  in  position  in 
the  probe  connection  support. 


PLATE  3  SEC'TIN G  BOUNDARY  LAYER  PRC3E 


Now  the  probe  ia  placed  in  the  flowing  medium  and  the 
function  switch  set  to  Bes.  I’eas .  (It  is  important  to 
note  here  that  if  the  function  switch  is  adjusted 
without  the  probe  being  immersed,  then  serious  damage 
to  the  probe  or  insulating  quartz  coating  may  result, 
hence  destroying  the  probe.)  Now  by  adjustment  of  the 
decade  resistance  box,  the  bridge  circuit  can  be 
balanced.  The  resistance  reading  on  the  decade 
resistance  box  corresponds  to  the  cold  resistance  of  the 
probe. 

The  decade  resistance  box  is  now  reduced  by  an  amount 
corresponding  to  the  probe-lead  resistance,  this, 
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however,  causes  the  meter  reading  to  change,  which 
should  then  be  brought  back  xo  cover  the  red  mid-scale 
mark  by  adjustment  to  the  Zero  Ohms  potentiometer. 

The  function  switch  should  now  be  re-set  to  'Stand  By'  . 

The  constant-temperature  anemometer  employs  a  sensor 
which  is  kept  at  a  pre-selected  excess  temperature  by 
means  of  an  electronic  control  circuit.  The  excess 
temperature  is  defined  by  the  so-called  overheating 
ratio:  'a ' . 

a  s=  R  -  Ro 

Ko  (4*1.) 

R  s  sensor  hot  resistance, 

Ro  =  sensor  cold  resistance  at  ambient  temperature. 

Sensor  hot  resistance  may  therefore  be  calculated  as 

B  *  Ro  (1  +  a)  (4.2.) 

In  this  investigation,  the  overheating  ratio  was  calculated 
to  be  1.1 , therefore  by  multiplying  the  resistance 
reading  on  the  decade  box  by  this  factor,  the  decade  box 

4 

can  be  re-set  so  that  the  overheating  ratio  is  provided 
automatically. 

Now,  by  turning  the  Function  switch  to  the  'operate’ 
position,  this  closes  the  servo  feed  back  loop  and  the 
anemometer  is  ready  for  operation. 


4.4  Balancing;  The  Bridge. 

Before  carrying  out  any  measurements  using  the  anemometer, 
it  must  be  checked  to  ensure  that  the  instrument  is  not 
oscillating  due  to  inbalance  of  the  bridge. 

The  frequency  behaviour  of  the  anemometer  can  be  determined 
by  producing  a  sudden  change  in  the  velocity  of  the  flow 
acting  on  the  probe.  However,  since  it  is  difficult  to 
produce  accurately  defined  sudden  velocity  changes, 
the  ideal  way  is  by  feeding  a  square-wave  signal  into 
the  bridge. 

By  connecting  an  oscilloscope  to  one  of  the  out  sockets 
on  the  55K01  and  with  the  sensitivity  set  at  approx. 

0.1  V/cm  it  is  possible  to  observe  the  impulse  response 
without  superimposed  oscillation.  However  it  is  necessary 
during  this  test  to  expose  the  probe  to  a  constant 
flow  whose  magnitude  should  at  least  equal  the  maximum 
velocity  occuring  during  the  measurement.  Therefore 
this  test  was  carried  out  in  the  calibration  tube  with 
the  globe  valve  at  inlet  set  for  maximum  flow,  approx 
lm/a . 

The  switch  on  -she  .Square  Wave  Generator  can  be  set  to  a 
desired  frequency  depending  on  the  oscilloscope  pattern. 
Both  test  frequency  and  oscilloscope  frequency  being 
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selected  to  give  an  easily  readable  pattern. 

It  is  necessary  to  have  a  good,  clear  trace  so  that  the 
changes  caused  by  operating  the  controls  give  a  clear 
indication  of  what  is  happening. 

The  Gain  switch  should  now  be  increased  until  the  res¬ 
ponse  signal  on  the  oscilloscope  shoves  damped  oscillations 
on  the  base  line.  By  increasing  the  HP  filter,  the 
oscillations  should  reduce.  Further  adjustment  should 
promote  oscillations  at  the  skirts  of  the  signal. 

The  oscillations  observed  can  however  be  compensated 
for  by  adjustment  of  the  Cable  Compensation  settings 
Q  and  L.  This  adjustment  should  firstly  start  with  I 
and  as  compensation  proceeds  it  may  be  necessary  to 
choose  different  settings  of  Q. 

The  Gain  and  HF  Filter  switches  can  now  be  increased 
and  the  compensation  procedure  repeated. 

The  above  should  be  carried  out  until  no  further 
compensation  can  be  accomplished,  The  Gain  and  HF  filter 
can  now  be  backed  off  until  the  oscillations  disappear. 

The  anemometer  is  now  set  for  optimum  frequency  response. 
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4.5  Linearity. 

In  order  to  obtain  results  that  folio w  a  linear  rela¬ 
tionship,  the  55MLO  is  connected  to  a  55D15  linearizer. 
By  suitable  adjustment  of  the  gain  and  suppression 
controls,  the  unit  can  be  adjusted  to  give  full  scale 
deflection  at  maximum  flow  velocity. 

The  full  scale  deflection  represents  3  Volts  LG  out¬ 
put  which  can  be  transferred  to  a  digital  voltmeter 
if  required. 


37 


4.6  Calibration  of  Boundary  Layer  Probe. 

Before  any  velocity  measurements  in  the  cnannel  can  De 
recorded,  it  is  necessary  to  calibrate  the  probe  over 
the  velocity  range  to  be  encountered.  This  is  done  by 
means  of  a  calibration  rig,  as  shown  in  Plate  4. 


ELATE  4  SHOWING  CALIBRATION  RIG 

The  rig  consists  of  a  2.15cm  dia.  calibration  tube, 

4.6m  long.  The  length  to  diameter  ratio  is  213,  there¬ 
fore  fully  developed  turbulent  flow  will  be  achieved  at 
the  pipe  exit.  (Por  fully  developed  turbulent  flow,  the 
length  to  diameter  ratio  must  be  greater  than  50. ) ^ 
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The  calibration  tube  is  connected  to  a  constant  head 
reservoir  via  a  globe  valve  and  hose  ays tec ,  the 
reservoir  being  supplied  with  water  pumped  from  the 
channel.  This  ensures  the  same  water  is  being  used 
during  both  calibration  and  testing. 

The  calibration  tube  is  attached  to  the  side  of  the 
channel  above  the  level  of  the  water,  and  the  probe 
inserted  at  the  exit  by  means  of  the  existing  carriage 
support . 

Calibration  is  carried  out  by  adjusting  the  globe  valve 
to  give  varying  flow  rates,  and  at  each  flow  condition 
3  separate  sets  of  voltage  readings  are  taken. 

Due  to  turbulance  and  flow  instabilities  the  recorded 
voltage  at  any  one  setting  is  not  static  but  fluctuates. 

The  degree  of  the  voltage  fluctuation  depending  on  zhe 
level  of  turbulance.  This  fluctuation  however,  renders 
the  analog  meter  on  the  linearizer  inadequate  for 
accurate  voltage  measurements.  Therefore,  the  linearizer 
is  connected  to  a  Solartron  Digital  Voltmeter  type  1M.  1420.2 
which  is  connected  into  a  paper  tape  data  logging  system. 

By  switching  the  Voltmeter  to  the  record  setting,  this 
automatically  samples  the  voltage  readings  at  intervals 
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of  0.33  seconds.  During  a  test  period  of  10  seconds, 
the  total  number  of  voltage  readings  recorded  is  30. 


The  voltage  results  obtained  can  now  be  summed  and  the 
* 

mean  value  calculated.  This  calculated  value  is  taken 
as  the  true  voltage  for  that  particular  valve  setting. 
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4.7  Calculation  of  Centre-Line  Velocity. 

Once  the  experimentally  recorded  voltage  and  flow  rate 
are  obtained,  it  is  necessary  to  calculate  the  actual 
velocity  of  flow  to  which  the  probe  is  subjected. 


The  probe  was  situated  along  the  centre-line  of  the 
calibration  tube,  therefore  the  centre-line  velocity 
must  be  calculated,  (this  being  the  actual  velocity  to 
which  the  voltage  refers.)  by  first  evaluating  Reynolds 
number.  Now  by  transferring  the  Re.  number  to  fig.  4.1 
or  4*2,  whereby  the  ratio  of  Average  velocity  /  Centre¬ 
line  velocity  can  be  obtained,  so,  the  Centre-line 
velocity  can  be  calculated  directly. 


Reynolds  No.  (Re)  = 


pva-1 

A 

Va.1 

V 


where  1 
and  $ 


dia.of  calibration  tube, 

A 

P 


(4.3) 

(4.4) 


The  Kinematic  Viscosity  (V>  )  varies  with  temperature, 
and  therefore  using  fig.  4.3  or  4.4  (Ref.  I)  this  may 
be  obtained  at  the  required  temperature. 


The  average  velocity,  va  is  calculated  by  dividing  the 
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Kinematic  viscosity,  mma/5 


PIG.  4.3  Kinematic  Viscosity  of  Fluids. 
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average  flow  rate  by  the  area  of  the  calibration  tube. 
During  calibration  of  the  probe,  the  fluid  mode  is 
always  turbulent  with  Reynolds  numbers  greater  than 
2300.  Therefore  due  to  the  shape  of  the  velocity 
profile  under  turbulent  conditions,  it  has  been 
assumed  that  precise  location  of  the  probe  along  the 
centre-line  of  the  calibration  tube  is  not  critical, 
and  that  for  a  small  distance  either  side  of  the  centre¬ 
line,  the  velocity  difference  will  be  approaching  that 
of  the  centre-line  velocity.  However,  location  of  the 
probe  was  carried  out  as  accurately  as  possible. 

It  should  be  noted  that  for  each  valve  setting,  the 
probe  was  located  in  exactly  the  same  position  to 
enable  uniformity  in  experimental  procedure. 

The  calculated  values  of  centre-line  velocity  and 
average  voltage  can  be  seen  in  tables  A3  and  A4 
Appendix  A. 


4.8  Probe  Calibration  Graphs . 

To  obtain  a  relationship  bet'ween  the  voltage  and 
centre-line  velocity,  it  is  necessary  to  graph  the 
results:  Tables  A3  and  A4  Appendix  A  which  refer  to 
probe  1  and  2  respectively. 

These  graphs  indicate  that  the  plotted  points  exhibit 
a  linear  relationship  governed  by  the  law  y  =  bx  +  a, 
where  'y*  and  'x'  are  the  variables,  and  'b'  and  ’a1 
are  the  slope  and  intercept  respectively. 

In  order  to  obtain  the  best  straight  line  to  be  drawn 
through  the  experiaentally  obtained  data,  it  is  necessary 
to  i;se  the  'Method  of  least  Squares'  analysis,  initially 
developed  by  Gauss. 

For  this  analysis  there  are  3  basic  equations  which  are: 


n  a  +  b£  X  = 

£  Y 

(4.5) 

a  £  X  +  b^X*  = 

£XY 

(4.6) 

n£XY  -  'SXSY 

-  -4  1 

(4.7) 

J (n*X*-  (  *X)x  KdSY1-  (SY?) 


where  (4.7)  indicates  whether  or  not  statistically  the 
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line  gives  a  good  fit  to  the  plotted  experimental 
data . 


In  the  previous  equations,  'a'  and  'b'  are  constants, 
and  'X'  and  'Y'  refer  to  the  experimentally  obtained 
results . 

Por  the  calculations  using  the  afore  said  equations, 
'X'  is  usually  designated  as  the  value  from  which  '  Y' 
relies  upon,  therefore,  in  this  investigation  the 
values  of  recorded  voltage  rely  upon  the  centre-line 
velocity,  thus  'X'  refers  to  the  centre-line  velocity, 
and  *Y'  to  the  recorded  voltage  readings. 

(4.5)  and  (4.6)  are  the  two  main  equations  for  estab¬ 
lishing  the  law. that  fits  the  best  straight  line, 
whereas  (4.7)  gives  an  indictation  as  to  whether  or 
not  statistically  the  straight  line  obtained  is  a 
good  fit  to  the  experimental  data.  The  value  obtained 
from  this  equation  should  be  approximately:  0.98  -  1 
for  good  fit. 

Por  the  two  Boundary  Layer  probes  used,  the  straight 
line  relationships  obtained  are: 

Probe  1 

v  =  0.301263V  +  0.1392472  ‘  (4.8) 
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Probe  2. 


v  =  0.303719V  *  0.2643  (4.9) 

where  v  =  velocity 

and  V  =  recorded  voltage  reading. 

4.9  Probe  Traversing. 

Once  the  probe  had  been  calibrated  and  calibration 
data  obtained,  i.e.  equations  (4.8)  and  (4.9)  for 
probe  1  and  2  respectively,  it  is  now  possible  to 
obtain  velocity  profiles  from  the  surface  under 
investigation . 

The  method  employed  here  required  the  probe  to  be 
traversed  across  the  flow  from  the  surface  in  question 
and  voltage  data  obtained  at  pre-determined  positions. 
Initially  probe  1  was  located  in  the  centre  of  the 
channel  at  a  position  of  0  Ft.  along  the  test  surface, 
(total  distance  from  channel  entry  9ft.  6inches.) 
using  the  carriage  support  assembly. 

The  probe  was  now  positioned  as  close  to  the  surface 
as  possible  without  actually  touching  it,  this  being 
so  as  not  to  damage  the  delicate  probe  sensing  element. 
The  vernier  measuring  system  attached  to  the  probe 
traversing  carriage  assembly  was  now  set  to  zero  and 


this  taken  as  the  probes  datura  position,  the  initial 
setting  being  within  G.5mm  from  The  surface. 

By  turning  the  function  selector  on  the  5 5^10  unit  to 
'operate*  this  meant  that  the  probe  was  working  and 
experiencing  a  cooling  effect  due  to  the  moving  fluid. 

The  voltage  required  to  balance  the  bridge  circuit  was 
indicated  on  the  digital  voltmeter. 

The  data  logger  could  now  be  switched  on  and  sampled 
voltages  recorded  for  a  test  period  of  10  seconds, 
whereby  the  total  number  of  results  obtained  was  30. 

After  this  test,  the  function  switch  was  returned  to 
the  'stand  by*  position  and  the  probe  removed  from  the 
flowing  regime.  Inspection  of  the  probe  could  now  be 
carried  out  to  ensure  no  damage  had  cccured. 

This  procedure  wa3  repeated  for  probe  positions  of: 
0,1,2,3,4,5,6,7,8,9,10,15,20,25,30,35,40,45,50,75,100, 

150  and  200mm  above  the  surface  in  question.  The 
recorded  voltages  at  each  probe  position  are  as  indicated 
in  appendix  B. 

The  reference  (l.l),(2.1)  etc.  refers  to  probe  traverse 
1  using  boundary  layer  probe  1  and  probe  traverse  1 
using  boundary  layer  probe  2  respectively . 
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4.10  Probe  Contamination* 

Before  any  measurements  were  taken  using  the  boundary 
layer  probe  it  was  gently  brushed  using  a  dilute  solu¬ 
tion  of  acetic  acid,  this  removed  any  impurities  from 
its  surface  which  may  alter  the  heat  transfer  charac¬ 
teristics  of  the  sensing  element. 

Because  the  probe  operates  on  the  amount  of  heat  removed 
from  the  element,  it  is  extremely  important  that  the 
conductivity  from  the  element  is  not  impeded  in  any  way. 

Conductivity  impedance  will  result  in  irregularities 
taken  hence  making  them  invalid. 
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5.1  Introduction 


The  36th  European  Mechanics  Colloquium,  which  was  held 
at  Imperial  College  between  the  17th  and  19th  of  April 
1972,  was  the  first  such  meeting  to  be  devoted  entirely 
to  laser  Anemometry.  This  colloquium  was  the  first  of 
its  kind  whereby  researches  from  all  over  Europe  inter¬ 
ested  in  Physics,  Electronics  and  Fluid  Mechanics  with 
the  sole  purpose  of  trying  to  communicate  their  theor¬ 
etical  knowledge,  exchange  experience  in  the  use  of  and 
demonstrate  contributions  made  in  the  principles  and  use 
of  Laser  Anemometry  as  a  means  for  velocity  measurements 
in  fluid  flow. 

5.2  Basic  Principles. 

The  optical  heterodyne  detection  of  the  Doppler  shift 
is  basic  to  the  operation  of  the  velocity  meter.  The 
fundamental  concept  of  a  frequency  shift  in  radiation 
received  from  a  moving  body  by  a  stationary  detector 
has  been  understood,  and  used  in  communications  and 
astronomy,  for  example,  for  many  years. 

One  of  the  first  suggestions  that  mixing  of  two  signals 

of  different  frequencies  to  generate  a  difference  or 

cu) 

beat  frequency,  was  made  by  A. T. Forrester  in  1347  end 
1948.  This  he  published  in  the  Physics  Rev.,  Volume  73 
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If  we  consider  a  situation  whereby  two  incident  mono¬ 
chromatic  plane  light  beams  are  superimposed  on  the 
surface  of  a  detector,  i.e.  photocathode,  and  we  assume 
their  angular  frequencies  are  &)iand£Oa.  The  electric 
fields  will  be: 

Ea  =  E,  cos  cb.t  (5.1) 

Eg  =  Eacosui3t  (5.2) 

The  output  current  from  the  photomultiplier  tube,  I, 
is  proportional  to  the  square  of  the  electric  field 
incident  upon  it,  i.e. 

I  -  (2*  +  E*f  x  K  (5.3) 

where  K  *  constant. 

However,  the  photomultiplier  will  not  fellow  frequencie 
greater  than  several  hundred  megahertz,  therefore  terms 
in  the  expansion  of  (5*3)  involving  u). ,  tOa,  or  + 
will  give  rise  only  to  a  d.c  current  proportional  to 
the  time  average  of  those  terms.  If  however,  ( 0l>.  —  a)*) 
/27T  is  below  about  lOCHIHz  there  will  be  an  a.c  signal 
of  frequency  lA—cA*  . 

With  the  more  conventional  light  sources,  the  optical 
heterodyne  process  was  only  marginally  effective  due 


3 

to  the  relatively  large  band  widths  (10  Hz).  The  low 
intensity  per  unit  bandwidth  gave  extremely  poor  signal- 
to-noise  ratios,  and  therefore  ruled  out  detection  of 
Doppler  shifts  smaller  than  10*Hz. 

However,  the  major  difficulties  encountered  above  using 
monochromatic  light  were  later  overcome  with  the  advent 
and  development  of  the  laser. 

The  laser,  (light  amplification  by  stimulated  emmissicn 
of  radiation)  emits  a  continuous  coherent  wave  of  plane 
polarized  light,  this  due  to  very  high  power  levels 
being  contained  within  a  very  narrow  bandwidth  makes 
the  laser  ideal  for  use  as  the  source  of  light. 
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5«3  The  Doppler  Effect. 

The  Doppler  effect  is  named  after  Christian  Doppler  'who 
gave  an  explanation  in  1842  of  the  higher  frequency  of 
a  train  whistle  as  it  approached  an  observer,  and  the 
lower  frequency  as  the  train  went  away  from  the  observer. 

An  observer  on  a  particle  moving  away  from  a  fixed 
source  of  light  would  see  the  light  at  a  lower  frequency 
than  the  source  frequency,  and  to  a  stationary  observer, 
the  light  scattered  from  the  particle  would  also  appear 
to  have  a  different  frequency. 

The  geometry  of  the  situation-  is  shown  FIG.  5.1. 

V 


FIG.  5.1. 

Light  of  frequency  f^  is  incident  upon  a  moving  scatt¬ 
ering  center.  An  oscillating  field  will  be  induced  in 
the  scattering  center  by  an  electric  field  of  the 
incident  radiation. 
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The  frequency  of  this  oscillation  e^  is  given  by  the 
relationship: 

/  v®n 
f  =  fi  (  1 - TT) 

where  C  =  velocity  of  light. 

f  =  frequency  of  oscillation, 
f^  =  frequency  of  incident  light, 
e^  =  unit  vector  in  incident  direction. 

V  =  velocity  vector  of  the  scattering  center. 

Also  C  =Xf  (5.5) 

where  X  =  wavelength  of  light. 

To  any  observer  at  the  position  of  the  photomultiplier 
tube  this  frequency  of  scattered  radiation  is  shifted 
because  of  the  Doppler  Effect. 

The  frequency  fg  (frequency  of  scattered  light)  observed 
at  the  photomultiplier  tube  is 

f0  =  f  J  1  +  (5.6) 

Where  e^  and  e0  are  the  unit  vectors  in  the  incident 
and  scattering  direction  respectively. 
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By  eliminating  f  from  equations  ( 5 •  4 )  and  (5.6)  we  have: 


1  - 


1  + 


vli) 

■"S-' 

Ve 


Ve*3, 

1  +  IT) 


Ve 

s  vei 

U - C”  " 


eies 


V*«c  C*  therefore  the  last  term  in  the  bracket  can  be 
ignored. 

from  equation  (5«5) 

f  1 
V  9  X 


therefore : 


fi  +  tV‘(*s  -  *i> 


^.(ls  -  e,) 


(5.7) 


where: 

fs  -  =  fD  =  Doppler  frequency.  (5*8) 


therefore  from  equation  (5*7)  we  can  see  that  the  Doppler 
frequency  is  directly  proportional  to  the  particle 
velocity. 
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5.4  Early  Developments. 


One  of  the  first  uses  of  Laser  Anemometry  for  the 

measurement  of  fluid  velocity  flow  was  by  Y.Yeh  and 

(#) 

H. Z. Cummins  (1964),  who  measured  velocity  in  a  fully 
developed  laminar  pipe  flow.  Foreman  (1965-1966)  later 
refined  the  techniques  and  applied  then  to  the  measure¬ 
ment  in  gases. 

/ 1*7) 

Following  on  from  the  initial  work  done,  R.J. Goldstein 
and  D.K.Kreid  (1967)  used  the  method  to  measure  the  flow 
developed  in  rectangular  ducts  and  also  to  study  turb¬ 
ulence  intensity. 

One  of  the  most  commanding  attributes  of  Laser  Anenom- 
etry  is  that  unlike  the  hot-wire,  ho-c-fila  anemometers 
calibration  is  not  necessary,  also  equation  (5*7)  indi¬ 
cates  that  the  Doppler  shift  fs  -  £ ^  is  directly  pro¬ 
portional  to  the  component  of  velocity  in  the  direction 
eg  -  e^.  This  does  give  one  major  advantage,  and  that 
is  that  the  anemometer  can  be  used  to  study  and  take 
measurements  in  situations  of  high  turbulance  without 
the  complication  of  non  linear  characteristics. 

The  earlier  deve?.oped  systems  did  however,  have  major 
drawbacks,  and  these  were  mainly  associated  with  the 
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types  of  optical  system  used. 
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The  units,  as  used  by  Yeh,  Curtains,  Foremen,  George  etc, 
all  required  several  optical  components  such  as  mirrors 
and  beam-splitters  with  the  resulting  situation  that 
alignment  was  critical,  hence  necessitating  the  equip¬ 
ment  being  mounted  on  large  lathe  beds. 


Goldstein  and  Xreid  refined  the  original  systems  to  the 
one  as  shown  PIG.  5-2. 


NEUTRAL  DENSITY 
FILTER 


FIG.  5.2  GOLDSTEIN -KREID  SYSTEM. 

This  system,  however,  had  several  disadvantages  in  that 

alignment  was  critical,  also  a  distinct  encumbrance  was 

encountered  when  adjustments  were  required,  i.e.  traverse 

across  a  pipe  channel  etc. 
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Rudd,  however,  refined  the  system  even  further,  where¬ 
by  the  initial  laser  beam  passed  through  a  le.n3  system 
to  produce  a  wide  spread  of  coherent  light.  This  light 
was  then  masked  to  produce  two  single  beams. 

The  major  fault  with  this  system  was  that  due  to  the 
masking  much  of  the  laser  power  was  lost,  thus  necess¬ 
itating  lasers  of  very  high  intensity  to  be  used,  e.g. 
25-50  mV/. 

Traversing  of  the  units  was  again  made  difficult  due 
to  the  complexity  of  the  system. 


5.5  The  Optical  System. 


In  all  such  systems  previously  discussed,  the  critical 
factor  is  the  optical  unit.  Basically  there  are  three 
main  types  which  are  shown  in  PIG.  5»3,  5*4,  and  5*5. 


FIG.  5«3  shows  two  right  angle  prisms  and  a  reflecting 
mirror.  Partially  reflective  surfaces  are  made  so  that 
about  4-57^  of  the  incoming  light  is  transmitted  through 
and  the  rest,  about  95 is  reflected  and  brought  parallel 
to  the  former  by  a  reflecting  mirror. 

PIG.  5*4  and  5*5  shows  alternative  design  using  three 
right  angled  prisms,  and  one  right  angled  prism  with  a 
rhomboid  respectively. 
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PIG.  5*5  is  typical  of  the  system  used  at  The  National 
Engineering  Laboratory,  although  the  laser  beam  entered 
from  the  base  and  not  as  indicated  in  the  PIG. 

PIG.  5*3  is  typical  of  the  unit  to  be  used  by  the 
author  in  the  preceding  investigations. 

The  main  advantage  of  this  type  of  beam-splitter  arrange¬ 
ment,  is  that  it  can  be  housed  in  a  unit  which  permits 
the  change  in  direction  of  velocity  measurement  with¬ 
out  disturbing  the  optical  unit.  The  Lisa  55L01  type 
system  is  basically  that  of  PIG.  5*3 »  and  the  whole 
optical  unit  can  be  rotated  about  its  own  axis,  thus 
always  ensuring  that  the  optics  are  aligned. 


5.6  Modes  of  Operation 


In  order  to  utilize  the  measuring  principle  "based  on 
the  Doppler  effect,  there  are  three  main  modes  of  oper¬ 
ation,  these  are*. 

1)  Ref erence-Beam  Mode .  FIG.  5*6. 

2)  Differential  Doppler  Mode.  FIC-.  5.7. 

3)  Dual-Beam  Mode.  FIG.  5.8. 


5.6.1. 

In  the  reference-beam  mode,  FIG.  5.6,  the  laser  light 
is  split  up  into  two  beams  which  are  directed  at  the 
place  of  measurement  in  the  flow. 


FIG.  5.6 


The  place  of  beam  intersection  forms  the  measuring 
volume.  The  intensity  of  the  reference-beam  can  be 
reduced  by  introducing  a  neutral  density  filter  to 
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optimise  the  signal  picked  up  from  the  photomultiplier. 
The  photomultiplier  receives  two  signals,  one  with 
frequency  shift  and  one  without.  The  photomultiplier 
then  emits  an  alternating  current  with  a  beat  frequency, 
fD  equal  to  that  of  the  Doppler  frequency,  equation 
(5.7) 

5.6.2. 

The  Diff erential-Doopler  Mode  PIG-.  5«7  is  also  known 
as  the  ’fringe  mode'  because  the  beams  form  a  fringe 
system  in  the  volume  of  intersection. 


PIG.  5.7 


This  system  is  best  used  where  the  intensity  of  scatt¬ 
ered  light  is  low  and,  therefore  the  scattered  light 
can  be  picked  up  over  a  wide  angle  since  the  differential 
Doppler  frequency  is  independent  of  direction  and  detec¬ 
tion. 
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In  this  mode  the  photomultiplier  picks  up  scattered 
light  signals  from  the  same  direction  from  two  incident 
beams.  Each  scattered  beam  therefore  has  a  frequency 
shift  relative  to  the  incident  beam  that  it  originates 
from. 

The  differential-Doppler  mode  can  also  be  used  with 
backscattering,  however,  even  with  high  concentrations 
of  scattering  particles,  the  system  necessitates  the 
use  of  an  extremely  powerful  laser.  e.g.25m’'V. 

In  this  mode,  the  equipment  is  located  on  one  side  of 
the  flow  only,  hence  the  name  ' backscatter * .  One  useful 
method  whereby  the  'backscatt'eTing*  mode  can  be  used 
is  when  the  equipment  is  mounted  inside  the  hull  of  a 
ship,  and  measurements  are  being  taken  of  the  velocity 
profile  of  the  ship,  relative  to  the  water  it  passes 
through. 

5.6.3. 

The  Dual-Beam  mode,  PIG.  5*8,  is  where  the  laser  is 
situated  on  one  side  of  the  measuring  volume  and  the 
Optical  unit  and  photomultiplier  are  on  the  other. 


PLOW 


In  this  mode  the  aperture  disks  supplied  with  the  opti¬ 
cal  unit  are  used  by  attaching  the  disk  required  to  the 
front  of  the  beam-splitter.  The  twin  holes  allow  only 
two  scattered  beams  to  pass  through  the  lens  which  are 
combined  into  one  beam  to  be  picked  up  by  the  photo¬ 
multiplier. 

The  photomultiplier  receives  two  light  signals,  each  of 
which  are  scattered  in  a  different  direction  from  the 
same  incident  beam.  Both  signals,  however,  will  have 
frequencies  shifted  with  respect  to  the  frequency  of  the 
incident  beam.  Again  the  output  from  the  photomultiplier 
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will  be  of  a  signal  equal  to  the  beat  frequency. 

Due  to  the  relative  magnitude  of  the  laser  power,  the 
dual-beam  mode  is  not  to  be  used  in  these  investigations. 
The  velocity  measurements  are  to  be  confined  to  Refer¬ 
ence-beam  mode  and  dif ferential-Doppler  mode  with 
forward  scattering. 

The  systems  as  described  previously,  can  also  be  modified 
whereby  the  introduction  of  extra  beam-splitters  and 
photomultipliers  give  the  added  advantage  that  both  two 
and  three  dimensional  flow  measurements  can  be  invest¬ 
igated. 


The  laser  anemometer  measures  the  Doppler  frequency, 


fp  obtained  from  a  frequency  shift,  due  to  particles 
in  the  flow.  The  advantage  of  this  system  is  that 

particles  may  therefore  he  as  small  as  the  wavelength 

-9  tV 

of  the  laser  light,  i.e.  632.8  x  10  m  (0. 0000006328m) 

in  diameter. 

With  particles  as  small  in  diameter  as  that  shown,  it 
is  assumed  that  they  have  the  same  velocity  as  that  of 
the  flowing  medium,  even  in  regions  of  very  high  turb¬ 
ulence. 

?or  all  the  modes  of  operation  to  he  used  in  this 
investigation,  the  flow  velocity  can,  therefore,  he 
calculated  from  the  equation: 

V  =  (5.9) 

2  sin  | 

where  V  =  local  flow  velocity, 
fD  =  Doppler  frequency, 

A  =  wavelength  of  laser  light, 
and  ©  =  Beam  intersection  angle. 


5.8  light  Scattering  Particles. 


The  presence  of  particles  in  the  fluid  under  investiga¬ 
tion  is  essential  to  laser  anemometry.  V/a ter  however, 
normally  contains  an  ample  concentration  of  suitable 
particles,  whereas  atmospheric  air  and  gases  do  not. 
Seeding  must  therefore  be  introduced  of  suitable  size 
and  quantity. 

In  general  the  particles  must  be  small  enough  to  be 
able  to  follow  high  frequency  turbulence  situations 
with  accuracy  better  than  1$,  also  the  concentration 
of  such  particles  is  of  grave  importance. 

The  particles,  although  small  enough  to  follow  the 
turbulence  fluctuations,  must  be  large  enough  to  cause 
detectable  scattering.  The  correct  concentration  amount 
varies  depending  on  the  mode  of  operation. 

In  forward  fringe  mode,  the  scattering  particles  are 
lees  than  that  for  the  reference-beam  mode,  whilst 
backscattering  requires  very  much  higher  concentrations. 

Due  to  the  existing  system,  i.e.  flow  channel, and  the 
conditions  as  previously  mentioned  with  respect  to  rust 
particles  etc.,  transported  around  the  system,  seeding 
of  the  flow  is  not  required. 


5.9  Equipment  Supplied. 

The  Laser  anemometer  to  be  used  to  verify  results  obt¬ 
ained  from  the  hot  film  anemometer  was  purchased  from 
Lisa,  and  represents  Lisa  55L  Laser  Anemometer  Mark  1 
system.  This  basically  consists  of  six  separate  items: 

1)  10m’<V  maximum  intensity  Hughes  Laser  class  3B.  with 
exciter . 

2)  55L01  Optical  Unit. 

3)  55102  Blow  Direction  Adapter. 

4)  55L10  Photomultiplier. 

5)  55L15  High  Voltage  Supply. 

6)  55L20  Doppler  Signal  Processor. 

These  will  be  discussed  separately. 


5.9.1. 

Laser. 

The  Laser  recommended  to  fulfill  our  requirements  was 
a  lOnW  maximum  intensity  Hughes  Helium-Nean  type  class 
3B.  The  usual  recommended  laser  intensity  for  operating 
in  the  forward  scatter  modes  is  5mW,  therefore,  the 
above  should  give  good  performance  results. 

The  Laser,  before  being  used,  was  checked  fcr  intensity 
and  was  found  to  be  6mW,  therefore  this  unit  is  not 
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powerful  enough  to  be  used  in  the  ba ckscattering  node 
of  operation. 

5.9  (2&3) 

Optical  Unit  including  Flow  Direction  Adapter. 

The  main  function  of  the  optical  unit  is  to  produce  two 
coherent  beams  of  light  from  the  initial  single  source. 
This  type  of  optical  unit  is  exceptionally  versatile 
and  can  be  used  in  all  the  modes  of  operation  as  dis¬ 
cussed  in  section  5.6,  except  2-Dimensional  as  this 
requires  a  supplementary  unit. 

Three  focal  length  lenses  are  also  supplied  which  give 
working  ranges  of  13cm,  30cm  and  60cm  respectively, 
thus  allowing  the  intersection  angle  of  the  two  laser 
beams  to  be  continuously  adjusted  within  a  1:3  range, 
e.g.  between  8  and  24  degrees  beam  angle  for  a  working 
distance  of  13cm. 

In  the  investigations  of  tfie  velocity  profile,  the  lens 
used  has  a  focal  length  of  30cm.  This  was  necessary  be¬ 
cause  of  the  channel  width,  i.e.  30.5cm.  The  lens  used 
also  had  minimum  and  maximum  beam  angles  of  3«4  to  10.3 
degrees  respectively. 
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The  biprism  and  adjustable  surface  mirror  inside  the 
bousing,  constitutes  the  beam-splitter.  The  mirror  can 
be  moved  25*5ma  in  a  direction  55°  to  the  system  axis 
by  means  of  a  micrometer  screw.  The  position  of  the 
mirror  is  indicated  in  a  small  window,  and  by  this 
adjustment,  the  beam  separation  can  be  varied  between 
18  and  54mm  respectively. 

The  two  beam  adjusters  located  on  top  of  the  unit,  en¬ 
able  perfect  intersection  of  the  beams  at  the  point  of 
measurement.  These  beam  adjusters  are  2  glass  wedges 
with  wedge  angles  of  3°  .These  are  positioned  next  to 
each  other,  but  their  tapers  are  at  90°  to  each  other. 

Now,  by  adjusting  the  micrometer  screws,  thi3  permits 
rotation  of  the  wedges,  ensuring  a  high  resolution  of 
beam-adjustment  for  perfect  beam  intersection. 

The  unit  also  comes  complete  with  an  aperture  disk  for 
measurement  in  the  dual-beam  mode  on  forward  scattering, 
however,  it  was  the  authors  intention  not  to  use  this 
mode  of  operation. 

PIG.  5«9  indicates  a  schematic  diagram  of  the  optical 
and  flow  direction  adapter. 
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PIG.  5-9 

The  flow  direction  adapter  is  connected  to  the  optical 
unit  and  contains  two  optical  Bragg  cells  mounted  in 
parallel.  These  cells  are  in  a  fixed  position  but  their 
angular  position  may  be  adjusted  to  allow  the  laser 
beams  to  pass  through  undisturbed. 


The  Bragg  cell  consists  of  a  block  of  acousto-optic 
material.  Along  one  side  of  each  of  these  blocks  is 
bonded  a  series  of  piezoelectric  transducers.  Upon 
the  application  of  an  RP-vcltage  across  the  tranducers 
a  train  of  plane  acoustic  waves  is  created  in  the  mat¬ 
erial  moving  through  the  cell  with  the  velocity  of 
sound. The  moving  waves  are  absorbed  by  a  damping  material 
on  arrival  to  the  other  side  of  the  block.  The  acoustic 
wave  fronts  in  the  cell  constitute  a  thick  moving  opti¬ 
cal  grating  which  diffracts  the  light  passing  through 
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the  cell.  One  cell  shifts  one  input  beam  by  a  reference 
f2  plus  an  increment  (fs/2)  and  the  other  by  a  refer¬ 
ence  minus  an  increment  (fs/2),  therefore  the  net 
frequency  shift  is  fQ,  This  frequency  shift  may  be 
used  to  measure  highly  turbulent  velocities,  and  since, 
the  Doppler  frequency  fluctuation  caused  by  the  turbul¬ 
ent  motion  do  not  change,  therefore  the  relative 
frequency  range  is  reduced  as  the  zero  shift,  fQ  is  in¬ 
creased. 

The  effect  of  the  frequency  shift  is  thus  to  produce  a 
detector  current  of  frequency  fQ  when  measuring  a  par¬ 
ticle  at  rest  in  the  measuring  volume,  whereas  the 
frequency  is  increased  if  the  particle  moves  to  one  side 
and  decreased  when  it  moves  to  the  other. 

Ey  means  of  the  switch  on  the  front  of  the  flow  direc¬ 
tion  adapter  driver,  seven  frequency  values  fQ  may  be 
selected,  each  of  these  frequencies  corresponds  to  the 
midscale  value  on  the  55L35  frequency  tracker  unit. 

If  a  digital  voltmeter  is  connected  to  the  output  from 
the  meter  unit  without  frequency  shift  the  scale  read¬ 
ing  on  the  analog  meter  corresponds  to  0-1QV  dc  in  all 
frequency  ranges.  However  with  frequency  shift  then  5V 
corresponds  to  zero  velocity,  also  the  sign  on  the  unit, 
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♦or  -  indicates  flow  direction 


5«9«4. 

Photomultiplier. 

The  photomultiplier  detects  the  Doppler  shifted  light 
scattered  from  the  point  of  measurement.  PIS.  5»10 
shows  a  diagramatic  view  of  the  unit. 


SHUTTER  VIEW  FINDER 


The  photomultiplier  is  located  in  the  main  housing  and 
is  protected  from  light  by  a  built-in  shutter  controlled 
from  a  knob  at  the  side. 

A  pin-hole  aperture  in  front  of  the  shutter  prevents 
undesired  light  from  hitting  the  HI.  tube  during 
measurement  thereby  optimizing  the  signal  quality.  Two 
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sprung  adjusting  screws  enable  the  pin-hole  aperture 
to  he  accurately  centred  on  the  scattered  light,  thus 
allowing  adjustment  for  maximum  possible  quality  and 
minimum  possible  ambient  unwanted  light. 

The  unit  consists  mainly  of  a  tube  section  and  focusing 
optics,  and  three  close-up  lenses  are  also  supplied  with 
focal  lengths  of  200,  333  and  600mm  respectively. 

It  is  important  that  the  correct  close-up  lens  is 
attached  to  the  photomultiplier  to  match  that  as  inst¬ 
alled  in  the  beam-splitter.  Correct  lens  installation 
ensures  optimum  performance  from  the  photomultiplier 
and  allowing  use  of  the  same  pinhole  and  aperture  for 
all  ranges  of  focal  lengths. 

Situations  which  occur  whereby  the  photomultiplier  is 
placed  differently  to  that  of  the  laser,  will  result 
in  the  pinhole  requiring  to  be  changed. 


5.9.5. 

High  Voltage  Supply. 


The  high  voltage  supply  is  the  driving  unit  for  the 
photomultiplier  and  supplies  a  continuously  adjustable 
D.C.  voltage. 


located  on  the  front  panel  are  two  meters  which  indicate 
photomultiplier  tube  anode  current  and  high  voltage 
supplied  respectively.  The  supply  can  be  adjusted  from 
0  -  2000V  dc,  and  has  an  accuracy  range  of  25# »  whereas 
the  anode  current  meter  range  is  0  -  lQOyuA. 

With  this  unit  care  must  be  taken  that  the  correct  high 
voltage  cable  is  connected  from  the  sxipply  to  the  rear 
of  the  photomultiplier  tube.  This  cable  has  modified 
ENC  sockets,  and  is  denoted  with  a  red  covering  at  its 
ends.  This  modified  ENC  connector  is  to  prevent  inter¬ 
changing  of  the  high  voltage  and  signal  cable  connec¬ 
tions  thus  causing  grave  damage  to  the  photomultiplier 
tube. 


5.9.6. 

Signal  Processor. 

The  Doppler  Signal  Processor  consists  of  the  following 
main  circuits: 

a)  55130  Preamplifier, 
to)  55L35  Frequency  Tracker, 
c)  55140  Meter  Unit, 
a)  Preamplifier. 

The  preamplifier  input  socket  is  connected  to  the  out¬ 
put  of  the  high  voltage  supply  unit.  As  the  name  implies 
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•the  unit  amplifies  and  filters  the  signal  received  from 
the  photomultiplier  tube.  Upon  setting  of  the  frequency 
range  on  the  frequency  Tracker,  the  built-in  filters 
in  the  preamplifier  are  automatically  switched  for  top 
and  bottom  signal  cutting. 

b)  The  Frequency  Tracker  can  be  both  manually  and  auto¬ 
matically  controlled.  The  manual  control  is  utilized 
during  the  initial  setting  of  the  equipment  and  then 

by  switching  to  'auto'  the  unit  tracks  the  Doppler 
frequency  automatically. 

This  unit  contains  the  electronics  for  optimising  the 
signal  to  noise  and  ratio  and  also  has  the  frequency 
range  selector  which  gives  seven  frequency  ranges  2.25 
K  Hz  to  15  K  Hz. 

c)  The  last  piece  of  equipment  is  the  Keter  Unit.  This 
has  2  ENC  output  sockets  located  in  parallel  and  the 
Keter  indicates  the  Doppler  frequency.  The  meter  gives 
an  analog  reading  of  accuracy  1  f*  f.s.d.,  however,  by 
connecting  one  of  the  BNC  sockets  to  a  digital  volt¬ 
meter  then  0-10V  refers  to  full  scale  deflection  in  all 
frequency  ranges.  This  voltage  therefore  provides  an 
instantaneous  indication  of  the  Doppler  frequency  and 
hence  the  flow  velocity. 
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DESIGN  CP  IDA  SUPPORT. 


6.1  Design  Considerations  of  ISA  Support  System. 

Before  the  equipment  as  described  in  sections  5.9  could 
be  used,  it  was  necessary  to  design  a  support  system 
whereby  the  laser,  beam-splitter  and  photomultiplier 
could  be  mounted. 

With  regards  to  the  existing  channel  and  bearing  in  mind 
that  it  was  hoped  to  take  velocity  measurements  from  the 
replica  of  the  ships  hull  to  the  mainstream  of  the  flow, 
as  with  the  CTA  unit,  some  means  of  traversing  had  to  be 
allowed  for. 

Initially  there  were  two  possible  design  routes: 

1)  The  support  for  the  units  could  be  attached  to  the 
existing  channel  main-frame.  This  was,  however,  dec¬ 
ided  against  because  of  the  restriction  imposed  from 
the  ancillary  channel  leveling  equipment.  Supporting 
the  units  in  this  manner  would  have  restricted  the 
traversing  along  the  channel  to  approximately 

1  meter. 

2)  The  second  alternative  was  to  mount  the  main  support 
onto  one  of  the  existing  traversing  carriages. 


Plate  5  indicates  the  final  designed  support  system 
drawings  of  which  can  be  seen  in  Appendix  F. 


PLATS  5  SHOTTING  LASER  SUPPORT 
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6.2  Design  of  IDA  Support  System. 

It  had  been  initially  intended  to  mount  all  three  units 
separately,  whereby  the  brass  support  rods  would  pass 
through  the  main  and  sliding  supports  respectively. 
These  would  then  be  clamped  into  position  using  a  cap 
screw  at  the  side.  Upon  final  receipt  of  the  equipment 
this  is  how  the  LDA  was  initially  set  up. 

Prom  this  initially  set  up  the  design  progressed  where¬ 
by  the  brass  rods  were  threaded  and  six  brass  adjusting 
nuts  were  made.  This  enabled  traversing  in  the  vertical 
plane  to  be  controlled  and  made  with  accuracy. 

The  critical  alignment  between  the  laser  and  the  beam¬ 
splitter  was  not  fully  realised  until  the  equipment 
was  working,  whereby  it  was  decided  that  if  the  laser 
was  attached  to  the  beam-splitter,  then  the  critical 
alignment  between  the  two  would  no  longer  be  necessary. 
Also  by  having  the  laser  and  beam-splitter  attached 
resulted  in  the  system  being  much  more  stable,  and 
during  traversing  the  unit  moves  as  one. 

Drawings  of  the  designed  frame  including  optical  unit 
back  plate  etc.,  are  as  noted  in  appendix  P. 
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To  enable  measurements  to  be  taken  in  both  the  vertical 

and  horizontal  plane,  the  beam-splitter  must  be  capable 

0 

of  being  rotated  through  20  .  However,  for  vertical 
measurements  this  presented  no  problem,  but  in  the 
horizontal  plane  the  laser  is  not  in  line  with  the 
central  axis  of  the  optical  unit.  Therefore  the  laser 
support  bar  had  to  be  cut  and  an  adjusting  bar  made  to 
account  for  the  eccentricity  between  the  two  units. 

In  the  investigation  to  be  carried  out,  it  was  envisaged 
that  only  velocities  in  the  horizontal  plane  would  be 
obtained,  therefore,  once  the  laser  /  beam-splitter  was 
set  then  this  would  not  be  moved  other  than  for  vertical 
traversing. 
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REFRACTION . 


7.1  Laws  of  Refraction. 

The  flow  velocity  can  he  calculated  from  the  equation 
(5.9)  of  section  5.7. 


i.e.  V  =  fp  x  X 

2  sin  ^ 

where  V  =  local  flow  velocity, 
fp  =  Doppler  frequency. 

-8 

X  =  wavelength  of  He-Ne  Laser  =  632.8  x  10  m. 

0  =  3eaa  intersection  angle,  degrees. 

However,  before  any  measurements  are  taken,  it  is 
necessary  to  take  into  consideration  the  fact  that  the 
laws  of  refraction  must  be  strictly  adhered  to. 

There  are  two  basic  laws  of  refraction,  and  these  are: 

1)  The  incident  ray,  refracted  ray  and  normal  are  all 
'coplaner'  at  the  point  of  incidence. 

2)  Por  any  two  substances. 

=  constant  ,na  (7.1) 

where  n  =  refractive  index  of  medium  (2)  with  respect 
to  (1). 


This  is  more  commonly  known  a3  ‘Snells  Law'  . 


PIG.  7*1  indicates  the  equation  (7.1) 


PIG.  7.1 


If  we  consider  PIG.  7*2 


Now  toy  considering  the  laws  of  refraction  and  from 
equation  (7.1) 


If  we  consider  a  toeam  of  light  passing  from  air  to 
water,  then  the  wavelength  of  the  light  will  slow  down 
as  it  travels  through  the  various  mediums. 

X«  *  nT  X  (7.2) 

where  Xi  =  wavelength  of  laser  light  in  glass, 
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and  X  = 

632.8  x  10 

m. 

also, 

X*  = 

(7.3) 

therefore 

» 

Xj  - 

fcilJri  > 

Xz  = 

(7.4) 

now, 

sin  1'  = 

□  6 

7T  S1D  o 
n,  2 

(7.5) 

and, 

sin  = 

n7  8iD  f 

therefore 

9 

sin  I*  = 

fe}{3r)'ta  1 

therefore , 


Prom  equations  (7*4)  and  (7.6)  it  has  been  shown  that 
both  the  angle  of  intersection  and  wavelength  of  laser 
light  are  independant  of  the  glass  sides  of  the  channel. 


In  order  to  establish  the  actual  position  of  the  bean 
intersection  within  the  flow,  this  can  be  done  by  meas¬ 
uring  the  width  of  the  two  beams  at  the  initial  point 
of  entry  on  the  glass  side  of  the  channel. 


therefore , 


tan  2*  = 


B 

TTT 


?vhere  a  =  half  the  width  of  the  channel, 
also, 


tan  ** 


(A-B) 

Txo 


(7.7) 


(7.8) 

(7.9) 


when  b  =  thickness  of  glass  side, 
therefore  A  =  3  +  (A-B) 
where  A  equal  to  the  distance  apart  of  the  two  laser 
beams  at  entry  into  the  channel. 


Therefore  it  is  possible  to  establish  the  intersection 
point  within  the  flowing  regime  by  fairly  accurate 
measurement  of  the  beam  distance  ’A'. 


Therefore  from  the  forgoing  conclusions  and  using  the 
initial  equation  (5*9) 


V  .  x  x 

2s  in 

This  equation  can  he  modified  to  take  into  account  the 
laws  of  refraction  whereby  V  can  be  calculated. 

therefore, 

I-n  ^ 

V  =  — - £ 

23in  i  (7.10) 


In  atmospheric  air  the  refractive  index  'n'  at  76Cma 
Hg  at  20®  C  tw 

=  1.000277 


If  however  the  flow  medium  is  water  as  in  these  inves¬ 
tigations,  then, 

n»  =  1.334  at  0°C 

The  temperature  dependence  of  'o'  in  water  is  approx. 

/jq) 

An  =  0.0001/*C 

Therefore  if  we  assume  that  the  water  temp,  is  10°C, 

j 

then , 
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na  =  1.334  +  10(0.0001) 
na  =  1.335 

where  n  =  new  refractive  index  at  10*0. 

Therefore  from  equation  (7.4)  the  wavelength  of  laser 
light  in  the  flowing  medium  can  he  calculated: 

Xa  a  2^20277  x  632<8  x  1Q-9 
=  474.13879  x  10*n. 
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BOUNDARY  LAYER  THEORY. 


PIG.  8.1  illustrates  a  real  fluid  flowing  past  a  surface 
or  solid  boundary. 


PIG.  8.1 


At  tbe  boundary  there  is  a  very  thin  layer  of  fluid 
that  adheres  to  tbe  surface  and  therefore  has  a  zero 
velocity  with  respect  to  tbe  surface. 

Some  distance  away  in  the  mainstream,  however,  there 
is  a  region  normal  to  the  surface  in  vhe  'y'  direction 
where  the  fluid  has  a  constant  velocity  'U'. 

Near  to  the  wall  there  is  a  region  normal  to  it  in 
which  there  is  a  velocity  gradient  where  'u'  varies 
from  zero  to  'U' ,  this  region  in  which  there  is  a 


velocity  variation  is  known  as  the  'boundary  layer', 
and  is  usually  designated  £  (delta). 

The  concent  of  a  boundary  layer  was  first  introduced 
in  1904  by  the  German  engineer,  Ludwig  Prandtl,  (1875- 
1953).  He  suggested  that  the  flow  may  be  considered  in 
two  parts : 

1)  That  in  the  boundary  layer  where  the  shear  stresses 
are  of  prime  importance. 

2)  Seyond  the  boundary  layer  where  (in  general) 
velocity  gradients  are  small,  therefore  the  effect 
of  viscosity  is  negligible. 

If  we  consider  the  boundary  layer  in  laminar  flow,  the 
fluid  can  be  said  to  move  in  layers  with  one  layer  tend¬ 
ing  to  slide  with  respect  to  its  neighbour. 

The  shear  stress  is  tf/m1 

therefore  'X.  =  f*. 

where  is  the  velocity  gradient  and^is  the  absolute 

viscosity. 

Shear  in  laminar  flow  is  caused  by  internal  friction 

r 

due  to  the  molecular  interaction  between  adjacent  layers. 
There  is  however,  no  transfer  of  fluid  between  these 
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layers 


Linear  Homentum. 

The  tern  'linear  momentum*  of  a  body  is  defined  as  the 
product  of  the  mass  of  the  body  multiplied  by  its  linear 
velocity.  In  laminar  flow,  the  momentum  transfer  between 
the  adjacent  layers  is  molecular,  however,  if  the  flow 
is  'turbulent'  there  is  molecular  friction  plus  an  addi¬ 
tional  interaction  due  to  the  momentum  transfer  of  fluid 
masses  between  adjacent  layers. 

The  term  'turbulance'  implies  'irregular  fluctuations', 
not  a  regular  motion. 

The  region  close  to  the  wall  or  surface  may  be  classed 
as  a  'turbulent  boundary  layer'.  Close  to  the  wall, 
however,  there  is  a  region  in  which  the  velocity  app¬ 
roaches  zero  and  the  action  between  layers  is  purely 
molecular  rather  than  turbulent.  Therefore,  although 
the  whole  layer  may  be  classed  as  turbulent,  the  small 
laminar  region  close  to  the  wall  is  termed  the  'laminar 
sublayer’.  This  laminar  sublayer  may  have  an  appreciable 
thickness,  or  it  may  be  so  thin  that  the  molecular 
viscous  action  is  extremely  small  in  comparison  with 
the  turbulent  mixing  of  the  central  region. 


8.2  Thickness  of  a  Boundary  Layer  &  (delta). 

Because  the  velocity  within  the  boundary  layer  increases 
to  the  velocity  of  the  main-stream  asymptotically  some 
arbitary  convention  must  be  adopted  to  define  the  thick¬ 
ness  of  the  layer. 

One  such  possible  definition  of  the  thickness  is  that 
distance  from  the  surface,  whereby  the  velocity  reaches 
99/*  of  the  freestream  velocity  'U'.  In  these  investiga¬ 
tions,  to  enable  this  distance  to  be  identified  clearly, 
the  boundary  layer  has  been  taken  to  extend  to  98/S  of 
the  freestream  velocity  'U*. 

8.3  Displacement  Thickness. 

If  we  consider  flow  over  a  flat  plate,  the  fluid  will 
shear  against  the  plate  due  to  no-slip  conditions 
causing  a  frictional  drag.  The  velocity  distribution 
will  show  a  smooth  drop-off  to  zero  at  the  wall.  If  we 
consider  the  boundary  layer  thickness  as  that  defined 
above,  which  is  the  thickness  of  the  shear  layer,  then 
streamlines  outside  the  shear  layer  will  deflect  amount 
b  'the  displacement  thickness'. 

Consider  FIG.  8.2. 
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PIG.  8.2. 

If  we  consider  the  velocity  'u*  parallel  to  the  surface 
at  a  distance  • y •  from  it,  then  the  volume  flow  per  unit 
width  through  an  element  y  in  two-dimensional  flow  is: 

u6y  (8.1) 

If  however  there  had  "been  no  boundary  layer  the  flow 
would  have  been: 

U  Sy  (8.2) 

Therefore  the  total  reduction  of  flow  caused  by  the 
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boundary  layer  is: 

r°° 

J(U-u)  dy  (8.3) 

Now  to  reduce  the  total  flow  of  a  frictionless  fluid 
by  that  of  equation  (8.3)  then  this  may  be  equated 
to  uS*, 


Therefore , 

Ui*  = 

&  *  jy  J  (U-u)  dy 

♦ 

JOO 

(1-  n)  dy  (8.4) 

Equation  (8.4)  is  the  formal  definition  of  the  boundary 
layer  displacement  thickness  &  and  bolds  true  for  any 
incompressible  flow,  whether  laminar  or  turbulent, 
constant  or  variable  pressure,  constant  or  variable 
temperature . 


i 


(U-u)  dy 
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8.4  Momentum  Thickness  Q  (theta). 


If  we  again  consider  the  previous  fig.,  the  fluid  passing 
through  an  element  of  the  boundary  layer  carries  momentum 
at  a  rate 


(pu  <$y)  u  per  unit  width. 

whereas  in  frictionless  flow  the  same  mass  of  fluid  would 
have  momentum 

(^u  <Sy)  U 


Therefore  for  constant  density  the  total  reduction  in 
momentum 


u)  u  ay 


(8.5) 


o 


Equation  (8.5)  represents  the  momentum  flow  under  fric¬ 
tionless  conditions  through  a  thickness,  0 


(8.6) 


Equation  (8.6)  is  the  defining  relation  for  a  second 
parameter,  the  momentum  thickness  0  . 


The  definition  holds  true  for  any  arbitrary  incompress¬ 
ible  boundary  layer,  but  if  we  are  not  assuming  that  of 

x 

a  flat  plate, 0is  not  equal  to  the  drag  divided  by  pU. 
Therefore,  this  means  that  momentum  thickness,  although 


useful  in  certain  empirical  correlations  is  by  no  means 
as  fundamental  a  quantity  as  &  . 

The  equations (8 . 4)  and  (8.6)  are  illustrated  by  FIG.  8.3. 

I 

O 

FIG.  8.3. 


9% 

It  should  be  noted  that  6  must  always  be  less  than  6 
and  the  ratio  b  /0  is  called  the  'shape  factor'  and 
usually  is  denoted  by  H. 

The  value  of  H  is  always  greater  than  unity  and  in  turb¬ 
ulent  flow  varies  from  about  1.3  to  2.5.  ^ 
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8.5. 

The  Two-Dimensional  Turbulent  Boundary  Layer  Equation* 

There  are  three  basic  laws  of  conservation: 

1)  Conservation  of  Mass  (Continuity) 

2)  Conservation  of  Momentum. 

3)  Conservation  of  Thermal  Energy. 


If  by  definition  a  boundary  layer  is  one  in  which  there 
are  large  lateral  changes  and  slow  longitudinal  changes 
in  flow  properties, 


and  is  «  u. 


and 


then  the  equation  for  continuity  and  Momentum  are: 


continuity , 

du. 

Zx. 

+ 

II 

0 

(8.7) 

Momentum, 

u  ^ 
bx. 

+  = 

u  iy  4 

dx. 

J- 

p  aa 

(8.8) 

By  integration  of  equation  (8.8)  across  the  boundary 
layer  this  yields  the  ordinary  differential  equation: 


is  ♦  (2 

a 


.  S  dU  _ 

+  H,D£  ' 


Is. 

pu * 


£i 

2 


(8.9) 


The  equation  (8.9)  is  known  as  the  Karman  momentum 


98 


■if 


integral  relation. 


At  the  minimum,  two  additional  relations  are  required 

1  \ 

to  supplement  the  Kerman  integral  relation,  e.g.  ( 8 .  S }  < 

This  was  first  demonstrated  by  Ludwig  and  Tiilmann  in 
1945,  who  correlated  experimental  skin— friction  data 
from  a  wide  variety  of  flow  conditions.  It  was  found 
that  an  effective  algebraic  correlation  existed  betweer 


tne  three  v 


aria bles  ©  ,  H  and  Cj  and  derived  the  empir¬ 


ical  relationship: 


_  .6- 2i6  .O.  eTSK 

Of  »  0.246  He,  .  10 

-0.2d8  -O.feOSH. 

Cf  =  C.  246  (^  )  .  10 


(8.10) 

(8.11) 


8.6  Use  cf  Momentum  and  displacement  Thickness. 

Once  the  velocity  distribution  has  been  obtained,  then  0 
and  6  can  be  calculated.  From  the  velocity  profile  for 
probe  traverse  (2.1)  ana  (2.2)  the  velocity  at  each 
distance  ' y '  from  the  base  was  obtained  and  divided  by 
the  freestream  velocity  *U' .  By  assuming  the  velocity 
distribution  follows  that  suggested  by  Prandtl  whereby: 


tj  =  (y/<$ ) 


(8.12) 


various  values  of  ft  were  chosen  whereby  the  value  oc- 
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r  * 

(V  , 

tained  from  ( y/& )  gave  the  beat  agreement  to  u/u  at 
that  particular  value  of  'y'.  The  results  of  which  can 
be  seen  in  table  £t  and  £1  in  appendix  E. 

It  is  now  possible  to  obtain  .accurate  values  of  9  and  d> 
by  integrating  the  velocity  distribution  from  0  -»  & 
using  the  summation  of  the  various  power  law  profile. 
(Calculations  of  which  are  on  the  following  pages.) 

Now  by  substitution  of  the  value  ©  ana  H  (shape  factor) 

into  the  empirical  relationship  proposed  by  Ludwig  and 

(3) 

Tillmann  in  1949,  whereby  Cf  may  be  obtained  directly: 

-  0.-263  -0.678H. 

Cf  *  0.246  (^  )  .10  (8.11) 


L  t 


■J.j* 


8.7. 

Calcula tioti  of  C ^  Using  the  Xiudwlcg—Tilxmau  .nmolr real 

Equation . 


To  calculate  the  Momentum  thickness  8  for  probe  traverse 
(2.1)  it  can  be  seen  from  table  E»  in  appendix  E.  that 
the  velocity  profile  follows  the  l/5th  power  law  to 
0.0425m.  and  the  l/4th  power  law  to  the  edge  of  the 
boundary  layer.  Therefore,  9  can  be  calculated  by  the 
evaluation  of  two  integrals,  from  equation  (8.6). 


o 

.0.0*1 5  lA  I*. 

01  =  J  Cfc)  ( '■  W  )  h 

e»  = 


eT  = 


eT  = 


6|  +  6; 


fc/A  -O.OS1 


0T  =  0.006232  m. 


The  displacement  thickness  o  can  now  be  calculated 
using  equation  (8.4) 
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4* 

6* 

4? 

4*. 


f(i-  C/s)r')'h 

0>o**S 

J(l-  (^)  )<% 

o  oft  t/ 

J(i- 


4' 


0-0*25 

=  &;  +  4 

# 

2 

_  4/J  ■ 

0-0*25 

•  [D  *  1 V" 

].  +  fr 

=  0.008484 

+  0.0104 

=  0.008712 

m. 

S/4.  O- 

b*  *  V*L 


S'/*  ^  0-052 
O.0*2S 


-  0.010172 


Shape 

Factor  H 

u 

<v» 

CD 

_  0.008712  _  t  5  r\nnc 

-  V:m'o2Tl  -  -1*39795 

Using 

equation 

(8.11) 

can  be  calculated 

(8.11) 

°f  - 

0.246 

-0.268  -0-578  H 

(u6/v>)  10 

°f  = 

0.246 

,  .  -0.425-/. 3<J7*S 

/  0.565  *  0,006232  j  .-s 

'  0. 0000013  /  • 

0.00334 

Xo 


/?  «  tOOO  K^/fn3 


“to 

=  <3.00334 

x  4  «  p  «  0.565* 

To 

=  0.533 

fJ/m1 

Ho 

11 

^iSj 

11 

JH  =,  O^Ogj.j/s 
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M.>*\ 


To  calculate  the  momentum  thickness  6  for  probe  traverse 
(2«2)  using  the  results  in  table  £2  appendix  E 
From  these  results  the  velocity  profile  follows  that 
as  shown 


'S.fc 

J  from  y  =  0  to  1.5  tarn* 

from  y  =  1.5  mm  to  17.5  mm. 

[g-)  from  y  =  17.5  mm  to  42.5  mm. 

IjJ  from  y  *  42.5  mm  to  52.5  mm. 

n\'^ 

\fj )  from  y  =  52.5  mm  to  65  mm. 

r* 

Therefore,  to  completely  evaluate  8  and  o  we  require 
the  addition  of  five  integral  equations  respectively: 

!•«< 

•■Hlf  f 

%  -  flV* (•'<?#*)  <L 

O.onS  ^ 

64=  /  <?*>  O-W'jj*. 

a.s«tS  0 

65=  rvo-H'/u 

9'«5W  0 

©m  =  8|  +  8*  +  83  ■+  ©4  +  ©5 


f«l  _  li  1L  1 

©I  =  [iC  ^  j4 


©■  =  0.0003616  m 


e2  =  l&  V*  7 
©2  =  0.0034708  m 

.  Wst 


k/<  %  _  O-OI’i" 

f5  ^  .5^.1 

U  3*  7  ^  J— 5 


o  r «-  u  cf  ^  I 

03  =  Lfcfc  5"'*6  '  76  ~P"**  j  0.6,75 

©3  =  0.0334258  -  O.0l(3e)0(,8  -  0>otl74(>8  -*  O>ooso7a2. 


.  m 

s  *j7/5 

a.O#2? 

7  7** . 

4  o.  oazS 

9  a  =  0.0023424 

.  _  *'«• 


©4  =  0.0414206-  0.0S442  44  -  0. 0325314  ■+  0.02S6I2S 
04  =  0« 0005723  m 

®5  ’  L5  s".  c  £»• 

6s  =  0.052  -  0. 0433333  -  0.03^8163  o.O^MSSl 

Os  =  0.0003055  m 


©T  =  0.0075526  m 
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Similarly: 


s;  =  s:  *  s:  ♦  st  *  s'*  s: 

<?-  ro-^u. 

0-0)7*  O 

a*- 

o-atfS 

Pd-WH* 

0.o$2S  V 


&*  -  - 
4/  -  [-3- 
<5s*  =  [d  - 

if’  -  [*- 
tf  -  [*- 


a**** 

si  t-  7 

M  S“*‘  J. 


_ 

•t  Oifrtf 

G  J’,/# 

-J  O'CO'f 

_ 

ft  u 

-tO.**!*" 

ft 

J  0.0)7$ 

«  *$'*  7 

0<J3f 

6  6"i  J 

0.0*2 S 

o.o^s- 

tmfc 

0-\fi 

£ 

1 - 

a.oSZ5 

=  0,0104699  m 


(5.  000  8507  m 

a  005S7I1  m 
0,  00352  1  m 

0.0006108  rn 

0.  0003163  m 
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8.8  Distribution  of  Velocity  in  Turbulent  Flow, 

Prandtl  is  credited  with  introducing  (in  1925)  the  con¬ 
cept  of  a  'mixing  length'.  For  molecular  motion  there 
la  defined  a  term,  mean  free  path  of  the  molecule,  which 
means  the  average  distance  a  molecule  travels  between 
collision.  Mixing  length  is  analogous  to  'mean  free 
path' .  This  concept  of  a  mixing  length  states  that  in 
the  turbulent  exchange  the  fluid  masses  displaced  per¬ 
pendicular  to  the  mean  flow  carry  momentum  over  a  cer¬ 
tain  length  perpendicular  to  the  flow. 

The  mixing  length  is  defined  as  the  average  distance 
perpendicular  to  the  mean  flow  covered  by  the  mixing 
particle.  The  turbulent  shear  stress  X  can  be  taken 
as  a  function  of, 

f(  f,du/dy,  1)  (8.13) 

where  1  =  mixing  length. 

du/dy  =  velocity  gradient. 
p  =  density. 

The  ratio  ^/p  has  units  of  mVsec*  and  usually  the 
expression  U,  =  (8.14) 

is  used  where  U,0  has  units  m/s. 


Prandtl  related  the  local  r.m.s.  fluctuations  to  the 
local  velocity  gradients  through  a  turbulent  length  scale. 


(8.15) 

(8.16) 

Where  i,  and  aTe  called  the  mixing  lengths  which  are 
analogous  to  the  molecular  mean  free  path  but  much  larger. 

In  practice  the  particle  does  not  move  a  distance  £  and 
then  suddenly  changes  velocity,  but  undergoes  a  gradual 
change.  However,  on  the  assumption  that  the  change  in 
velocity  Au  experienced  by  the  particle  in  moving  a 
distance  1  in  the  y  direction  is  1  du/ay  and  propor¬ 
tional  to  both  u1  and  v1  as  described  previously  in 
equations  (8.15)  and  (8.16). 

Now  %  (turbulent  shear)  =  ^-r"C—  =  -p u*  v'  (8.17) 

Area 

JPrandtl  then  deduced  that  u1  v1  «*  1*  (*\a/ay)X 

Now  by  absorbing  the  coefficient  of  proportionality  into 
the  l2  and  accounting  for  the  sign  of  with  ^u/ay 
using  equation  (8.17) 


and , 


*  ti  ay 
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Equation  (8.18)  is  known  as  the  Prandtl  mixing  length 
theory.  It  is  however,  dimensionally  correct  and  real¬ 
istic  physically.  The  mixing  length  1  is  strongly 
dependant  upon  the  distance  from  the  wall  y  and  the 
following  empirical  approximations  have  been  suggested. 

Very  near  the  wall.  ,  Suggested  by 

1  ~  y 

(inner  layer)  (Reichardt  1951) 

Fairly  near  the  wall. 

1  =  Ky  (Karman  1930) 

(overlap  layer) 

where  K  =0.41 

Far  from  the  wall. 

Inconstant.  (Prandtl  1942) 

(outer  layer) 

Karman  in  1930  introduced  a  further  concept  by  means 
of  similarity,  he  suggested  that  the  mechanism  of 
turbulent  friction  can  be  fully  characterized  by  a 
single  length  1  only  for  the  case  in  which  the  flow 
pattern  in  the  neighbourhood  of  every  point  in  similar. 

Karman  presented  the  following  expression  in  terms  of 


the  mean  flow  where: 


(8.19) 


where  K  =  The  internationally  accepted  constant  0.4-1. 
This  analogy  however,  does  have  loop  holes,  one  of 
which  is  with  regards  to  pipe  flow. 


The  velocity  profile  across  a  turbulent  boundary  layer 
can  be  split  into  three  regions;  the  'inner',  or  wall 
layer  where  viscous  shear  dominates,  the  'outer',  or 
wake  layer  where  turbulent  shear  is  dominant  and  an 
'overlap'  layer  where  both  types  of  shear  are  signifi¬ 
cant. 


Prandtl  deduced  that  the  mean  velocity  must  depend 
upon  the  wall  shear  stress,  the  fluid  properties  and 
the  distance  y  from  the  wall,  therfore  the  inner  law 
can  be  expressed  as: 


u  =  f(Te,  p,  %  )  (8.20) 

and  therefore, 

u/Uo  =  fCiJ’Vv  )  (8.21) 

where  m/a 
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For  the  outer  layer  Karaan  deduced  that  the  wall  tends 
to  act  oerely  as  a  source  of  retardation  reducing  the 
local  velocity  u  below  the  freeotream  velocity  U  inde¬ 
pendently  of  viscosity  yii  but  dependant  upon  the  wall 
shear  stress  and  the  distance  y  over  which  its  effect 
has  diffused.  Therefore  the  outer,  or  velocity  defect 
law, 


U-U  =  ) 

(8.22) 

U-u/U*  =  g(  y/J  ) 

(8.23) 

where  f  and  g  refer  to  function  off  respectively. 

For  the  overlap  region  it  is  possible  to  merge  the  two 
functions  together  because  these  two  laws  are  valid. 

This  is  the  overlap  region! 

uAo  =  A))  =  u/u®  -  g(yA)  (8.24) 

For  the  equation  above  to  be  true  it  can  only  be  so  if 
both  f  and  g  are  logarithmic  functions, 

where,  u/o.  Ue(^)  +Bjinner  variables  (8.25) 

and  U-u/u«  =  ”<  I<ge(^)  +A|Outer  variables  (8.26) 
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In  the  equations  (8.25)  and  (8.26)  K,A  and  B  are  purely 
dimensionless  constants. 

where  K  =  0.41  Von  Karma' ns  constant. 

B  =  5*0  Coles  constant  (1955) 

and,  A  varies  depending  on  whether  or  not  the  flow 
is  boundary  layer  flat  plate  or  Fipe  or  channel  flow, 

therefore  for  a  smooth  flat  plate, 

>  *  5-0 

u/u.  •  )  *  5.0  (8.27) 

The  constant  A  varies  from  2.35  for  boundary  layer  flow 
(flat  plate)  to  0.65  for  Fipe  or  channel  flow.  If  in 
the  case  A  =  0  this  would  mean  that  the  outer  layer  has 
completely  vanished  so  the  lo^rithmic  overlap  extends 
all  the  way  to  y  -S  .  therefore  for  A  =  2.35  means 
that  there  is  a  small  overlap  layer  and  for  A  =  0.65 
means  that  a  negligibly  small  overlap  exists. 

The  validity  of  the  inner  law  can  be  seen  by  plotting 
u/u0  against  loge(  ^57)  or  (^7)  for  probe  traverse  (2.1) 
and  (2.2)  graphs  in  appendix  £. 

From  these  graphs  the  inner  region  extends  between 


36  <  (  )<  90  and  the  logarithnic  overlap  would 

aeem  to  extend  from  90 <  (  )  <  270. 

The  outer  layer  turns  upwards  and  extends  from  (y^r*  )  > 

270,  whereby  for  traverse  (2.2)  a  strong  adverse  press¬ 
ure  gradient  exists. 
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8.9  Law  of  the  Wake. 

With  reference  to  the  graphs  Es  and  Efe  in  appendix  E  the 
profiles  clearly  differ  in  their  outer  behaviour  for 
probe  traverse  (2.1)  the  profile  is  fairly  linear, 
whereas  (2.2)  indicates  strong  adverse  gradients. 


Clauser  noticed  that  the  outer  law  should  be  modified 
to  include  pressure  gradient,  dP/dx,  and  he  also  re¬ 
placed  &  ,  by  8  ,  arriving  at  the  parameter, 

I3"  4*  £  (8.28) 

After  considerable  experimental  procedure,  Clauser 
showed  that  a  boundary  layer  with  variable  pressure 
gradient  but  with  6  constant  is  in  turbulent  equilibrium 
in  the  sense  that  all  the  gross  properties  of  that 
boundary  layer  can  be  scaled  with  a  single  parameter. 


Clauser  determined  the  most  relevent  thickness  parameter 
for  equilibrium  to  be  the  defect  thickness,  A 

=  A  (8.29) 


where,  A  =  /  Tu 


where , 


is  a  measure  of  the 
local  skin  friction. 
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Since  the  skin  friction  varies  with  x  it  follows  that 
H  is  not  constant  in  an  equilibrium  boundary  layer. 

This  can  be  seen  from  the  previous  calculations, 

probe  traverse  (2.1)  H  =  1.39795 

probe  traverse  (2.2)  K  =  1.386264 

Prom  further  work  by  Alber(l968)  G  was  adopted  as  a 
form  factor  rather  than  H  and  a  curve-fit  relation 
proposed  by  Nash  (1965)  related  G  to  the  Clauser  para¬ 
meter  (3  ,  thus  the  shape  factor  was  related  approximat¬ 
ely  to  the  pressure  gradient. 

One  of  the  major  weaknesses  to  Clauser  approach  was 
that  curves  produced  by  him  were  not  easily  readable, 
however,  Coles  (1956)  noted  that  deviations  of  the 
velocity  above  the  overlap  layer  when  normalised  by 
the  taaximum  velocity  deviation  (at  y  =£  )  would  be  a 
single  wake-like  function  of  y/cT  only. 

Coles  therfore  proposed  the  following  deviation, 


IL*  -  2.5  In  y*  -  5.5 

as  t  W  (^/<0 

(8.30) 

U*  -  2.5  l«  S*  -  5.S 

where  u4  =  u/u. 

(8.31) 

(8.32) 

and  U*  =  U/tu 
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(8.33) 

W  is  known  as  the  wake  function  and  is  normalised  to 
he  zero  at  the  wall  and  have  a  value  of  2  at  y  =  $  . 
Thus  Coles  proposed  the  following  curve-fit  approxima¬ 
tion  for  the  wake  function: 

W(  )  *  2sin*  (n/2.  (8.34) 

Therefore  ths  equation  (8.25)  was  modified  to  include 
this  function  resulting  in 

i-  (8-35) 

When  II  is  called  the  wake  parameter  and  is  related  to 
the  outer -variable  overlap  constant  of  (8.26)  by  H 
a  KA/2  where  for  a  flat  plate  II  ^0,5. 

Coles  wake  function  V /('V/£  )  has  since  been  modified  to 
the  empirical  modified  version. 

w(^)  -  f «%)’-  (8.36) 

However,  as  previously  stated  at  y  =  &  this  function 
has  the  value  2,  and  at  y»  O  ,  equals  0.  Therefore  near 
to  the  wall,  e.g.  y  =  0.004m.  the  wake  function  has 
the  value  0.045  for  probe  traverse  (2.1)  and  0.029 
for  probe  traverse  (2.2).  Therefore,  due  to  its  rela- 


116 


sidering  the  equation  (8.35). 


The  preceding  computations  are  valid  only  for  smooth- 
walled  channels.  Kougnness  of  the  surface  has  very 
little  effect  in  laminar  flow  but  influences  turbulent 
flow  due  to  the  laminar  region  being  fairly  thin. 

In  laminar  flow  the  viscous  region  extends  all  the  way 
to  the  edge  of  the  boundary  layer. 


The  average  roughness  height  k  is  such  that  k  = 
and  will  effect  the  inner  law  but  not  the  outer  law. 
Therefore,  from  the  graphs  Es  and  £4  appendix  E  for  slopes 
plotted  with  different  values  of  u0  then  the  graph 
which  compares  to  the  smooth  wall  profile  is  displaced 
an  amount  therefore  equation  (8.35)  can  be  rewritten 
to  take  into  account  this  shift, 

a*'  t  loge(^)  +  S.0  -  AB+  £  W(%)  (8.37) 

For  probe  traverse  (2.1)  the  graph  which  lies  parallel 
to  that  for  a  smooth  plate  is  plotted  using  u„  =  0.0241 
the  displacement  &8  =  1.658. 

Ua  *  K  lo^gC^)  +5.0  -1.658  ■+  ^  VlC^/s) 


(8.38) 


k. 


For  probe  traverse  (2.2)  the  graph  plotted  using 

Uo  a  0.024  gives  best  agreement  with  that  of  the  smooth 

plate  and  the  displacement  AB  =  2.118. 

1  =  i  (?£)  *  5.0  -  2.  us  4  £  (8.39) 

Therefore  the  displacement  AS  would  have  been  zero  had 
the  plate  been  smooth,  thus  this  value  can  be  related 
to  the  roughness  of  the  surface. 

(3} 

Prandtl  and  Schlichting  analysed  various  compiled 
roughness  configurations  obtained  by  Clauser  and  evolved 
the  approximate  roughness  representation, 

A6  *  k  \o$e  -  3.0  (8.40) 

A 6  *  t  [o&  -?.o  (8.41) 

AB  £  k  0. 3  ( ^ )  (8.42) 
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8.11  Calculation  of  usin;:  the  Inner  Law. 

Prom  toe  previous  section  using  equation  (8.38)  for 
probe  traverse  (2.1), 

E.  =  i  *S.c  -  '■  “8  neglecting  the 

wake  function. 

the  value  of  u0  that  satisfies  this  equation  is  0.0241, 
therefore,  from  equation  (8.14). 


to  =  (0.0241 f  x  P 

°t  =  iTU1 

=  (-glfT24.1..?*  =  0.00364 

r  £  -  u* 

using  equation  (8.42)  and  relating  £>8  to  the  roughness 
height 

=  K  °*M  V  ) 

1.658  =  2.439  logfc0.3(  ^  ) 

therefore  the  value  of  k  which  satisfies  the  above 
equation 


0.000355m 


For  probe  traverse  (2.2) 

St  =  *  log.(  ^  )  +  5.0  -  2.118 
Uo  c 

and  uo  =  0.024  m/s 

C*.  =  t°i °g£L  =  0.00361 

f  r  - 

as  previous, 

A6  =  2.439  log*  G . 3 (  ^  ) 

2.118  =  2.439  log t  0.3(  ^  ) 
therefore  the  value  of  k  which  satisfies  this  value  is 


0.00043a 


OPERATION  0?  LEA  SYSTEM. 


p 


9»1  Connecting  Up  Of  Laser  Electrical  Equipment. 

Cnee  the  three  units , (laser ,  beam-splitter  and  photo¬ 
multiplier)  are  attached  to  the  support  system  they  can 
be  connected  up  electrically  by  following  the  procedure 
as  laid  down  below. 

1)  Connect  the  High  Voltage  cable  supplied  to  the 
rear  of  the  High  Voltage  supply  Unit  (55L15)  and 
Photomultiplier  (55L10)  respectively. 

2)  Connect  the  *  signal  out'  on  the  photomultiplier 
to  the  'in*  socket  on  the  high  voltage  supply 
unit. 

3)  Connect  the  'out'  socket  on  the  high  voltage 
supply  unit  to  the  'in'  socket  on  the  front  of 
the  55L30  Preamplifier. 

4)  Connect  one  of  the  output  sockets  on  the  55L40 
Meter  unit  to  a  dc  digital  voltmeter.  (If  possible 
the  digital  voltmeter  sh>uld  have  the  facility 

of  an  exponential  time  constant  variation,  this 
will  allow  the  Doppler  frequency  to  be  averaged 
over  a  particular  time  interval). 

5)  The  equipment  can  now  be  switched  on  as  it  is 
ready  for  use. 
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6)  Note:  the  laser  should  not  be  switched  on  until 
measurements  are  to  be  taken. 

7)  The  photomultiplier  shutter  must  be  kept  closed 
so  as  not  to  damage  the  photo-tube. 


9.2. 

Adjusting  IDA  Equipment  t,o  Obtain  Doppler  Frequency. 

As  soon  as  the  position  of  the  laser/beam-splitter , 
relative  to  the  side  of  the  channel,  has  been  establi¬ 
shed  by  consideration  of  section  7  (laws  of  refraction) 
the  photomultiplier  can  be  adjusted  for  maximum  scatt¬ 
ered  light  entering  the  lens.  This  is  done  by  looking 
down  the  view  finder  and  adjusting  the  two  screws  on 
the  Pin-hole  aperture  for  maximum  effectiveness.  The 
shutter  on  the  photomultiplier  can  now  be  opened  using 
the  knob  placed  on  the  end  of  the  housing,  whereby  this 
allows  the  scattered  light  to  hit  the  R1  tube. 

The  high  voltage  control  knob  should  be  increased  until 
the  Anode  Current  Meter  reads  approximately  50  ^A, 
this  should  be  at  about  400  Volts,  depending  on  the 
amount  of  scattering. 

Now  the  SET  LEVEL  knob  on  the  Preamplifier  should  be 
adjusted  so  the  meter  reads  approx,  mid-scale,  which 
can  be  readjusted  when  the  correct  Doppler  frequency 
range  has  been  found. 

The  tuning  selector  switch  on  the  Frequency  tracker 
should  be  set  to  MAN  and  the  tuning  knob  turned  fully 
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clockwise.  Adjust  the  SET  LEVEL  and  ADJUST  THRESHOLD 
to  the  middle  positions  and  also  set  the  IE  bandwidth 
at  2°/c  of  frequency  range,  the  drop  out  detector  can 
now  be  adjusted  and  frequency  range  set  at  15.  MHz.  The 
Doppler  frequency  meter  should  show  full  scale  deflec¬ 
tion. 

Now  turn  the  tuning  knob  slowly  anti-clockwise  until 
the  IP  level  meter  indicates  a  peak,  if  no  peak  is 
obtained,  then  adjust  the  Frequency  range  to  the  next 
lower  value.  The  tuning  procedure  should  now  be  repea¬ 
ted  . 

The  frequency  selector  should  be  decreased  until  the 
IP  level  meter  peaks  and  the  Doppler  frequency  meter 
indicates  the  Doppler  frequency  of  the  signal.  This 
will  be  indicated  on  the  dc  voltmeter  whereby  0-10V 
indicates  meter  full  scale  deflection  in  all  frequencies. 

If  however,  the  IP  level  meter  does  not  peak,  this  may 
mean  that: 

a)  the  flow  velocity  is  zero. 

b)  there  is  poor  beam  intersection. 

c)  scattering  of  the  particles  is  too  low. 

d)  Photomultiplier  has  been  incorrectly 
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focused  on  the  measuring  volume, 
e)  the  flow  may  be. highly  turbulent. 

After  manual  setting,  set  the  tuning  selector  to  AUi'G 
and  the  Doppler  frequency  is  automatically  tracked. 

The  SET  LEVEL  knob  can  now  be  readjusted  so  the  I? 
level  meter  reads  mid-scale. 

If  the  velocity  range  varies  quite  considerably  over 
the  measuring  volume,  then  the  Doppler  frequency  may 
go  out  of  frequency  range,  the  tracker  will  drop  out 
of  lock  and  resetting  is  required. 

If  the  velocity  of  the  fluid  is  known  approximately, 
then  the  correct  frequency  range  can  be  set  initially. 
By  traversing  across  the  flowing  regime  the  Doppler 
frequencies  can  be  obtained  at  various  positions  from 
the  surface  and  hence,  using  equation  (7-10),  the 
velocity  can  be  calculated. 

However,  due  to  the  late  arrival  of  the  equipment,  no 
time  was  available  to  complete  this  work. 


Plates  7  and  8  indicate  the  intersection  of  the  two 
laser  beans  as  they  pass  through  the  channel  for  neas 
ureiaents  in  the  vertical  and  horizontal  plane  respec¬ 
tively. 
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9«3  Safety  Precautions 


Because  of  the  light  source  being  what  it  is,  that  of 
a. laser*  it  is  extremely  important  that  certain  safety 
precautions  are  taker . 

firstly,  the  laser  beam  oust  never  be  directly  looked 
into  because  the  harmful  radiation  will  damage  the 
delicate  cells  in  the  eye* 

Secondly,  if  the  beam  hits  the  eye  directly,  then  due 
to  the  eyes  focusing  ability  this  magnifies  the  power 
of  the  laser  approximately  1QQ000  times  over  a  minute' 
area  killing  the  sensitive  cells  that  the  focused  beam 
hits*  This  magnified  power  is  of  the  order  of  6C0  Watts. 

Therefore,  because  the  working  region  is  at  approxim¬ 
ately  eye  level,  it  would  be  beneficial  to  provide  some 
form  of  'platform'  along  the  working  region  of  the 
channel  ensuring  that  any  future  operators  are  at  least 
working  above  the  line  of  action  of  the  leaser* 

Due  to  the  two  laser  beams  passing  in  and  out  of  the 
channel,  this  means  that  there  are  four  possible 
reflective  surfaces  to  which  the  beams  pass  through. 

When  operating  the  unit  it  was  noticed  that  due  to 


reflection,  there  was  a  region  approx.  12inches  either 
side  of  the  initial  entry  point  of  the  two  beams, 
whereby  reflected  laser  light  emitted, 

Therefore,  to  prevent  the  author  being  hit  in  the  eye, 
a  timber  screen  was  attached  to  the  beam-splitter  which 
extended  this  spread  and  shielded  him  from  the  stray 
beams. 
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9«4  Future  Developments. 

It  is  envisaged  that  this  project,  using  the  LBA  will 
be  continued.  However,  the  support  system  designed  was 
only  that  of  a  temporary  set  up  in  order  to  verify  tut 
-CTA  results,  and  it  is  hoped  that  a  future  sophisticated 
support  system,  incorporating  automatic  -traversing  etc. 
will  be  built.  This  future  support  will  allow  the 
laser/beam-splitter  and  photomultiplier  to  be  positioned 
at  various  angles  whereby  the  two  laser  beams  can  be 
focused  on  the  surface,  hence  accurate  velocities  at 
the  surface  can  be  obtained. 

With  the  present  system,  the  closest  one  can  measure 
from  the  surface  is  approx.  0.5om.  However,  there  are 
certain  set-backs  using  the  existing  channel,  and  these 
are : 

1)  Because  the  existing  glass  sides  are  plate  glass 
then  there  is  quite  a  large  amount  of  refraction 
of  the  two  laser  beams  at  inlet  and  entry. 

2)  If  the  beam-splitter  is  to  be  set  at  various 
angles,  then  alignment  of  the  two  glass  sides 
will  be  critical. 

3)  Because  of  the  beam-splitter  set  at  an  angle, 

then  this  may  mean  that  one  beam  enters  the  channel 
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before  the  other,  this  v/ill  result  in  the  two 
beams  being  out  of  alignment  which  may  be  out  of 
the  range  of  the  beam  adjusters. 

4)  If  the  units  are  to  be  used  other  thao  in  the 
horizontal  .plane,  then  a  working  region  in  the 
channel  will  have  to  be  decided  upon,  and  the 
existing  glass  sides  should  be  replaced  with 
Schlieren  glass  to  obtain  best  possible  results. 
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DISCUSSION 


10.1  Equipment  Used. 

From  the  previous  sections,  3.1.1  to  3.1.4,  it  is  quite 
evident  that  there  are  only  two  types  of  equipment  suit¬ 
able  for  close  boundary  layer  work.  These  are  the  Con¬ 
stant  Temperature  Anemometer  and  the  laser  Doppler 
Anemometer. 

All  other  types  of  equipment  cause  large  scale  disturb¬ 
ance  to  the  flow  invalidating  any  results. 

The  Propeller  type  current  meter,  however,  is  a  very 
useful  piece  of  equipment  because  it  does  give  a  quick 
indication  of  the  velocity  in.  a  freestream.  This  was 
used  to  establish  that  there  is  a  central  core  within 
the  channel  where  the  flow  is  uniform  potential.  The 
proof  of  this  can  be  seen  from  appendix  D  FIG.  Dl. 
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10.2  Difficulties  Observed  Using  the  CTA. 


This  equipment,  because  of  its  relative  size,  makes  it 
ideal  for  close  boundary  layer  work,  there  is  a  minimum 
flow  disturbance,  both  between  the  probe  and  the  sur¬ 
face  and  around  the  probe.  However,  there  is  some  dis¬ 
turbance  no  matter  how  slight,  which  therefore  causes 
the  validity  of  results  obtained  very  close  to  the 
surface,  to  be  treat  with  caution. 

The  CTA  relies  upon  the  amount  of  heat  conducted  away 
from  the  sensing  element,  therefore  when  working  very 
close  to  steel  surfaces  there  will  be  a  higher  conduc¬ 
tion  rate  than  when  working  close  to  plastic  surfaces, 
thus  again  results  obtained  depending  on  the  environ¬ 
ment  must  be  treat  with  suspicion . 

It  should  be  remembered  that  the  concept  of  the  bound¬ 
ary  layer  probe  being  so  minute",  is  that  velocities 
can  be  obtained  as  near  to  the  surface  as  possible. 

During  the  velocity  traverse  using  probes  1  and  2  the 
linearizer  was  constantly  connected  to  an  oscilloscope 
so  that  the  turbulence  fluctuations  could  be  observed. 
Prom  the  results  for  probe  1  and  2  calibration,  the 
voltage  results,  although  taken  through  a  varying 
range  of  flow  rates  with  Re  varying  between  10615  to 


2840  for  probe  1  and  13530  to  3438  for  probe  2  are  fairly 
steady  with  no  large  scale  fluctuations.  However,  upon 

observation  of  results  obtained  for  the  various  probe 

£ 

traverses  where  Re  is  in  the  region  1.79  x  10  to  2.053 
x  10  there  are  large  voltage  fluctuations.  The  Re  trans- 
ition  range  for  a  flat  plate  is  approx.  2  x  1G*  depending 
upon  roughness,  therefore  it  has  been  assumed  that  the 
values  above  are  within  the  Re  transition  range.  This 
turbulence  intensity  was  observed  on  the  oscilloscope  and 
it  was  also  noticed  that  the  voltage  recording  equipment 
did  not  seem  to  keep  pace  with  the  response  from  the 
probe.  Therefore  it  was  decided  to  build  a  small  RC  inte¬ 
grator  with  a  time  constant  of  approx.  10  seconds  using 
solid  state  resistors  and  capacitors.  V/hen  connected  to  a 
standard  type  ac/dc  digital  voltmeter  this  then  integrated 
the  voltage  signals  from  the  linearizer  giving  a  single 
value  of  mean  integrated  voltage  and  not  a  fluctuating  one. 
After  comparing  the  results  from  this  unit  with  that  ob¬ 
tained  from  the  data  logging  system  the  two  sets  of  voltage 
did  compare  quite  favourably.  Also  a  Solartron  digital 
voltmeter  was  used  with  the  facility  for  varying  the  expo¬ 
nential  time  constant  and  the  results  obtained  from  that 
compared  with  the  results  from  the  data  logger.  One  other 
possibility  however  exists,  and  that  is,  the  probes  record¬ 
ing  equipment  cannot  keep  pace  with  the  frequency  response 
from  the  probe  due  to  the  turbulent  intensity.  If  this  is 
so  then  improvements  in  the  response  of  the  unit  will  be 
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required.  It  should  be  noted  that  all  voltage  readings 
recorded  on  the  data  logger  during  a  test  were  accepted 
as  being  required  by  the  bridge  for  balance  even  though 
some  readings  seemed  either  too  high  or  too  low  compared 
with  the  rest.  As  mentioned  in  section  4.10,  Probe  Con¬ 
tamination,  before  being  used  the  probe  was  gently  bru¬ 
shed  using  a  dilute  solution  of  acetic  acid.  This  was  to 
remove  impurities  from  the  probes  operating  surface.  How¬ 
ever,  it  was  also  observed  that  if  the  probe  was  submerged 
and  set  to  operate,  the  recorded  voltages  could  only  be 
guaranteed  for  a  period  of  some  20  seconds,  whereby  the 
probe3  sensing  characteristics  would  start  to  reduce 
thereby  giving  slowly  decreasing  voltage  results.  The 
wandering  of  the  probes  voltage  would  eventually  end  up 
at  zero  after  about  45-60  seconds.  Therefore,  it  was  nec¬ 
essary  to  limit  the  test  time  to  approx.  10  seconds, 
whereby  30  voltage  readings  v/ere  recorded.  If  the  probe 
was  now  withdrawn  from  the  flow  and  brushed  with  clean 
water,  re-submerged  and  switched  to  operate,  the  probe 
would  once  again  be  operating  satisfactorily.  Therefore, 
it  is  quite  evident  that  some  form  of  insulation  was 
building  up  on  the  probe,  which  was  most  probably  due  to 
air  bubbles  forming  on  the  surface  caused  by  nucleation 
effects  due  to  the  operating  temperature  of  the  probe. 

The  perpetual  brushing  of  the  delicate  probes  sensing 
element  inevitably  meant  that  the  odds  against  damaging 
the  probe  were  reducing.  Probe  1  was  however,  damaged 

135 


J 


in  this  way.  Before  Probe  2  could  be  used  this  meant  re¬ 
setting  of  the  probes  operating  resistance  ana  frequency 
response  \v£3  necessary,  and  then  recallibraticn  Thus  re¬ 
ducing  the  overall  efficiency  in  the  method  of  velocity 
measurements.  As  mentioned  in  section  2*3,  the  system 
incorporates  a  filtration  system. This  system  cannot, 
however,  be  guaranteed  to  filter  all  the  water  in  the 
system  and  particles  of  rust  and  debris  still  pass  throug 
the  channel.  Probe  2  was  damaged  at  probe  traverse  (2.3) 
by  one  such  deposit  hitting  the  probe. 

The  calibration  of  the  probes  was  carried  out  in  a  glass 
tube,  where  the  flow  was  such  that  the  Re  numbers  indi¬ 
cated  fully  turbulent  flow.  If  however,  the  probes  had 
then  been  used  to  obtain  velocity  measurements  whereby 
the  Re  numbers  indicated  laminar  flow  the  equations 
(4.8)  and  (4.9)  could  not  have  been  used  with  accuracy. 
Therefore  this  presents  another  problem  using  the  CTA, 
and  that  is  that  it  must  only  be  used  in  the  same  type 
of  flow  it  is  calibrated  in  ,  otherwise  results  would  be 
suspect . 

Therefore,  from  these  observations  and  the  difficulties 
encountered,  it  would  indicate  that  the  hot  film  probe 
is  not  very  suitable  a  method  for  velocity  measurements 
in  a  substance  that  has  such  large  conductivity  changes, 
and  varying  flow  regimes  unless  calibration  in  the  flew 
types  has  been  carried  out. 
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10.3  Velocity  Profiles. 

The  graphs  as  plotted  for  probe  traverses  (1.1),  (2.1), 
(2.2)  and  (2.3)  do  however  give  fairly  good  interpreted 
velocity  profiles.  To  establish  the  thickness  of  the 
boundary  layer  this  has  been  taken  as  to  extend  to  38$ 
of  the  freestream  velocity,  and  for  probe  traverses 

(2.1)  and  (2.2)  was  52mm  and  65mm  respectively. 


10.4  Evaluation  of  C^. 

Initially  was  calculated  from  the  Ludwieg  -  Tillman 
empirical  equation  (8.11),  but  first  it  was  necessary 
to  calculate  the  displacement  thickness  £  ,  and  the 

momentum  thickness  ©  .  For  the  calculation  of  6  and  © 
the  Prandtl  power  law  was  used  wherby  u/U  =  (y/<S  )n 
where  a  is  some  power  (usually  1/7  for  smooth  surfaces). 

The  results  table  SI  and  E2  are  for  probe  traverses 

(2.1)  and  (2.2)  respectively  and  indicates  that  the 
best  values  of  a  which  satisfy  (y/<5  )"  =  u/J. 

For  rough  surfaces  the  value  of  a  reduces  to  less  than 
1/7  and  these  results  illustrate  that  this  is  so,  the 
value  of  a  varies  from  l/4  to  1/5.6  for  the  profiles 

(2.1)  and  (2.2). 
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Once  the  shape  of  the  velocity  profile  was  obtained  the 
equations  (8.4)  and  (8.6)  could  be  evaluated  and  Cf 
obtained  directly.  The  calculations  of  which  are  in 
sections  8.7. 

The  equation  (8.11)  is  an  empirical  equation  and  is 
formerly  accepted  as  an  addition  to  boundary  layer 
computation.  However,  its  accuracy  is  only  of  the  order 
of  t  10;* *  due  to  the  fact  that  the  pov/er  law  dees  not 
give  a  good  approximation.  However,  from  these  calcul¬ 
ations,  one  power  law  was  not  used  and  the  values  of  © 
and  6  were  calculated  from  an  addition  cf  power  laws 
that  gave  the  best  profile,  therefore  the  accuracy  of 
these  results  should  be  well  within  t  1($. 

The  next  method  of  calculating  Cf  was  using  the  log¬ 
arithmic  inner  law.  By  plotting  graphs  FIG.  21  and  22 
of  velocity  Vs  logc  y  we  can  see  that  there  is  a  small 
region  where  the  graph  is  straight.  However,  this 
particular  point  is  not  so  well  defined. 

Equation  (8.27)  is  the  inner  law  equation  for  a  smooth 
flat  plate  and  graphs,  FIG.  E3>  E4,  E5  and  E6  appendix 
E  illustrate  the  inner  regions  of  probe  traverse  (2.1) 
and  (2.2).  FIG.  21  and  E2  were  plotted  using  different 
values  of  friction  velocity  u«. 
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The  values  of  Cf  obtained  using  the  Ludwieg  -  Tillman 
equation  for  probe  traverses  (2.1)  and  (2.2)  were 
0.00334  ana  0.00323  respectively.  Therefore  by  calcul¬ 
ating  X  from  equation  (8.9),  Ho  the  friction  velocity 
can  be  calculated  from  (8.14),  therefore  an  indication 
as  to  the  range  of  Uo  could  be  obtained  for  traverse 
(2.1)  and  (2.2),  u,  was  calculated  to  be  0.023  and 
0.0227  respe ctively. 

From  FIG.  S3  and  E4  the  value  of  u0  whereby  when  the 
inner  region  velocities  and  distance  y  are  plotted 
u/u*  Vs  log€ (yuo/u ) ,  that  gives  the  best  parallel  line 
to  that  for  a  smooth  flat  plate  is  u0  =  0.0241  ana 
u„  =  0.024  for  traverse  (2.1)  and  (2.2)  respectively. 

The  resulting  calculated  values  of  Cf  are  therefore, 
0.00364  and  0.00361  as  noted  in  section  (8.11). 

These  values,  obtained  using  the  inner  region  do  how¬ 
ever  differ  approximately  10#  from  those  obtained  from 
the  ludwieg  -  Tillman  equation,  but  the  two  sets  of 
results  are  however  comparable. 

FIG.  E5  and  E6  indicate  3  sets  of  plotted  points  with 
u«  varying  from  0.022  to  0.026.  Each  profile  plotted 
is  for  y  varying  from  0  to  S  . 
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Prom  these  plotted  profiles  it  is  easier  to  distinguish 
the  inner  and  outer  region  showing  the  overlap  region 
inbetween . 

The  profiles  plotted  with  U®  =  0.0241  for  (2.1)  and 
Uo  =  C.024  for  (2.2)  are  however ,  displaced  downward 
and  amount  as  which  is  as  suggested  by  Scholz  (1S55) 
This  dov/nward  shift  increases  for  larger  values  of 
roughness  k. 

Coles  modified  his  original  smooth  plate  equation  to 
take  account  of  this  downward  shift  whereby  A6  was 
called  the  roughness  function.  This  roughness  function 
can  only  be  obtained  from  experiment  and  modifies  the 
equation  for  the  inner  law  to  that  of  (8.38)  and  (8.39) 
for  (2.1)  and  (2.2)  respectively. 

The  value  of  A8  ,  the  roughness  function  has  also  been 
(22^ 

designated  where  Au.  is  the  difference  in  velocity 

between  a  rough  and  smooth  surface,  at  the  surface,  due 
to  the  roughness  profile. 

Sad  the  surface  been  smooth  then  h>6  or  (A«yu*)  would 
have  been  zero. 
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The  wake  function  has  been  ignored  due  to  it  having 
very  little  effect  because  of  the  closeness  of  the 
inner  region  to  the  surface  in  question. 

The  inner  region  extends  between  36  <  (yu*/i))  <  90  and 
the  logarithmic  overlap  would  seem  to  extend  from  90 < 

( yud /v  )  <  270. 

Various  interpretations  of  this  region  have  been  used 
and  Hiller  vised  the  inner  law  for  a  flat  plate  Eq.(8.27) 
and  using  the  region  100  4  (yu«,A>)  <  30C  for  his  cal¬ 
culation.  He  then  wrote  a  computer  program  that,  for 
each  experimental  point  within  that  region  computed  a 
value  of  u  that  gave  least  residual  to  the  equation 
(8.27).  Results  for  u  from  each  experimental  point  in 
the  defined  range  100  -  300  were  then  averaged  to  give 
their  arithmetic  mean  value. 

In  these  investigations  it  would  seem  that  the  region 
100  -  300  is  outside  the  inner  region,  but  does  corres¬ 
pond  to  the  logarithmic  overlap.  Prom  PIG.  E3,  E4,  E5 
and  E6  we  can  see  that  the  profiles  plotted  for  increas¬ 
ing  u0  rotate  in  a  clockwise  manner.  Therefore,  for  the 
overlap  region  to  be  parallel  to  that  for  a  smooth  plate 
would  have  resulted  in  u.  being  in  the  order  of  0.032 
which  would  have  increased  Of  to  0.006,  whereby  this 
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value  is  approximately  IOO7S  greater  than  that  calcul¬ 
ated  from  the  Ludwieg  -  Tillman  equation 
Therefore  it  must  be  assumed  that  the  inner  region, 
shown  in  FIG.  25  and  26  is  tnat  relating  to  this  part¬ 
icular  profile. 

« 

The  average  rougutiess  height,  k,  has  been  calculated 

(3) 

using  the  Prandtl  -  Schlichting  approx,  solution  and 
was  found  to  be  O.OGC355m  and  0. 00043m  for  probe  trav¬ 
erse  (2.1)  and  (2.2)  respectively. 

These  results  could  however  be  verified  using  surface 
measurement  techniques. 

One  difficulty  in  obtaining  surface  friction  coefficients 
is  that  the  surface  is  not  a  constant  roughness.  There¬ 
fore  due  to  the  flow  regime  being  that  of  tranquil  mode 
the  flow  as  it  passes  over  the  surface  is  influenced 
by  the  flow  that  is  in  front  of  it.  Therefore,  it  would 
be  more  beneficial  if  the  surface  in  question  was  that 
of  a  constant  roughness  which  would  ensure  that  effects 
due  to  uneven  roughness  were  alleviated. 

One  of  the  major  problems  associated  with  the  LBA  is 
that  equation  (7«10)  depends  upon  the  wavelength  of 
laser  light  in  the  flowing  medium.  This  however,  is 
governed  by  the  refractive  index  which  is  temperature 
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dependant.  Therefore,  if  correct  velocity  aeasurements 
arc.-  to  he  obtained,  xhen  it  must  be  ensured  that  the 
correct  value  of  wavelength  is  calculated  at  the  correct 
water  teaperature. 
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CONCLUSION 


11.1  Conclusion. 


One  of  the  major  problems  in  obtaining  the  surface  fric¬ 
tion  coefficient  during  this  project,  was  accurately 
recording  the  velocity  in  the  boundary  layer.  The  prob¬ 
lems  and  disadvantages  have  been  discussed  in  the  pre¬ 
ceding  sections  in  detail.  However,  had  the  time  been 
available  then  these  could  have  been  checked  using  the 
Laser  Doppler  Anemometer. 

Evaluation  of  Cf  from  the  inner  la?/  does  require  that 
measurements  of  velocity  are  made  as  accurately  as 
possible. 
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APPENDIX  A 


VOLTAGE  RESULTS  PROSE  (l)  CALIBRATION. 

(1) 

+0275  +C273  +0272  +0270  +0268  +0275  +0270  +0273  +0270 

+0271  +0274  +0273  +0271  +0268  +0264  +0271  +0272  +0270 

+0269  +0264  +0258  +0263  +0260  +0259  +0258  +0259  +026] 

+0262  +0260  +0256  +0253  +0253  +0254  +0249  +0256  +0252 

+0249  +0255  +0254  +0255 

(Average  Value  =  2.616  Volts.) 

+0202  +0269  +0267  +0273  +0271  +0272  +0273  +0271  +0266 

+0266  +0268  +0268  +0267  +0269  +0268  +0264  +0269  +0268 

+0274  +0272  +0268  +0264  +0264  +0269  +0272  +0269  +0266 

+0266  +0268  +0267  +0267  +0269  +0269  +0269  +0268  +0271 

+0270  +0269  +0264  +0260 
(Average  Value  =  2.6665  Volts.) 

+0064  +0266  +0266  +0264  +0262  +0265  +0265  +0267  +0262 

+0260  +026i  +0264  +0264  +0265  +0263  +026l  +0263  +0265 

+0265  +0262  +0263  +0262  +0267  +0262  +0262  +026i  +0260 

+0261  +0260  +0265  +0267  +0263  +026i  +026i  +0259  +0261 
+0261  +0261  +0263  +00 

(Average  Value  =  2.6289  Volts.) 


(2) 

+0226  +0229  +0225  +0227  +0223  +0226  +0228  +0227  +0227 

+0221  +0219  +0224  +0223  +0227  +0224  +0225  +0225  +0221 

+0227  +0227  +0220  +0223  +0226  +0223  +0222  +0225  +0226 

+0221  +0223  +0224  +0225  +0201  +0201  +0199  +0097  +0198 
+0199  +0195  +0197  +0200  +0196 
(Average  Value  =  2.1863  Volts.) 

+0201  +0231  +0227  +0228  +0223  +0228  +0228  +0226  +0227 
+0228  +0228  +0226  +0226  +0226  +0224  +0229  +0227  +0220 
+0226  +0223  +0224  +0228  +0226  +0222  +0228  +0227  +0224 

+0228  +0226  +0226  +0222  +0224  +0225  +0223  +0221  +0222 

+0218  +0217  +0219  +001 

(Average  Value  =  2.2441  Volts.) 

+022  +0228  +0224  +0224  +0225  +0230  +0228  +0220  +0222 
+0226  +0224  +0224  +0224  +0221  +0225  +0228  +0224  +0228 
+0224  +0217  +0223  +0220  +0223  +0222  +0228  +0225  +0223 
+0223  +0224  +0223  +0225  +0225  +0227  +0223  +0225  +0224 
+0223  +0225  +000 
(Average  Value  =  2.2424  Volts.) 
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(3) 


+0005 
+0179 
+01 82 
+01 82  +0178 
+0178  +0179 


+0183  +01 83  +01 83 
+0l8i  +0183  +0176 
+01 81  +0179  +0179 
+01 81  +0178 
+0178  +0001 


+0178  +0177 
+0181  +0180 
+0182  +0182 
+0179  +0175 


(Average  Value  =  1.7503  Volts.) 


+0178  +0180  +0179 
+oi 80  +0183  +0180 

+0178  +0182  +0l8l 
+0l8l  +0176  +0179 


+0007  +0186  +O186  +O188  +0i 85  +0186  +0189  +0185  +O183 

+0189  +O186  +0189  +0183  +0189  +0186  +0186  +0183  +01 84 

+0186  +0187  +0185  +0187  +O187  +0182  +O183  +0l82  +0181 

+0187  +O187  +0136  +O183  +01 84  +0184  +0188  +0181  +0186 
+0186  +0181  +0181  +0180 
(Average  Value  =  1.8505  Volts.) 

+0186  +Ol87  +O186  +O189  +O186  +0186  +O186  +0186  +O186 

+0187  +0187  +0187  +0188  +01 86  +0185  +0187  +0187  +0186 

+0187  +0189  +0188  +0190  +0190  +0186  +0186  +0188  +0192 
+0187  +0188  +0188  +0186  +0186  +Oi87  +0188  +0188  +0185 
+0183  +0182  +O186  +0001 

(Average  Value  =  1.869  Volts.) 


(4) 

+0139  +01 40  +0142  +0137  +0l4i  +01 4i  +01 4i  +01 4i  +0i44 
+01 4i  +01 4o  +01 42  +0137  +01 4o  +0i40  +01 39  +0142  +0138 

+01 4i  +0137  +01 41  +0137  +0139  +0142  +01 4o  +0139  +0135 

+0137  +0142  +0135  +0133  +0135  +0142  +0138  +0139  +01 4i 

+01 4o  +0142  +0137  +0100 

(Average  Value  =  1.374  Volts.) 

+0043  +0142  +01 44  +01 4i  +O138  +01 4o  +0139  +0136  +0i37 

+01 4i  +0137  +0137  +0136  +0133  +0139  +0135  +0136  +0137 

+0137  +0129  +0138  +0138  +0131  +0140  +0133  +0134  +0135 

+0139  +0137  +0138  +0140  +0138  +0138  +0137  +0136  +0136 
+0138  +0137  +0138  +0001 
(Average  Value  =  1.3737  Volts.) 

+0137  +0137  +0139  +01 4i  +0138  +0136  +0138  +0139  +0136 

+0136  +0137  +0137  +0140  +0138  +0138  +0136  +0134  +0138 

+0137  +0136  +0139  +0138  +0137  +0139  +01 4o  +0139  +01 38 
+0138  +0138+0135  +0139  +0139  +0138  +0137  +0137  +0138 
+0136  +0137  +0138  +0000 

(Average  Value  =  1.3764  Volts.) 
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(5) 

+2080  +0080  +0081  +0076  +0079  +0080  +0080  +0080  +0077 
+0078  +0078  +0080  +0080  +0079  +0079  +0079  +0078  +0077 

+0083  +0079  +0081  +0080  +0078  +0080  +0079  +0082  +0030 

+0078  +0080  +0078  +0078  +0081  +0078  +0077  +0080  +0078 
+0076  +0077  +0082  +000 
(Average  Value  =  0.791  Volts.) 

+006  +0078  +0080  +OO81  +0076  +0080  +0078  +003l  +0078 
+0080  +0078  +0075  +008 1  +0082  +0079  +0079  +0079  +0078 
+0079  +0080  +0078  +0073  +0079  +0077  +0075  +0078  +0079 

+0077  +0078  +0078  +0080  +0078  +0074  +0077  +0079  +0078 

+0080  +0080  +0078  +0000 

(Average  Value  =  0.785  Volts.) 

+000  +008 1  +0080  +0077  +0082  +0078  +0081  +0076  +0078 

+0081  +0078  +0073  +0078  +0078  +0076  +0078  +0080  +0079 

+0079  +0079  +0077  +0078  +0077  +0076  +0079  +0080  +0079 

+0079  +0078  +0076  +0075  +0077  +0078  +0079  +0079  +0078 
+0077  +OO80  +0079  +20 
(Average  Value  =  0.784  Volts.) 


(6) 

+00^9  +0051  +0050  +0050  +0051  +0049  +0049  +0049  +0051 

+0047  +0048  +0052  +3050  +0050  +0052  +0053  +0050  +0049 

+0051  +0050  +0048  +0046  +0049  +0047  +0050  +0051  +0052 

+0050  +0048  +080 
(Average  Value  =  0.497  Volts.) 

+0040  +0047  +0047  +0047  +0049  +0049  +0051  +0047  +0049 
+0050  +0048  +0048  +0048  +0052  +0050  +0050  +0049  +0052 
+0049  +0050  +0050  +0048  +0049  +0050  +0053  +3048  +0052 

+0049  +0046  +0048  +0050  +0048  +0053  +0047  +0049  +0043 

+0049  +0048  +0049  +0051 
(Average  Value  =  0.489  Volts.) 

+3849  +0048  +0047  +004q  +0050  +0049  +0047  +0046  +0048 

+0049  +0049  +0049  +0050  +0050  +0051  +0048  +0049  +0049 

+0051  +0050  +0049  +0048  +0049  +0049  +0049  +0051  +0048 

+0049  +0049  +0049  +00494+0047  +0050  +0050  +0050  +0051 

+0049  +0049  +0049  +0800 

(Average  Value  =  0.490  Volts.) 
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CALCULATIONS  FOR  CALIBRATION  OF  PROBE  1 


SET  (1)  (1)  29^0  =  2?0>642  cc  /B 

iirr  =  266*667  cc  /s 

^  xiTT  =  269-^9  cc  /a 
Average  value  of  (l),(2),(3)  =  268.893  cc  /a 

Diameter  cf  calibration  tube=  2.15  cm 

=  0.0215  m 

Area  =  AD* =3. 1415927x2. 15x2. 15 

T - ? - 

=  3.631  cm* 

Average  mean  velocity  =  A^ra ig.e ,  vaXp e  of,  flowrate 

268.893 

"  3. 631" 

=  74.055  cm/a 

*  0.74055  m/a 

Reynolda  No.  (Re)  a 

Where  1  =  diameter  of  calibration  tube. 

From  fig.  4.3  at  8°C  t>=  0.0000015  ma/s 

b6  Va x  0.0 21 5  _  i  rt 

Be  “  O.OCCoOl!?  14333.333  va 
From  above  v»  «  0.74055  m/s 
Re  *  14 333.333  x  0.74055  =  10614.55 
From  fig.  4.2  with  Re=  10614.55 

Eatio  of  Average  Velocity  ,  0.7915 
Centre-Line-Velocity 

Centre-Line-Velocity  =  0'1t915^  =  0.936  m/a 


SET  (2) 


270.642  cc  /s 


(1)  2848  = 

TS7J 

(2)  2485  =  227.982  cc  /s 

IGT? 

(3)  1765  =  229.221  cc  /a 

”777 

Average  value  of  (l),(2)  and  (3)  =  229*582  cm 

Average  mean  velocity  va  =  229*582 

3*631 

=  63*228  cm  /s 

va  =  0.63228  m/s 

Reynolds  Humber  (Re)  =  14333.333  va 

=  9062.74 

Prom  FIG.  4*2  with  Re  =  9062*74 

ratio  of  va/vc  =  0,7865 

Centre-line  velocity  vc  =  0.63228 


SET  (3) 

(1) 

1555  = 
5TT7 

194*375 

cc 

/a 

(2) 

1795  = 

9*2 

195*109 

cc 

/s 

(3) 

1590  = 

*"57T 

196.296 

cc 

/s 

Average  value  of  (1),(2)  and  (3)  *  195*26  cm3 

Average  deans  velocity  va  =  195*26 

3*631 

=  53*776  cm  /s 

va  =  0* 53776  m/s 

Reynolds  Number  (Re)  =  14333*333  x  va 

=  7707*87 

Prom  PIG.  4*2  with  Re  =  7707*87 

ratio  of  va/vc  =  0.7835 


Centre-line  velocity  vc 


0.53776 

”7777535 


SET  (4) 

(1) 

1595 

TUT? 

= 

156.373 

cc 

/s 

(2) 

1865 

I2TU 

= 

155.417 

cc 

/s 

(3) 

1714 

11.1 

= 

154.415 

cc 

/s 

Average  value 

of  (1), 

(2)  and 

(3) 

=  155 

.41 

cm 

Average  means  velocity  va  =  155»41 

T7t>31 


=  42.798  cm  /s 

va  =  0.42798  m/s 


SET  (5) 


(D  1455 

TT72 


102.465  cc  /s 


(2)  1555  =  102.980  cc  /s 

I5TT 

(3)  1490  =  103.472  cc  /s 

14. 4 

Average  value  of  (l)f(2)  and  (3)  =  102.973  cm3/s 

Average  mean  velocity  ya  =  102.973 

'■17631“ 

=  28.359  cm  /s 

ya  =  0,28359  m/a 

Reynolds  Number  (Re)  =  14333.333  x  va 

=  4064.824 

Prom  PIG.  4.1  with  Re  =  4064.824 

ratio  of  va/vc  =  0.745 

Centre-line  velocity  vc  =  0. 28359 

0.745" 


vc  s  0.3807  m/a 
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SET  (6) 


(1)  1280 

1777 


72.316  cc  /s 


(2)  1320  =  71.739  cc  /s 

TH7T 

(3)  1400  =  71795  cc  /s 

1975 


Average  value  of  (l),(2)  and  (3)  =  71*95  cm/s 

Average  mean  velocity  va  =  71.95 

T3il 


=  19.816 


va  a  0.19816  m/3 


Reynolds  Number  (Re)  =  14333*333  x  va 


=  2840.22 


From  FIG.  4.1  with  Re  *  2840.22 


ratio  of  va/vc  =  0.7025 


Centre-line  velocity  vc  =  0.19816 

"077025 


=  0.282  m/s 
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CALCULATIONS  FOR  PROBB  1  CALIBRATION  GRAPH 


Using  the  Method  of  Least  Squares  as  described  in  section 
4.8  and  using  equations  (4«5)>  (4.6)  and  (4.7). 


X 

x1 

Y 

a 

Y 

XY 

0.9360 

0.87609600 

2.6371 

6.95429610 

2.46832560 

0.8039 

0.64625521 

2.2243 

4.94751049 

1.78811477 

0.6864 

0.47114496 

1.8233 

3.32442289 

1.25151312 

0.5598 

0.31337604 

1.3747 

1.88980009 

0.76955106 

0.3807 

0.14493249 

0.7870 

0.61936900 

0.29961090 

0.2820 

0.07952400 

0.4920 

0.24206400 

0.13874400 

3.6488 

2.53132870 

9.3384 

17.9774 6288 

6.71586545 

6a  +  3.6488b  =  9*3384  (4.5.1) 

3.6488a  +  2.5313287b  =  6.71586545  (4.6.1) 

(4.5.1)  x  3-6488  AND  (4.6.1)  x  6 
21.8928a  +  13.31374144b  =  34.07395392  (4.5.2) 

21.8928a  +  15.18797220b  =  40.2951927  (4.6.2) 


(4.6.2)  -  (4.5.2) 


b  =  3.319355819 
b  =  3.319356 


Substitute  b  into  equation  (4.5.1) 

6a  =  9.3384  -  3.6488  x  3.319356 

a  =  -0.462211 


Equation  for  the  best  straight  line  is: 

V  =  3.319356v  -  0.462211 
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therefore 


v  =  0.301263V  +  0.1392472 

where  v  =  velocity 
and  V  =  voltage. 


To  check  whether  or  not  statistically  the  equation  gives 
a  good  fit  to  the  experimental  data: 


(4.7)  = 


nf.  XY  -  SXtY 


J(  mix’-  (sxfn  rt? -  (Vi>Y 


=  0.999793 

Therefore  the  above  equation  does  give  a  good  fit  to  the 
experimental  data. 


■BESIiSTBAI&H2i:LIK2 


>=  0.301263V  +  0.1392472 


n 


VOLTAGE  RESULTS  PROBE  (2)  CALIBRATION. 


(A) 

+004  +0334  +0335  +0333  +0336  +0331  +0331  +0327  +0328 
+0329  +0329  +0327  +0331  +0334  +0328  +0327  +0327  +0330 
+0330  +0329  +0331  +0326  +0325  +0326  +0326  +0329  +0328 

+0328  +0329  +001 


+0334  +0335  +0336  +0335  +0333  +0330  +0332  +0335  +0333 

+0331  +0333  +0332  +0332  +0334  +0334  +0331  +0332  +0330 
+0332  +0331  +0327  +0328  +0325  +0329  +0328  +0329  +0328 
+0327  +0330  +0328  +000 


+0329  +0329  +0327  +0325  +0321  +0326  +0319  +0319  +0319 
+0317  +0319  +0319  +0318  +0316  +0309  +0308  +0309  +0307 
+0303  +0302  +0305  +0305  +0313  +0319  +0320  +0318  +0319 
+0315  +0315  +0316  +001 
(Average  Value  =  3*255  Volts.) 


(B) 

+0208  +0204  +0205  +0206  +0203  +0205  +0208  +0202  +0207 
+0205  +0204  +0202  +0205  +0204  +0202  +0202  +0204  +0199 
+0200  +0202  +0197  +0198  +0199  +0202  +0195  +0194  +0198 
+0201  +0199  +0000 


+0210  +0211  +0205  +0210  +0210  +0211  +0204  +0208  +0203 

+0204  +0206  +0203  +0206  +0205  +0204  +0209  +0207  +0205 

+0205  +0207  +0206  +0206  +0206  +0201  +0204  +0202  +0204 

+0198  +0201  +0200 


+000  +0204  +0208  +0205  +0200  +0208  +0207  +0207  +0205 
+0207  +0204  +0206  +0202  +0207  +0204  +0203  +0204  +0204 
+0207  +0206  +0205  +0204  +0205  +0207  +0204  +0201  +0206 
+0205  +0205  +0000 
(Average  Value  =  2.042  Volts.) 
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(C) 

+015  +0153  +0156  +0152  +0 1 5 1  +014$  -K5] 49  +0l47  +0150 
+01 49  +0150  +01^9  +0150  +0151  +01  49  +0148  +0146  +Ol46 
+0147  +0149  +0149  +01 51  +0152  +0150  +0152  +0152  +01 46 
+0146  +0147  +0001 


+0005  +0154  +0152  +0150  +0150  +01 5 1  +0152  +0156  +0150 
+0154  +0151  +01 48  +0147  +0145  +0151  +0155  +01 48  +01 46 
+01 48  +0154  +0150  +0150  +01 46  +01 50  +0146  +01 47  +01 47 
+0146  +0147  +000 


+0148  +0152  +0150  +0150  +0151  +0151  +0149  +0150  +01 47 
+0150  +0148  +0150  +0145  +01 48  +0147  +Ol46  +0151  +01 49 
+0150  +0149  +0152  +0150  +01 48  +01 48  +01 49  +0152  +01 44 
+0147  +0150  +0000 
(Average  Value  =  1.494  Volts.) 


(D) 

+0127  +0124  +0127  +0122  +0112  +0113  +0111  +0109  +0109 
+0111  +0111  +0113  +0113  +0109  +0112  +0111  +0112  +0111 
+0114  +0I08  +0110  +0113  +0114  +0110  +0113  +0111  +0109 

+0113  +0113  +0100 


+0124  +0124  +0124  +0123  +0120  +0123  +0123  +0120  +0120 
+0123  +0127  +0121  +0120  +0121  +0120  +0123  +0119  +011 8 
+0124  +0118  +0121  +0121  +0117  +0120  +0121  +0123  +0123 
+0122  +0119  +002 


+0123  +0119  +0126  +0126  +0121  +0127  +0124  +0124  +0125 
+0126  +0124  +0120  +0122  +0122  +0119  +0120  +0124  +0121 
+0120  +0123  +0126  +0127  +0123  +0123  +0122  +0122  +0126 
+0123  +0121  +0100 
(Average  Value  =  1.188  Volts.) 
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(E) 

+208  +0080  +0083  +008 1  +008 1  +0082  +0083  +0078  +0080 
+0084  +0084  +0079  +0083  +0079  +0082  +0083  +0082  +0078 
+008l  +0081  +0080  +0079  +0077  +008l  +0084  +0082  +008l 
+0080  +0078  +000 


+0080  +0078  +0080  +0080  +0082  +0079  +0080  +0082  +0082 

+0080  +0079  +0080  +0079  +0078  +0077  +0079  +008 1  +0073 

+0081  +0080  +0080  +OO81  +0079  +0081  +0082  +008 l  +0078 

+0083  +0080  +001 


+2083  +0083  +0083  +0084  +0084  +0082  +0083  +0078  +0082 

+OO81  +0080  +0080  +0082  +0079  +0083  +0079  +0081  +0082 

+0082  +OO81  +0079  +0080  +0080  +0079  +0079  +0081  +0082 

+0080  +OO81  +001 
(Average  Value  =  0.807  Volts.) 

(P) 

+0806  +0037  +0036  +00354+0036  +0037  +0036  +0038  +0035 
+0033  +0037  +0034  +0036  +0036  +0036  +0038  +0036  +0035 
+0036  +0035  +0034  +0036  +0035  +0036  +0037  +0037  +0034 
+0037  +0036  +0300 


+003  +0038  +0035  +0038  +0034  +0036  +0035  +0037  +0037 
+0035  +0038  +0035  +0036  +0035  +0035  +0039  +0036  +0035 
+0036  +0036  +0036  +0034  +0033  +0035  +0036  +0036  +0034 
+0036  +0034  +0030 


+0037  +0035  +0036  +0037  +0036  +0034  +0035  +0035  +0039 

+0035  +0036  +0037  +0034  +0035  +0036  +0034  +0036  +0037 

+0036  +0035  +0035  +0037  +0035  +0038  +0035  +0037  +0034 

+0035  +0036  +080 
(Average  Value  =  0.356  Volts.) 
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BOUNDARY  LAYER  PROBE  (2)  CALIBRATION 


POSITION 

OP  GLOBE 

VALVE 

PLOW  RATE 

cc 

TILTS 

sec 

AVERAGE 

FLOW  RATS 

cc  /s 

HEAR 

VELOCITY 

m/s 

2570 

7.6 

A 

2852 

8.2 

342.646 

0.9438 

2770 

8.1 

2210 

8.4 

-  B 

1750 

6.3 

273.112 

0.75227 

1680 

6.5 

8.4 

C 

1  1 

8.5  . 

209.284 

0.57646 

wSBM 

8.4 

1650 

8.75 

D 

1590 

8.65 

186.338 

0.513257 

1605 

8.6 

■ra%n 

■Hi 

..  E 

1 

■ 

141.809 

.0.39061 

m 

890 

10.1 

P 

1155 

13.5 

87.087 

0.23988 

1270 

14.5 

TABLE  A2. 
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CALCULATIONS  POP.  CALIBRATION  OP  PROSE  2. 


SET  (A) 

(1) 

2570 

=  338.158  cc  /s 

T7Z 

(2) 

2852 

=  347.805  cc  /s 

• 

8.2 

(3) 

2770 

=  341.975  cc  /s 

-87T 

Average 

value  of  (1) 

,(2)  and 

(3)  =  342.646 

Area  of 

calibration 

tube  = 

3.631  cm 

Average  mean  velocity  =  342.646 

'3.631 

=  94.38  cm  /s 

»  0.9438  m/s 

Reynolds  No.  (Re)  =  pVa.1  =  Ygl 

/L  V 

Prom  PIG.  4.3  at  8*  C  V  =  0.0000015  mVs 

Re  =  14333-333  Va 

Prom  above  Va  =  0.9438  m/s 

Re  =  13527*79 

Prom  PIG.  4.2  with  Re  =  13527.79 

Ratio  of  VaAc  ^  0.795 

Centre-line  velocity  =  0.9438  =  1.187  m/s 

“0.795 

Calculations  for  (B) , (C) , (L) , (E)  and  (P)  are  similar  to 
the  above, results  of  which  are  in  table  A4. 
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AD-AOfll  394  TEESSIDE  POLYTECHNIC  MIDDLESBROUSH  (ENGLAND)  DEPT  0— ETC  F/6  13/10 
DETERMINATION  OF  SURFACE  FRICTION  COEFFICIENT  of  A  REPLICA  OF  S— ETC(U) 
APR  79  D  C  AUSUSTSON  N0001 4*70-9-0059 

UNCLASSIFIED  NL 


CALCULATIONS  FOR  PROBE  2  CALIBRATION  GRAPH. 


Using  the  Kstb.od  of  Least  Squares  as  described  in  section 
4.8  and  using  equations  (4.5),  (4.6)  and  (4.7). 


X 

2 

X 

Y 

* 

Y 

XY 

1.1870 

1.40896900 

3.255 

10.595025 

3.8636850 

0.9490 

0.90060100 

2.042 

4.169764 

1.9378580 

0.7329 

0.53714241 

1.494 

2.232036 

1.0949526 

0.6567 

0.43125489 

1.188 

1.411344 

0.7801596 

0.5100 

0.26010000 

0.807 

0.651249 

0.4115700 

0.3268 

0.10679824 

0.356 

0.126736 

0.1163408 

4.3624 

3.64486554 

9.142 

19.186154 

8.2045660 

6a  +  4.3624b  =  9.142  (4.5.3) 

4.3624a  +  3.64486554b  =  8.204566  (4.6.3) 

(4.6.3)  x  4.3624  AND  (4.6.3)  x  6 

26.1744a  +  19.03053376b  =  39.8810608  (4-5.4) 

26.1744a  +  21.86919324b  =  49.227396  (4.6.4) 


(4.6.4)  -  (4.5.4) 

2.83865948b  =  9.3463352 

b  =  3.292517213 
b  =  3.292517 

Substitute  (b)  into  equation  (4.5.3) 

6a  ♦  4.3624  x  3.292517213  =  9.142 

a  =  -0.8702128484 
a  =  -0.870213 
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Equation  for  the  heat  straight  line  ia: 


V  =  3.292517v  -  0.870213 


therefore 


v  =  0.303719V  +  0.2643 

where  v  =  velocity 
and  Y  =  voltage. 


To  check  whether  or  not  statistically  the  equation  gives 
a  good  fit  to  the  experimental  data: 


(4.7)  = 


n  £XY  -  $X1Y 


y(  n$x"-  (SXf  )(  n*Y‘-  (*Y)V 


=  0.9895 


Therefore  the  line  obtained  ia  a  good  fit  to  the 
experimental  data. 
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BESTiSTHA  ISKEi  L  IS 


BOUNDARY  LAYER  PROBE  (1)  CALIBRATION 


CENTRE 


PCS-  AVERAGE  REYNOLDS 
ITION  MEAN  NTIT™ 

OP  VELOCITY  r<U1'^ER 
VALVE 


va/vc 


VELOCITY 


AVERAGE 

VOLTAGE 

READING 


B 

0.74055 

10614.55 

0.7915 

0.936 

2.6371 

2 

0.63228 

9062.74 

0.7865 

0.8039 

2.2243 

3 

0.53776 

7707.87 

0.7835 

0.6864 

1.8233 

B 

0.42798 

6134.45 

0.7645 

0.5598 

1.3747 

5 

0.28359 

4064.82 

0.745 

0.3807 

0.787 

6 

0.19816 

2840.22 

0.7025 

0.282 

0.492 

' 

TABLE  A3. 

BOUNDARY  ] 

IiAYER  PROBE  (2)  CALIBRATION 

POS¬ 

ITION 

OP 

VALVE 

AVERAGE 

MEAN 

VELOCITY 

a/s 

REYNOLDS 

NUMBER 

(He) 

va/vc 

CENTRE 

LINE 

VELOCITY 

m/s 

AVERAGE 

VOLTAGE 

READING 

(V) 

B 

0.9438 

13527.79 

0.795 

1.187 

3.255 

B 

0.75227 

10782.54 

0.792 

0.949 

2.042 

C 

0.57646 

8262.59 

0.7865 

0.7329 

1.494 

B 

0.513257 

7356.68 

0.7815 

0.6567 

1.188 

E 

0.39061 

5598.74 

0.766 

0.51 

0.807 

P 

■in 

0.23988 

3438.28 

0.734 

0.3268 

0.356 

TABLE  A4. 
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VOLTAGE  READINGS  PROBE  T3AVE5SE  (1.1) 


+0072  +0081  +0071  +0063  +0079  +0073  +OO81  +0085  +0068 
+0071  +0069  +0082  +0089  +0069  +0063  +0073  +0072  +OO81 

+0060  +0092+0093  +0070  +0075  +0070  +0071  +0039  +0095 

+0083  +0078  +0088  +0090  +0084  +000 

+0075  +0048  +00 S 7  +0074  +0084  +0058  +0082  +0088  +0080 

+0057  +0077  +0067  +0088  +0077  +0079  +0073  +0083  +0080 

+0084  +0079  +0085  +0070  +0064  +0078  +0067  +0087  +0073 
+0092  +0074 

+0087  +0090  +0082  +0064  +0069  +008l  +0084  +0074  +0086 

+0080  +0083  +0064  +0065  +0087  +0067  +0076  +0072  +0082 

+0087  +0087  +0074  +0065  +0080  +0102  +0091  +0071  +0075 

+0093  +0087  +0079 
(Average  Value  =  0.7741  Volts.) 

+2080  +0098  +0097  +0076  +OO63  +0075  +0087  +0095  +0080 

+0084  +0090  +0102  +0089  +0093  +0078  +0084  +0079  +0075 

+0082  +0079  +0074  +OO81  +0088  +0082  +0078  +0089  +0080 
+0068  +0081  +001 
(Average  Value  =  0.83  Volts.) 

+0085  +0070  +0067  +0070  +0086  +0091  +0093  +0084  +0083 
+0084  +OO81  +0088  +0063  +0087  +OO81  +0096  +0078  +0071 
+0087  +0092  +0085  +0085  +0086  +0081  +0062  +0076  +0073 
+0069  +0082  +001 
(Average  Value  =  0.8072  Volts.) 

+0087  +0087  +0077  +0039  +0084  +0088  +0073  +0090  +0084 
+0079  +0091  +0091  +0087  +0090  +0066  +0092  +0092  +0087 
+0090  +0085  +0087  +0080  +0066  +0073  +0078  +0069  +0082 
+0089  +0078  +0073  +0093  +0090  +0086  +0088  +0075  +0072 
+0085  +0086  +0093  +0002 
(Average  Value  =  0.8351  Volts.) 

+0070  000+0007  +0088  +0030  +0084  +0079  +0089  +0086  +0092 
+0087  +0091  +0088  +0084  +0073  +0082  +0081  +0067  +0067 
+0085  +0095  +0091  +0091  +0075  +0082  +0083  +0089  +0080 
+0079  +0088  +0079  +0091  +0091  +0097  +0089  +0081  +0093 
+0087  +0085  +000 
(Average  Value  =  0.8430  Volts.) 
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0091  +0085  +0087  +0095  +0085  +0078  +0083  +0097  +0085 
+0090  +0088  +0091  +0096  +0096  +0092  +0087  +0101  +0098 

+0050  +0084  +0087  +0081  +0091  +0082  +0077  +0098  +0089 

+0091  +0086  +0092  +0000 

(Average  Value  =  0.891  Volts.) 

+0070  +0082  +0088  +0083  +OO81  +0086  +0079  +0094  +0086 

+0079  +0086  +0087  +0079  +0065  +0085  +0084  +0075  +0094 

+0090  +0087  +0076  +0087  +0088  +0082  +0090  +0088  +0083 

+0088  +0094  +0002 
(Average  Value  =  0.84  Volts.) 

+0077  +0082  +0098  +0094  +0087  +0094  +0096  +0091  +0082 
+0088  +0093  +0088  +0102  +0087  +0093  +0095  +0097  +0093 
+0096  +0092  +0100  +0092  +0083  +0089  +0099  +0094  +0090 
+0096  +0094  +0096  +0001 
(Average  Value  =  0.9194  Volts.) 

+0096  +0094  +0087  +0088  +0091  +0091  +0096  +0099  +0094 
+0097  +0094  +0094  +0097  +0089  +0094  +0092  +0100  +0090 

+0094  +0092  +0092  +0093  +0104  +0096  +0097  +0087  +0085 

+0100  +0088  +0089  +0085  +0086  +0090  +0091  +002 
(Average  Value  *  0.9271  Volts.) 

+2090  +0103  +0102  +0101  +0102  +0096  +0097  +0090  +0092 
+0102  +0102  +0098  +0104  +0102  +0100  +0099  +0095  +0089 
+0091  +0096  +0098  +0102  +0092  +0090  +0083  +0092  +0092 

+0102  +0099  +0092  +0093  +0100  +0098  +0097  +OO89  +0095 

+0091  +0093  +0091  +0000 

(Average  Value  a  O.9590  Volts.) 

+0097  +0100  +0096  +0092  +0089  +0095  +0097  +0092  +0101 

+0099  +0098  +0101  +0097  +0090  +0093  +0095  +0096  +0102 

+0100  +0099  +0088  +0091  +0087  +0098  +0098  +0098  +0097 

+0097  +0098  +0094  +0092  +0098  +0106  +0106  +0099  +2001 
(Average  Value  =  0.9646  Volts.) 

+2000  +0095  +0097  +0106  +0097  +0102  +0094  +0096  +0101 

+0103  +0103  +0094  +0098  +0090  +0096  +0093  +0092  +0086 

+0097  +0090  +0090  +0090  +0084  +0090  +0099  +0093  +0079 

+0091  +0096  +0092  +0095  +0089  +0094  +0092  +0091  +0088 
+0091  +0091  +0090 
(Average  Value  *  0.9356  Volts.) 
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+0109  +0102  +0109  +01 06  +0110  +01 08  +0104  +01 04  +0103 
+0102  +0103  +0103  +0090  +0091  +0094  +0096  +0105  +0112 
+0108  +0106  +0104  +0105  +0101  +0098  +01 06  +0108  +0105 
+0100  +0104  +0102  +0103  +0103  +0102  +0105  +0103  +0105 
+0108  +0105  +0105  +000 
(Average  Value  =  1.0351  Volts.) 

+0108  +0109  +0108  +0112  +0110  +0107  +0109  +01 06  +0103 
+0105  +0105  +0105  .0106  +01 08  +01 08  +0108  +0109  +0110 
+0114  +0112  +0111  +0110  +0113  +0111  +0113  +0112  +0111 
+0111  +0112  +0113  +0112  +0115  +0115  +0115  +0ll6  +0113 
+0113  +0114  +0115  +0000 
(Average  Value  =  1.1044  Volts.) 

+0002  +0105  +0105  +0101  +0100  +0102  +0105  +0112  +0111 
+0109  +0108  +0105  +0109  +0110  +0106  +0105  +0105  +0105 

+0103  +0105  +0108  +0107  +0108  +0105  +0105  +0109  +01 04 

+0106  +0104  +0103  +0105  +0107  +0107  +0106  +0105  +0106 

+0107  +0104  +0102  +0104  +0000 

(Average  Value  =  1.0572  Volts . ) 

+0000  +0109  +0108  +0111  +0107  +0109  +0111  +0108  +0108 
+01 06  +0108  +0108  +0107  +0110  +0108  +0110  +0107  +0111 

+0110  +0113  +0110  +0109  +0107  +0110  +0109  +0109  +0110 

+0114  +0111  +0109  +0107  +0109  +0108  +01 06  +0104  +0108 
+0112  +0109  +0110  +0000 
(Average  Value  =  1.0892  Volts.) 

+0109  +0108  +0108  +01 08  +0109  +01 04  +0102  +0101  +0107 
+0107  +0105  +0109  +0107  +0105  +0103  +0107  +0108  +0107 

+0110  +0111  +0108  +01 08  +0107  +0109  +0107  +01 06  +01 04 

+0107  +0108  +0110  +0109  +0109  +0110  +0100 
(Average  Value  =  1.0697  Volts.) 

+007  +0108  +0113  +0107  +oio6  +oio8  +0107  +0110  +0112 
+0109  +0111  +0109  +0107  +0106  +0107  +0111  +011 4  +0112 
+0109  +0110  +0106  +0107  +0114  +0110  +0109  +0107  +0112 
+0109  +0107  +0107  +0110  +2000 
(Average  Value  =  1.0914  Volts.) 

+0002  +0110  +0111  +0115  +0113  +011 4  +0111  +0109  +011 6 
+0105  +0105  +0108  +0112  +0113  +0119  +0111  +0109  +0108 

+0109  +0106  +0107  +0110  +0115  +0111  +0108  +0109  +0107 

+0109  +0111  +011 4  +0110  +0113  +0108  +0107  +0109  +0108 
+0105  +0111  +0113  +0111 
(Average  Value  *  1.1026  Volts.) 
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+0102  +0099  +0104  +0100  +0108  +0105  +0103  +01 06  +01 06 
+0106  +0101  +0101  +0101  +0104  +0103  +0105  +0102  +0106 
+0112  +0108  +0110  +0107  +0109  +0108  +0108  +0113  +0105 
+0108  +01 06  +0110  +0104  +0110  +0102  +0099  +0104  +0104 
+0102  +0103  +0110  +0110 
(Average  Value  =  1.0535  Volts.) 

+2009  +0107  +0103  +0105  +01 08  +0105  +0101  +0103  +0110 
+0108  +0115  +0110  +0100  +0103  +0100  +01 06  +0111  +0110 

+0105  +0105  +0111  +01 06  +0105  +0109  +01 06  +01 06  +0103 

+0108  +0107  +0105  +0111  +0109  +0112  +0105  +0102  +0105 

+0100 

(Average  Value  =  1.0621  Volts.) 

+0111  +0112  +0120  +0118  +0117  +0113  +011 4  +011 4  +01 1 4 

+0115  +011 6  +0113  +0121  +0118  +0113  +0120  +0120  +0111 

+0108  +0108  +0107  +0109  +0110  +0111  +0112  +0112  +0n5 
+0113  +0110  +0114  +0112  +0112  +0111  +0113  +0111  +0113 
+0110  +0110  +0112  +0100 
(Average  Value  =  1.1283  Volts.) 

+0008  +0132  +0124  +0134  +0130  +0130  +0123  +011 8  +0130 
+01 40  +0133  +0127  +0124  +0l24  +0129  +0127  +0119  +011 6 

+0127  +0121  +0121  +0123  +0129  +01 2l  +0120  +0120  +0125 

+0120  +0119  +0119  +011 6  +011 4  +0119  +0117  +011 6  +011 8 

+0117  +01 l6  +0117  +0100 

(Average  Value  =  1.2566  Volts.) 


PROBE  TRAVERSE  (1.1) 

NOTES  : 

START  0?  TRAVERSE  WATER  TEMP.  12  CENTIGRADE. 

END  OP  TRAVERSE  WATER  TEMP.  12.5  CENTIGRADE. 
CHANNEL  SLOPE  0.3  inches. 

BAIL  PLOW  TUBE  VERTICAL  MANOMETER  READING  3Ft. 

WATER-ON -MERCURY . 
HEIGHT  OP  WATER  IN  CHANNEL  20 5am. 

PROBE  RESISTANCE  READING  6.63  Ohms. 

PRECEDDTG  VOLTAGE  READ  DIGS  TAKEN  AT,  0,1, 2, 3, 4, 

5,6,7,8,9,10,15,20.25,30,35,40,45,50,75,100,150, 
200  mm  ABOVE  MODEL  SURFACE. 

PROBE  POSITION  -  CENTRE  LDIE  MODEL,  OFt  PROM 
LEADD^G  EDGE. 


VOLTAGE  HEADINGS  P305E  TRAVERSE  (2.1) 


+3859  +0064  +0043  +0047  +0052  +0046  +004l  +0045  +0058 
+0067  +0063  +0070  +0051  +0027  +0047  +0047  +0053  +0055 
+0031  +0066  +0034  +0045  +0059  +0022  +006l  +0044  +0051 
+0031  +0061  +000 
(Average  Value  =  0.493  Volts.) 

+0075  +0058  +0075  +0047  +0098  +0077  +0078  +0066  +0101 
+0100  +0092  +0082  +0076  +0081  +0072  +0096  +0092  +0071 
+0109  +0098  +0096  +0097  +0119  +0082  +0070  +0083  +0088 
+0103  +0097  +000 
(Average  Value  =  0.879  Volts.) 

+0094  +0110  +0074  +0098  +0090  +0094  +0100  +0104  +0113 

+0072  +0110  +0081  +0075  +0119  +0080  +0103  +0095  +0084 

+0078  +0118  +0103  +0111  +0099  +0099  +0103  +0089  +0097 

+0109  +0120  +0100 
(Average  Value  =  0.974  Volts.) 

+007  +0103  +0091  +0106  +0107  +0110  +01 11  +0120  +0108 
+0119  +0095  +0107  +0108  +0116  +0103  +0084  +0099  +0125 

+0110  +0125  +0105  +0109  +0121  +0120  +0113  +01 08  +0100 

+0073  +0091  +0000 

(Average  Value  =  1.043  Volts.) 

+0092  +0107  +0116  +0102  +0090  +0085  +0107  +0111  +0120 
+0124  +0092  +0089  +0092  +0108  +0087  +0115  +0093  +0098 
+0117  +0095  +0105  +0115  +0113  +0088  +0115  +0112  +0099 
+0115  +0110  +0100 
(Average  Value  =  1.03733  Volts.) 

+0096  +0117  +01 4i  +0132  +011 6  +0120  +0124  +0130  +0126 
+0111  +0109  +0127  +0095  +0122  +0108  +0109  +0125  +0113 
+0111  +0125  +0117  +0132  +0147  +0118  +0127  +0149  +0133 
,  +0108  +0135  +0100 
(Average  Value  =  1.02077  Volts.) 

+0130  +0126  +0130  +0121  +0126  +0124  +011 6  +0120  +0097 
+0112  +0115  +0130  +0105  +0098  +0111  +0115  +0117  +011 4 
+0100  +0125  +0121  +0108  +0133  +0142  +0134  +0129  +0104 
+0135  +0137  +0000 

(Average  Value  =  1.1886  Volts.) 
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Jj* 


+0132  +0128  +0100  +0095  +0094  +0107  +0105  +0088  +0117 
+0120  +0120  +0123  +01 06  +0110  +011 6  +0129  +01 21  +0129 
+0130  +0121  +0091  +0119  +0134  +0117  +0125  +0124  +0115 
+0120  +0121  +0101 
(Average  Value  =  1.1526  Volts.) 

+0135  +0131  +0120  +0121  +0115  +0111  +O116  +0ll4  +0115 

+009)  +0100  +0115  +0129  +0106  +0105  +0105  +0120  +0117 

+0)1 6  +0)12  +0113  +0088  +0101  +0113  +0127  +0126  +011 6 
+0136  +0117  +0102 
(Average  Value  =  1.1443  Volts.) 

+0133  +0129  +0135  +0126  +0115  +0143  +0102  +0105  +0130 
+0137  +0134  +0129  +011 8  +0109  +0106  +0098  +0129  +0105 

+0120  +0105  +01 08  +0099  +0119  +0108  +0l 1 4  +0105  +0093 

+0119  +0107  +0101 
(Average  Value  =  1.1603  Volts.) 

+0137  +011 4  +011 6  +0128  +0119  +0123  +0120  +0117  +0152 

+0128  +0111  +0124  +0115  +0113  +0121  +0126  +0127  +0125 

+0109  +0107  +0134  +0109  +0126  +0150  +0143  +0i43  +0129 

+0142  +0131  +0130 
(Average  Value  =  1.25633  Volts.) 

+0151  +0138  +0i40  +0132  +0133  +01 58  +0127  +0129  +0132 
+0)48  +0154  +0137  +0143  +0154  +oi 6o  +0156  +0162  +O166 
+0154  +0134  +0149  +01 6l  +0128  +0155  +0142  +Oi47  +Oi43 
+0162  +0155  +0000 
(Average  Value  =  1.4655  Volts.) 

+0158  +0152  +0156  +0160  +0139  +O160  +0150  +0156  +0151 

+0154  +0154  +0134  +0157  +0157  +0l6l  +0146  +0160  +0157 

+0155  +0149  +0153  +0149  +0129  +0136  +Ol64  +0150  +O131 

+0123  +0122  +0000 
(Average  Value  =  1.4365  Volts.) 

+0178  +0156  +0179  +0172  +0168  +0162  +0163  +0174  +0182 

+0158  +0157  +0154  +0172  +01 63  +01 65  +01 67  +01 63  +O162 

+0152  +0182  +01 66  +0155  +0152  +0155  +01 60  +0170  +01 69 
+0173  +0166  +0170 
(Average  Value  =  1.6543  Volts.) 
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+0156  +0l6l  +0171  +0159  +0158  +0157  +0152  +0l64  +0154 
+0161  +0l64  +0166  +0157  ->0154  +0158  +0139  +0130  +0158 

+0156  +0157  +0156  +0148  +0145  +0154  +0144  +0139  +O163 

+O161  +01 44  +016 
(Average  Value  =  1.5486  Volts.) 

+0102  +0170  +O162  +0170  +0167  +0172  +0l64  +0173  +O180 
+0171  +0178  +0166  +O161  +0167  +0158  +O181  +0173  +0171 

+0182  +0177  +0180  +Ol65  +0157  +O166  +0176  +0173  +O171 

,+0177  +0179  +0102 
(Average  Value  =  1.664  Volts.) 

+0183  +0176  +O167  +0175  +0185  +01 88  +O182  +0170  +0178 
+0185  +0175  +0174  +0173  +O180  +0183  +01 84  +0168  +0165 
+0166  +0172  +0162  +0160  +0165  +0167  +0173  +O167  +0166 
+0170  +O167  +O161 
(Average  Value  =  1.729  Volts.) 


+0176  +0177  +0173  +0174  +0165  +0176  +0170  +0178  +0188 
+0177  +0192  +0193  +0188  +0187  +0179  +O182  +0180  +0176 
+0171  +0178  +0166  +0177  +0178  +0183  +0189  +0188  +0183 
+0178  +0174  +0101 

(Average  Value  =  1.7917  Volts-.) 


+010  +0182  +0183  +0178  +0175  +0183  +0187  +0181  +0176 
+0189  +0193  +0185  +0182  +0175  +0177  +0177  -*-0172  +0176 
+0171  +0170  +0174  +0178  +0180  +0175  +0175  +0175  +0175 
+0176  +0172  +0171  +0174 
(Average  Value  =  1.779  Volts.) 


+0177  +0180  +0178  +0180  +0178  +0178  +0176  +0182  +0176 
+0176  +0175  +0177  +0182  +0180  +0174  +0172  +0177  +0181 
+0173  +0172  +0167  +0170  +0176  +0177  +0186  +0180  +0178 
+0179  +0182  +0177  +0174  +017K+0172  +0101 
(Average  Value  =  1.7676  Volrs.) 


+0001  +0180  +0177  +0190  +01 85  +0171  +0183  +0184  +0196 

+0184  +0185  +0187  +0187  +0185  +0187  +0183  +0139  +0185 

+0180  +0133  +0186  +0188  +0137  +0179  +01 34  +0186  +0187 

+0184  +0185  +0001 
(Average  Value  *  1.8454  Volts.) 
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+0101  +0187  +0136  +0187  +0191  +0191  +01 83  +0193  +0190 
+0192  +01 83  +0190  +0196  +0195  +0189  +0194  +01 39  +01 3? 
+0195  +0206  +0197  +0201  +0200  +0195  +0185  +O188  +0193 
+0197  +0001 

(Average  Value  =  1.9185  Volts.) 

+0093  +0192  +0189  +0191  +Ol37  +0192  +0200  +0197  +0204 

+0190  +01 36  +0190  +0186  +0193  +0202  +01 90  +0197  +0194 

+0204  +01 39  +0186  +0186  +0194  +Ol66  +0148  +0146  +0148 

+0153  +0152  +000 
(Average  Value  =  1.84  Volts.) 


NOTES  : 

START  OP  TRAVERSE  WATER  TEMP.  10*  CE11 TIGRADE . 

END  OP  TRAVERSE  WATER  TBIP.  10.5*  CENTIGRADE. 
CHANNEL  SLOPE  0.3 inches. 

DALI  PLOW  TUBE  VERTICAL  MANOMETER  READING  7Ft . 

WATER-ON -MERCURY . 
HEIGHT  OP  WATER  IN  CHANNEL  240mm. 

PROBE  RESISTANCE  READING  6.73  Ohms. 

PRECEDING  VOLTAGE  READINGS  TAKEN  AT,  0,1 ,2, 3, 4, 
5,6,7,8,9,10,15,20,25,30,35,40,45,50,75,100,150, 
200  mm  above  model  surface. 

PROBE  POSITION  -  CENTRE  LINE  MODEL,  2Ft  FROM 
LEADING  EDGE. 
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VOLTAGE  READINGS  PBCBE  TRAVERSE  (2*2) 


+0000  +0093  +0073  +0068  +0070  +0086  +0049  +0057  +0039 
+006l  +0063  +0059  +0050  +0074  +0060  +0056  +0058  +0054 

+0049  +0044  +0052  +0039  +0070  +0055  +0070  +0056  +0043 

+0057  +0044  +000 
(Average  Value  =  O.589  Volts.) 

+0090  +0084  +0079  +0066  +0112  +0089  +0082  +0096  +0101 

+0084  +0083  +0105  +0090  +0087  +0106  +0088  +0104  +0076 

+0079  +0072  +0106  +0071  +0086  +0081  +0084  +0078  +0099 

+0086  +0102  +000 
(Average  Value  =  0.885  Volts.) 

+0002  +0111  +0115  +0088  +0093  +0086  +0103  +0107  +0123 
+0102  +0082  +0094  +0078  +0093  +0088  +0110  +0103  +0086 

+0087  +0086  +0096  +0085  +0113  +0099  +0089  +0083  +0098 

+0073  +0092  +201 
(Average  Value  =  O.948  Volts.) 

+0100  +0125  +0091  +0104  +0083  +0098  +0108  +0094  +0089 

+0115  +0120  +0113  +0121  +0112  +0109  +0117  +0100  +0107 

+0105  +0097  +0102  +0069  +0113  +0108  +0087  +0103  +0112 

+0105  +0097  +0000 
(Average  Value  =  O.999  Volts.) 

+0101  +0106  +0116  +0105  +0098  +0090  +0117  +0115  +0124 

+0111  +0104  +0098  +0119  +0095  +0099  +0112  +0127  +0103 

+0098  +0118  +0106  +0091  +0083  +0114  +0095  +0119  +0117 
+0097  +0091  +0100 
(Average  Value  =  1.056  Volts.) 

+0116  +0076  +0120  +0092  +0121  +0096  +0117  +0122  +0114 
+0088  +0107  +0095  +0088  +0108  +0102  +0097  +0083  +0120 
+0111  +0106  +0103  +0114  +0120  +0115  +0090  +0107  +0102 
+0125  +0132  +01 

(Average  Value  =  1.065  Volts.) 

+0118  +0114  +0112  +0127  +0136  +0122  +0126  +0134  +0113 

+0117  +0115  +0114  +01 40  +0113  +0118  +0104  +0136  +0101 

+0099  +0114  +0104  +0143  +0117  +0103  +0105  +0113  +0126 

+01 40  +0130  +0136 
(Average  Value  =  1.197  Volts.) 
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1  0130  ;  108  +0121  +0084  +0136  +0126  +0126  +0119  +0132 
+0116  +0111  +0124  +011 6  +0117  +0121  +0113  +0119  +0110 

+0106  +0101  +0117  +0115  +0134  +0112  +0121  +0136  +0138 

+0113  +0101 

(Average  Value  =  1.1804  Volts.) 

+0121  +0107  +0117  +0131  +0133  +0128  +0108  +0130  +0107 

+0121  +0091  +0129  +0119  +0101  +0108  +0123  +0134+0131 

+0134  +0147  +0124  +0121  +0ll6  +01 1 6  +0130  +0112  +0104 
+0118  +0100  +000 
(Average  Value  =  1.1935  Volts.) 

+0004  +011 6  +0124  +01 4i  +0137  +01 40  +0129  +0123  +0093 

+0132  +0146  +0147  +0138  +0122  +0117  +0123  +011 6  +0144 

+0127  +0139  +0132  +0117  +0113  +0117  +0133  +0110  +0123 

+0086  +0110  +200 
(Average  Value  =  1.2482  Volt3.) 

+0112  +0122  +01 4i  +0134  +0130  +0093  +0087  +01 46  +0126 
+0107  +0105  +0111  +0136  +0126  +0133  +0115  +0099  +0109 
+01 06  +0115  +0098  +0105  +0122  +0133  +01 4o  +0134  +0123 
+0125  +0148  +0000 
(Average  Value  =  1.2004  Volts.) 

+0154  +0134  +0146  +01 64  +0148  +0152  +0137  +0136  +0145 
+0134  +0146  +0154  +0117  +0136  +0147  +01 4o  +01 64  +0130 
+0136  +0153  +0139  +0136  +0127  +0120  +0125  +0143  +0129 
+0152  +0143  +0101 
(Average  Value  =  1.396  Volts.) 

+0176  +0170  +0154  +0165  +0143  +0158  +01 41  +01 43  +0150 

+0144  +0155  +0163  +0146  +0145  +0137  +01 69  +01 64  +0135 

+0162  +0171  +0158  +0156  +0151  +0131  +0128  +0151  +0146 

+0158  +0159  +0133  +0129  +0100 
(Average  Value  =  1.5133  Volts.) 

+01 41  +0137  +01 49  +0176  +0165  +0169  +0155  +0157  +01 61 

+0156  +0158  +0165  +01 65  +01 42  +0138  +0158  +01 45+01 42 

+01 64  +0159  +0162  +0168  +0168  +0155  +0160  +0165  +0160 

+0152  +0156  +0001 
(Average  Value  =  1.5683  Volts.) 


177 


+0101  +0153  +01 60  +01 64  +0170  +O169  +0158  +0153  +01 60 
+0155  +0148  +0153  +0160  +0160  +0l6l  +0157  +0144  +Ol65 
+0172  +0159  +01 48  +0171  +0157  +0150  +0156  +0169  +0157 
+0153  +0l6l  +0000 

(Average  Value  =  1.5669  Volts.) 

+0173  +0173  +0164  +0176  +0180  +0173  +0 1 61  +0174  +0167 

+0162  +0169  +0146  +0174  +0183  +0175  +0178  +0173  +0175 

+0179  +O163  +0168  +Ol64  +0154  +0l45  +0l7i  +0163  +0170 

+0162  +01 64  +0001 

(Average  Value  =  1.6824  Volts.) 

+0000  +01 64  +01 64  +0165  +0160  +0161  +0159  +0175  +0163 
+0167  +0162  +0172  +0170  +01 64  +01 64  +01 64  +01 54  +0169 
+0167  +0165  +0163  +O163  +0l65  +O166  +Ol56  +0159  +0152 
+0163  +0162  +0000 

(Average  Value  =  1.635  Volts.) 

+O161  +0163  +O161  +0167  +01 64  +0180  +O168  +01 64  +0167 
+0169  +0173  +O166  +0160  +0167  +0i7l  +0174  +0178  +0170 
+0175  +0165  +0171  +0168  +O169  +0173  +O168  +0175  +O172 

+0171  +0178  +0174  +0l82  +O166  +0173  +O161  +0100 

(Average  Value  *  1.7015  Volts.) 


+0101  +0176  +0175  +0179  +0l8l  +01 83  +0180  +0181  +0i85 
+0172+0163  +0170  +0155  +0177  +0170  +0158  +0177  +0179 
+0178  +0168  +0170  +O173  +0176  +01 3i  +0179  +0178  +0175 
+0171  +O183  +013 

(Average  Value  =  1.722  Volts.) 


+0177  +0180  +0183  +01 85  +0177  +0176  +0i82  +0185  +Oi85 
+0185  +0179  +0178  +O183  +O180  +0182  +0l85  +0182  +0179 
+0182  +0182  +0181  +0177  +0171  +0171  +0176  +0176  +0180 
+0176  +0180  +0101 

(Average  Value  =  1.772  Volts.) 

+0190  +0l85  +0186  +Ol84  +O186  +O188  +0194  +0191  +0195 

*0191  +0196  +0191  +0189  +0183  +0182  +0194  +0192  +0182 
+0183  +O188  +0185  +0193  +0192  +0190+0179  +0193  +0i9i 
+0190  +0190  +0002 

(Average  Value  =  1.8873  Volts.) 
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+0096  +0196  +0193  +0207  +0204  +01 89  +0195  +0196  +0193 
+0186  +0184  +0192  +0195  +0195  +0194  +0194  +0198  +0196 
+0195  +0191  +0l8l  +0194  +O189  +0202  +0193  +0193  +0193 
+0189  +0195  +0100 
(Average  Value  =  1.8727  Volts.) 

+0211  +0212  +0210  +0204  +0205  +0212  +0213  +0205  +0200 

+0195  +0202  +0200  +0197  +01 87  +0192  +0210  +0201  +0190 

+0199  +0201  +0193  +0197  +0203  +0204  +0208  +0201  +0219 

+0210  +0209  +0200 
(Average  Value  =  2.03  Volts.) 


NOTES  : 

START  OP  TRAVERSE  WATER  TEMP.  11  CERTIORATE. 

END  OF  TRAVERSE  WATER  TEMP.  12  CENTIGRADE. 

CHANNEL  SLOPE  0.3  inches. 

DALL  PLOW  TUBE  VERTICAL  MANOMETER  READING  7  Ft. 

WATER-ON-MERCURY. 
HEIGHT  OF  WATER  IN  CHANNEL  240nm. 

PROBE  RESISTANCE  READING  6.76  Ohms. 

PRECEDING  VOLTAGE  READINGS  TAKEN  AT,  0,1, 2, 3, 4, 
5,6,7,8,9,10,15,20,25,30,35,40,45,50,75,100,150, 
200  am  ah'ove  model  surface. 

PROBE  POSITION  -  CENTRE  LINE  OP  MODEL,  4  Ft.  FROM 
LEADING  EDGE. 
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+0042  +0071  +0062  +0065  +0058  +0079  +0042  +0054  +0047 

+0071  +0047  +004i  +0034  +0057  +0071  +0C54  +001 6  +0065 

+0076  +0034  +0075  +0068  +0053  +0052  +0077  +0072  +0049 
+0046  +0045  +0059  +000 

(Average  Value  =  0.5607  Volts.) 

+0075  +0083  +0070  +0113  +0091  +0075  +0094  +0121  +0096 

+0079  +0074  +0083  +0064  +0094  +0083  +0103  +0094  +0107 

+0106  +0087  +0084  +0092  +0074  +0092  +0122  +0111  +0095 
+0097  +0102  +00 

(Average  Value  =  0.9176  Volts.) 

+0130  +0104  +0116  +0084  +0075  +0102  +0082  +0131  +0121 

+0130  +0113  +0106  +0128  +0125  +0140  +0128  +0110  +0127 

+0118  +0135  +0094  +0101  +0133  +0100  +0094  +0123  +0116 
+0121  +0084  +0118 

(Average  Value  =  1.1297  Volts.) 

+0094  +0101  +0112  +0125  +0121  +0104  +0138  +0105  +0098 

+0120  +0129  +0125  +0124  +01 4l  +0115  +0139  +0088  +0086 

+0095  +0116  +0114  +0130  +0117  +0112  +0131  +01 48  +0120 

+0112  +0095  +0001 

(Average  Value  =  1.1569  Volts.) 

+0130  +01 44  +0146  +0129  +0093  +0139  +0127  +01 44  +0159 

+0114  +0114  +0103  +0116  +0152  +0129  +0110  +0120  +0130 

+0114  +0168  +0139  +0125  +0134  +0133  +01 65  +0122  +O139 

+01 44  +0115  +00 

(Average  Value  =  1.3093  Volts.) 

+0126  +0113  +0122  +0126  +0126  +0130  +0109  +0156  +0131 

+0113  +0116  +0112  +0151  +0128  +0136  +0106  +0116  +0125 

+0118  +0110  +0082  +0134  +0118  +0132  +0119  +0120  +0145 

+0132  +0107  +0000 

(Average  Value  =  1.2273  Volts.) 

+0001  +0107  +0126  +0121  +0151  +0149  +0119  +0129  +0150 

+0145  +0123  +0121  +0125  +0129  +0097  +0144  +0117  +0115 

+0134  +0137  +0111  +0149  +0151  +0147  +0121  +0131  +0l40 
+0130  +0120  +01 41 

(Average  Value  =  1.3035  Volts.) 


+01 4l  +0138  +0130  +0129  +0115  +0120  +0130  +0130  +0126 

+01 4l  +0143  +0128  +0121  +0135  +0129  +0117  +0127  +0125 
+0135  +0130  +0152  +01 44  +0127  +01 51  +0135  +0l40  +0150 
+0118  +0138  +001 

(Average  Value  =  1.3259  Volts.) 

+040  +01 48  +0124  +0137  +0129  +0172  +0159  +Ol64  +0147 
+0125  +0136  +0124  +0139  +0117  +0134  +0156  +0126  +0128 
+0144  +0145  +0139  +0128  +0152  +0162  +0115  +0122  +0152 
+0174  +0135  +0101 

(Average  Value  =  1.3910  Volts.) 

+01 4i  +0116  +0129  +0160  +01 48  +0166  +01 46  +0139  +0150 
+0165  +01 4i  +0154  +0132  +O161  +0136  +0147  +0131  +O178 
+0154  +0127  +0146  +01 49  +01 40  +0153  +0159  +01 43  +0149 
+0157  +0184  +0100 

(Average  Value  =  1.467  Volts.) 

+'0077  +0149  +0137  +0158  +0156  +0146  +0133  +01 46  +01 47 
+0170  +0194  +0175  +0177  +0166  +01 44  +0150  +0169  +0155 

+0147  +0170  +O186  +0175  +0143  +0134  +0169  +0169  +0172 

+01 68  +01 4i  +0100 

(Average  Value  =  1.541  Volts.) 

+003  +0172  +0146  +0176  +0169  +0168  +01 64  +0170  +0189 
+0185  +0186  +0187  +0156  +0160  +0151  +0157  +01 64  +0151 

+0155  +0181  +O161  +0170  +0l6l  +0176  +0156  +0l80  +0183 

+0182  +0160  +0001 

(Average  Value  =  1.6843  Volts.) 

+0190  +0183  +0171  +0180  +01 64  +0162  +0183  +0195  +0165 
+0190  +0167  +0l62  +0171  +0182  +0186  +0167  +0165  +0i65 

+0175  +0160  +0173  +0185  +0190  +0188  +0177  +O186  +0177 

+0180  +0179  +0167  +0172  +O161 

(Average  Value  =  1.7557  Volts.) 


+0197  +0196  +0203  +021 4  +0197  +0190  +0210  +0186  +0155 

+0176  +0l6l  +0174  +0192  +0193  +0189  +0194  +0183  +0189 

+01 64  +0170  +0199  +0211  +0203  +0202  +0196  +O188  +0194 

+0187  +0192  +0180  +0171  +0100 

(Average  Value  =  1.8890  Volts.) 
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PROBE  TRAVERSE  (2.3) 
NOTES  i 


START  OE  TRAVERSE  WATER  TEMP.  13  CENTIGRADE. 

END  OP  TRAVERSE  WATER  TEMP.  13*5  CENTIGRADE. 

CHANNEL  SLOPE  0.3  inches. 

DALL  PLOW  TUBE  VERTICAL  MANOMETER  READING  7  Ft. 

WA  TER-ON  -MERCURY . 

HEIGHT  OP  WATER  ABOVE  SURFACE  IN  CHANNEL  236  cm. 
PROBE  RESISTANCE  READING  6.76  Ohms. 

PRECEDING  VOLTAGE  READINGS  TAKEN  AT,  4, 5, 6, 7, 8, 

9,10,11,12,13,14,19,24,29  mm  above  surface  of  model. 
PROBE  DAMAGED  AT  34  mm. 

PROBE  POSITION  -  CENTRE  LINE  OF  MODEL,  6  Ft. 

PROM  LEADING  EDGE. 
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CALCULATED  VELOCITIES  FOR  PROBE  TRAVERSE  (1.1 


POSITION  F30K  LEADING  EDGE  0  Ft. 


POSITION  OF 
PROBE  ABOVE 
SURFACE 
(□a) 

AVERAGE 

VOLTAGE 

READING 

(Volts) 

CALCULATED 
VELOCITY 
(USING  EQ.4.6) 
(a/s) 

TRUE 

FLUID 

VELOCITY 

(a/s) 

0 

0.7741 

0.3725 

0.2332 

1 

0.8300 

0.3893 

0.2500 

2 

0.8072 

0.3825 

0.2432 

3 

0.8351 

0.3908 

0.2516 

4 

0.8430 

0.3932 

0.2540 

5 

0.8910 

0.4077 

0.2684 

6 

0.8400 

0.3923 

0.2531 

7 

0.9194 

0.4163 

0.2770 

8 

0.9271 

0.4186 

0.2793 

9 

0.9590 

0.4282 

0.2890 

10 

0.9646 

0.4299 

0.2906 

15 

0.9356 

0.4211 

0.2819 

20 

1.0351 

0.4511 

0.3118 

25 

1.1044 

0.4720 

0.3327 

30 

1.0572 

0.4578 

0.3185 

35 

1.0892 

0.4674 

0.3281 

40 

1.0697 

0.4615 

0.3223 

45 

1.0914 

0.4681 

0.3288 

50 

1.1026 

0.4714 

0.3322 

75 

1.0535 

0.4566 

0.3174 

100 

1.0621 

0.4593 

0.3200 

150 

1.1283 

0.4792 

0.3400 

200 

1.2566 

0.5178 

0.3786 

TABLE  Cl 


CALCULATED  VELOCITIES  FOR  PR 032  TRAVERSE  (2,1) 


POSITION  PP.OH  LEADING  EDGE  2  Ft. 


POSITION  0? 
PROEE  ABOVE 
SURFACE 
(mm) 

AVERAGE 

VOLTAGE 

READING 

(Volts) 

CALCULATED 
VELOCITY 
(USING  EQ.4.9) 
(o/s) 

TRUE 

FLUID 

VELOCITY 

(cj/s) 

0 

0.493 

0.4140 

0.1500 

1 

0.879 

0.5313 

0.2670 

2 

0.974 

0.5601 

0.2958 

3 

1.043 

0.5811 

0.3168 

4 

1.037 

0.5793 

0.3150 

5 

1.208 

0.6312 

0.3669 

6 

1.189 

0.6254 

0.3611 

7 

1.153 

0.6145 

0.3502 

8 

1.144 

0.6118 

0.3475 

9 

1.160 

0.6166 

0.3523 

10 

1.256 

0.6458 

0.3815 

15 

1.466 

0.7096 

0.4453 

20 

1.437 

0.7008 

0.4365 

25 

1.654 

0.7667 

0.5024 

30 

1.549 

0.7348 

0.4705 

35 

1.664 

0.7697 

0.5054 

40 

1.729 

0.7894 

0.5251 

45 

1.792 

0.8086 

0.5443 

50 

1.779 

0.8046 

0.5403 

75 

1.768 

0.8013 

0.5370 

100 

1.846 

0.8250 

0.5607 

150 

1.919 

0.8472 

0.5828 

200 

1.840 

0.8232 

0.5588 

TABLE  C2. 
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CALCULATED  VELOCITIES  FOR  PROBE  TRAVERSE  (2.2 


POSITION  FROM  LEADING  EDGE  4  Ft. 


POSITION  OF 
PROBE  ABOVE 
SURFACE 
(mo ) 


AVERAGE 

VOLTAGE 

READING 

(Volts) 


0.5890 
0.8850 
0.9480 
0.S990 
1.0560 
1.0650 
1.1970 
1.1804 
1.1935 
1.2482 
1.2004 
1.3960 
1.5133 
1.5683 
1.5669 
1.6824 
1.6350 
1.7015 
1.7220 
1.7220 
1.8873 
1.8727 
2.0300 


CALCULATED 
VELOCITY 
(USING  EQ.4.9) 
(m/s) 


0.4432 

0.5331 

0.5522 

0.5677 

0.5850 

0.5878 

0.6279 

0.6228 

0.6268 

0.6434 

0.6289 

0.6883 

0.7239 

0.7406 

0.7402 

0.7753 

0.7609 

0.7811 

0.7873 

0.7873 

0.8375 

0.8331 

0.8809 


TRUE 

FLUID 

VELOCITY 

(n/s) 


0.1789 

0.2688 

0.2879 

0.3034 

0.32C7 

0.3235 

0.3636 

0.3585 

0.3625 

0.3791 

0.3646 

0.4240 

0.4596 

0.4763 

0.4759 

0.5H0 

0.4966 

0.5168 

0.5230 

0.5230 

0.5732 

0.5688 

0.6165 


CALCULATED  VELOCITIES  FC-R  PROBE  TRAVERSE  (2.3) 
POSITION  FRON.  LEADING-  EDGE  6  Ft . 


POSITION  OF 
PROBE  ABOVE 
SURFACE 
(mm) 

AVERAGE 

VOLTAGE 

READING 

(Volts) 

CALCULATED 
VELOCITY 
(USING  EQ.4.9) 
(m/s) 

TRUE 

FLUID 

VELOCITY 

(m/s) 

4 

0.5607 

0.4346 

0.1703 

5 

0.9176 

0.5430 

0.2787 

6 

1.1297 

0.6074 

0.3431 

7 

1.1569 

0.6157 

0.3514 

8 

1.3093 

0.6620 

0.3977 

9 

1.2273 

0.6370 

0.3728 

10 

1.3035 

0.6602 

0.3959 

11 

1.3259 

0.6670 

0.4027 

12 

1.3910 

0.6868 

0.4225 

13 

1.4670 

0.7099 

0.4456 

14 

1.5410 

0.7324 

0.4680 

19 

1.6843 

0.7759 

0.5116 

24 

1.7557 

0.7976 

0.5332 

29 

1.8890 

0.8380 

0.5737 

TABLE  C4. 
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FbsrriONl  OF  PtoBE  SURFACE  1  Cmm'i 


•jV  . 


flpp  FBE3S?REA1£:  TELOCITY  :=:^O.L'56  5 : :  m/s.~ 
CTG:EjO?_'3CIBT:DA?.Y~^^ 

Z.  BOUNDARY  LAYER  TILOCIIY' =.:JD.  98  X  .ci.565-.  -- 


THICKNESS" OFTBCUNDARY  .LAYER 


e-i 


.PRESSTREAll  VELOCITY,  = 


'  G  .  565  g/a  . 


EDC-D  0?  BOUNDARY  LAYER  OCCURS: AT  98?-  PR 


BOUNDARY  rLAYER  VELOCITY ;  =:;p;.  98  ar ‘O'.  56 


PROPELLOR  TYPE  CURRENT  I*IETER 


CROSS  CHANNEL  VELOCITY  MEASUREMENTS 


— -  PLOW  DIRECTION 

POSITION  PROM  LEADING  EDGE  2  Ft. 

TOTAL  DISTANCE  ALONG  CHANNEL  11  Ft  6  ins. 


DALL  TUBE 

DEPTH  OF 

POSITION  OF 

PERIOD 

TOTAL 

REVS 

VERTICAL 

PROPELLOR 

PROPELLOR 

OF 

REVS/ 

PER 

MANOMETER 

ACROSS 

TEST 

TEST 

SEC 

READING 

(Ft) 

(mm) 

CHANNEL 

(mm) 

(sec) 

(n) 

VELOCITY 

USING 

EQUATION 


0.4324 

0.4426 

0.4443 

0.4460 

0.4443 

0.4443 

0.4426 

0.4426 

0.4168 


0.5243 

0.5481 

0.5651 

0.5651 

0.5685 

0.5651 

0.5668 

0.5651 

0.5651 

0.5464 

0.5243 


0.5789 

0.6043 

0.6196 

0.6179 

0.6264 

0.6230 

0.6264 

0.6162 


0.6060 


Ft" ‘WAT  jSH-ON  -MERCURY' . 


VELOCITIES  OBTAINED  FROM  PROFILE  FOR  PROBE  TRAVERSE 
(2.1)  PIG.  C3  APPENDIX  C. 


u  = 

0.565  m/s 

S  =  52  mm. 

DISTANCE 

PROF. 

SURFACE 

y 

(mm) 

VELOCITY 

PROM 

PROPILE 

u 

(m/s) 

u 

U 

(y  A)% 

-!♦ 

— o - 

“07330 

0.266 

-00 

0.000 

_ 

1 

0.258 

0.457 

0.000 

0.454 

— 

2 

0.294 

0.520 

0.693 

0.521 

— 

3 

0,317 

0.561 

1.099  ' 

0.565 

— 

4 

0.335 

0.593 

1.386 

0.598 

— 

5 

0.350 

0.619 

1.609 

0.626 

— 

6 

0.362 

0.641 

1.792 

0.679 

— 

7 

0,375 

0.664 

1.946 

0.669 

— 

8 

0.385 

0.681 

2.079 

0.687 

— 

9 

0.395 

0.699 

2.197 

— 

10 

0.404 

0.715 

2.303 

0.719 

— 

15 

0.438 

0.775 

2.708 

0.779 

— 

20 

0.465 

0.823 

2.996 

0.826 

— 

25 

0.486 

0,860 

3.219 

0.863 

— 

30 

0.504 

0.892 

3.401 

0.895 

— 

35 

0.519 

0.91S 

3.555 

0.923 

— 

40 

0.533 

0.943 

3.689 

0.948 

_ 

45 

0.543 

0.961 

3.807 

0.971 

0.964 

50 

0.550 

0.973 

3.912 

0.992 

0.990 

52 

0.554 

0.980 

3.951 

1.000 

— 

TABLE  El 


r 


i 


VELOCITIES  OBTAINED  FROM  PROFILE  FOR  PROBE  TRAVERSE 
(2.2)  FIG.  C5  APPENDIX  C. 


u 

=  0.565  m/s 

S 

=  65  mm. 

DISTANCE 

FROM 

SURFACE 

y 

(mm) 

VELOCITY 

FROM 

FROFILE 

u 

(m/s) 

TJ 

iogt(y) 

(y/4? 

(y/S?* 

mmsmm 

0.179 

0.316 

-CO 

0.000 

- 

— 

i 

0.268 

0.474 

0.000 

0.434 

0.474 

— 

2 

0.286 

0.506 

0.693 

0.458 

- 

— 

3 

0.304 

0.538 

1.099 

0.540 

- 

— 

4 

0.321 

0.568 

1.386 

0.573 

— 

5 

0.335 

0.593 

1.609 

0.598 

— 

— 

6 

0.349 

0.618 

1.792 

0.620 

- 

— 

7 

0.358 

0.634 

1.946 

0.640 

— 

— 

8 

0.369 

0.653 

2.079 

0.657 

- 

9 

0.380 

0.673 

2.197 

0.673 

- 

— 

10 

0.387 

0.685 

2.303 

0.687 

— 

— 

15 

0.425 

0.752 

2.708 

0.746 

- 

- 

20 

0.454 

0.804 

2.996 

0.790 

0.810 

— 

25 

0.476 

0.842 

3.219 

0.826 

0.843 

_ 

30 

0.492 

0.870 

3.401 

0.857 

0.871 

— 

35 

0.506 

0.896 

3.555 

0.864 

0.896 

— 

40 

0.517 

0.915 

3.689 

0.908 

0.916 

— 

45 

0.526 

0.930 

3-807 

0-929 

— 

— 

50 

0.535 

0.947 

3.912 

0-948 

- 

— 

55 

0.540 

0.956 

4.007 

0.967 

- 

0.959 

60 

0.548 

0.970 

4.094 

0.984 

- 

0.980 

65 

0.554 

0.980 

4.174 

1.000 

— 

— 

•TABLE  E2. 
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PROBE  TRAVERSE  (2.1] 


INKER  LAV/  REGION  2.3  and  4  bb. 


■'  i 

>*m0* 

Dist¬ 

ance 

y 

(mm) 

_ 

u 

(m/s) 

u 

u. 

(yu./^) 

log  (yu.A>) 

n 

0.294 

14.700 

30.769 

3*426 

0.020 

0.317 

15*350 

46.153 

3*832 

B 

0.335 

16.750 

61.538 

4.129 

mm 

0.294 

14.000 

32.308 

3*475 

0.021 

mm 

0.317 

15*095 

48.462 

3*881 

0.335 

15*952 

64*615 

4.169 

mm 

0.294 

13*364 

33*846 

3*522 

0.022 

mm 

0.317 

14.409 

50.769 

3*927 

■■ 

0.335 

15*227 

67*692 

4.215 

m 

0.294 

12.783 

35*385 

3.566 

0.023 

Mm 

0.317 

13*783 

53*077 

3.972 

KK 

0.335 

14*565 

70.769 

4*259 

mm 

0.294 

12.250 

36.923 

3*609 

0.024 

if 

.0.317 

13*208 

55*385 

4*014 

KB 

0.335 

13*958 

73*646 

4*302 

m 

0.294 

11.760 

38.462 

3*649 

0.025 

Mm. 

0.317 

12.680 

57*692 

4.055 

HH 

0.335 

13*400 

76*923 

4*343 

mm 

0.294 

11.308 

40.000 

3.689 

0.026 

0.317 

12.192 

60.000 

4.094 

B 

0.335 

12.885 

80.000 

4.382 

TABLE  E3 


PROBE  TRAVERSE  (2.2) 


INKER  LAV/  RES  ION  3,4  and  5  mm. 


(m/s) 

Dist¬ 

ance 

y 

( mm) 

u 

(m/s) 

u 

u«, 

(yUeA>) 

log  (yu«4>) 

0.304 

14.476 

48.462 

3.881 

0.021 

0.321 

15.266 

64.616 

4.168 

0.335 

15.952 

80.769 

4.392 

0.304 

13.818 

50.769 

3.927 

0.022 

0.321 

14.591 

67.692 

4.215 

0.335 

15.227 

84.615 

4.438 

0.304 

13.217 

53.077 

3.972 

0.023 

0.321 

13.957 

70.769 

4.259 

0.335 

14.565 

88.462 

4.483 

0.304 

12.667 

55.385 

4.014 

0.024 

0.321 

13.375 

73.846 

4.302 

0.335 

13-958 

92.308 

4.525 

0.304 

12.160 

57.692 

4.055 

0.025 

0.321 

12.840 

76.923 

4.343 

0.335 

13-400 

96.154 

4.566 

0.304 

11.692 

60.000 

4.094 

0.026 

0.321 

12.346 

80.000 

4.382 

0.335 

12.885 

100.000 

4.605 

3 

0.304 

11.259 

62.308 

0.027 

4 

0.321 

11.809 

83*077 

5 

0.335 

12.407 

103.846 

TABLE  E4. 
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PROBE  TRAVERSE  (2.1i 


Dist¬ 

ance 

(nm) 

u0  =  0 

.022 

u*  =  0 

.024 

u.  =  0 

.026 

u 

u„ 

'•8.0?) 

u 

Uo 

u 

Uo 

1 

11.727 

2.829 

10.750 

2.916 

9.923 

2.996 

2 

13.364 

3.522 

12.250 

3. 609 

11.308 

3.689 

3 

14.409 

3.927 

13.208 

4.014 

12.192 

4.094 

4 

15.227 

4.215 

13.958 

4.302 

12.885 

4.382 

5 

15.909 

4.438 

14.583 

4.525 

13.462 

4.605 

6 

16.455 

4.621 

15.083 

4.708 

13.923 

4.787 

7 

17.045 

4.775 

15.625 

4.862 

14.423 

4.942 

8 

17.500 

4.908 

16.042 

4.995 

14.808 

5. 075 

9 

17.955 

5.026 

16.458 

5.113 

15.192 

5.193 

10 

18.364 

5.131 

16.833 

5.218 

15-538 

5.298 

15  • 

19-909 

5.537 

18.250 

5.624 

16.846 

5.704 

20 

21.136 

5.825 

19.375 

5.911 

17.885 

5.991 

25 

22.091 

6  •  048 

20.250 

6.135 

18.692 

6.215 

30 

22.909 

6.230 

21.000 

6.317 

19.385 

6.397 

35 

23-591 

6.384 

21.625 

6.471 

19.962 

6.551 

40 

24.227 

6.518 

22.208 

6.605 

20.500 

6.685 

45 

24.682 

6.635 

22.625 

6.722 

20.885 

6.803 

50 

25.000 

6.741 

22.917 

6.828 

21.154 

6.908 

TABLE  E5. 
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PROBE  TRAVERSE 


(2.2). 


Dist¬ 

ance 

y 

(ora) 

Uo  =  0 

.022 

u.=  0 

.024 

Uo  =  0 

.026 

u 

u. 

u 

Uo 

Ml") 

u 

U. 

Mt) 

1 

12.182 

2.829 

11.167 

2. 916 

10.308 

2.996 

2 

13.000 

3.522 

11.917 

3.609 

11.000 

3.689 

3 

13.818 

3.927 

12.667 

4.014 

11.692 

4.C94 

4 

14.591 

4.215 

13.375 

4-302 

12.346 

4.382 

5 

15-228 

4.438 

13.959 

4-525 

12.885 

4. 605 

6 

15.864 

4.621 

14.542 

4. 70S 

13.423 

4.788 

7 

16.273 

4.775 

14.917 

4.862 

13.769 

4.942 

8 

16.773 

4. 908 

15.375 

4.995 

14.192 

5.075 

9 

17.273 

5.026 

15.833 

5.113 

14.616 

5.193 

10 

17.591 

5.131 

16.125 

5.218 

14.885 

5.298 

15 

19.318 

5.537 

17.708 

5.624 

16.346 

5.704 

20 

20.636 

5.825 

18.917 

5.911 

17.462 

5.991 

25 

21.636 

6.048 

19.833 

6.135 

18.308 

6.215 

30 

22.364 

6.230 

20.500 

6.317 

18.923 

6.397 

35 

23.000 

6.384 

21.083 

6.471 

19-462 

6.551 

40 

23.500 

6.518 

21.542 

6.605 

19.885 

6.685 

45 

23-909 

6.635 

21.917 

6.722 

20.231 

6.802 

50 

24.318 

6.741 

22.292 

6.828 

20.577 

6. 908 

55 

24.546 

6.836 

22.500 

6.923 

20.769 

7.C03 

60 

24.909 

6.923 

22.833 

7.010 

21.077 

7.090 

TABLE  E6. 
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